Bulletin of Engineering Geology and the Environment最新文献

{"title":"Evolution characteristics of unfrozen water content and damage for saturated sandstone during freezing–thawing cycle process based on the in-situ nuclear magnetic resonance","authors":"Shuailong Lian,&nbsp;Lianyang Zhang,&nbsp;Yanlin Zhao,&nbsp;Qiuhong Wu,&nbsp;Can Du,&nbsp;Wen Wan","doi":"10.1007/s10064-025-04330-y","DOIUrl":"10.1007/s10064-025-04330-y","url":null,"abstract":"<div><p>The accurate determination of unfrozen water content is significant to evaluating the freezing process of pore water and to revealing the freezing damage mechanism. In this study, an in-situ low filed nuclear magnetic resonance (LF-NMR) testing system equipped with a low temperature thermal cycle system were designed to conduct freezing–thawing cycle tests. The pore water freezing process of three sandstone samples (Coarse-grained sandstone, Medium-grained sandstone, Fine-grained sandstone) was studied by using the LF-NMR technology. The results indicate that the free water content of the three sandstone samples decreases sharply as the temperature decrease, the free water signal intensity and peak region decrease significantly from 0℃ to -5 ℃, whereas bound water that exists in small pores needs a much lower temperatures to freezing. The larger the equivalent average pore size, the faster the water–ice conversion rate. The more movable water the sample contains, the greater the frost heaving force generated inside the sample. At the temperature range of 0 ~ -10 ℃, the freezing-heaving strain increasing sharply, and gradually become gentle as the temperature declining. Coarse-grained sandstone with the highest porosity has the largest freezing-heaving strain, while the Fine-grained sandstone with the lowest porosity has the smallest freezing-heaving strain. The freezing damage of Coarse-grained sandstone is more serious than the other two sandstone. P-wave velocities and scanning electron microscope (SEM) test also indicate that Coarse-grained sandstone has more serious damage than Fine-grained sandstone subjected to freezing–thawing cycle.</p></div>","PeriodicalId":500,"journal":{"name":"Bulletin of Engineering Geology and the Environment","volume":"84 6","pages":""},"PeriodicalIF":3.7,"publicationDate":"2025-05-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143938620","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Evolution law of phreatic water leakage caused by underground coal mining in the semiarid region of Northwest China","authors":"Shiliang Liu,&nbsp;Ao Wang,&nbsp;Wenhui Zhang,&nbsp;Dong Zhou,&nbsp;Yaowen Wu,&nbsp;Shanlin Wang,&nbsp;Yusheng Zheng,&nbsp;Deqiang Mao","doi":"10.1007/s10064-025-04312-0","DOIUrl":"10.1007/s10064-025-04312-0","url":null,"abstract":"<div><p>Currently, studies on the permeability evolution characteristics of overlying aquiclude protective layers caused by coal mining focus on single lithological protective layers and assume the permeability coefficient remains constant. However, these studies fail to consider the variation characteristics of the combination protective layer structure and permeability coefficient. Therefore, an analytical method is proposed to study coal seam leakage under mining conditions in the blown-sand beach region based on the structure and permeability coefficient of the combination protective layer. First, the stress path of the overlying combination aquiclude under coal mining disturbance is comprehensively considered. Based on this, triaxial loading and unloading seepage creep experiments are conducted with different proportions of overlying combination aquiclude. The analytical relationship between the permeability coefficient of the samples and loess proportion, stress level, and soil depth in the stress recovery stage is determined, leading to the establishment of a creep permeability coefficient evolution model for the overlying combination aquiclude of the coal seam under the stress path of coal mining. Second, the creep permeability coefficient evolution model is integrated with a fusion algorithm in COMSOL numerical simulation software. Numerical simulations are then performed to examine the evolution law of phreatic leakage during coal seam mining and recovery, revealing a relationship curve in which leakage gradually decreases over time before stabilizing in the post-mining recovery stage. Finally, based on mathematical and statistical methods, a phreatic leakage evolution model is developed for both mining and post-mining stages to provide a theoretical basis for environmental protection.</p></div>","PeriodicalId":500,"journal":{"name":"Bulletin of Engineering Geology and the Environment","volume":"84 6","pages":""},"PeriodicalIF":3.7,"publicationDate":"2025-05-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143944358","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Energy release and related sensitivity analysis of anisotropic CJBs under compression conditions without and with lateral pressure","authors":"Bin Gong,&nbsp;Yongyi Wang,&nbsp;Yongjun Zhang,&nbsp;Xiaoshuang Li","doi":"10.1007/s10064-025-04317-9","DOIUrl":"10.1007/s10064-025-04317-9","url":null,"abstract":"<div><p>The cumulative energy released by micro-cracks when columnar jointed basalts (CJBs) attain peak strength (micro-crack energy index (MCEI)) under compression is a crucial foundation for understanding the mechanical behavior of CJBs and informing engineering strategies for monitoring and reinforcement. Beginning with the construction of detailed CJB images owning varied settings, the meso-scale damage mechanics, stochastic strength principle, and continuum mechanics are integrated. Employing the rock failure process analysis method improved by digital image correlation, the related heterogeneity numerical models are established using the generated images, facilitating a comprehensive scrutiny of the force–deformation features, fracture process, and energy change of CJBs exposed to compression conditions without and with lateral pressure (CCWOLP and CCWLP). Then, the effects of different factors on the micro-crack quantity index (MCQI) and MCEI are investigated, and the sensitivity analysis is conducted to clarify the impact relationships and recognize pivotal variables. The results show that under the CCWOLP and CCWLP, the increased rock homogeneity and larger column diameter lead to the intensified stress concentration at joints and accelerate the occurrence of the MCEI. Conversely, the higher residual strength coefficient of joints can delay the MCEI. The specimen, with the greater joint elastic modulus / no secondary joint set, may release MCEI later, and the MCEI value may be lower, hinging on the specific column dip angle. Under the CCWOLP, the dominant sensitive variables for the MCEI comprise joint constitutive correlation, column diameter, column irregularity degree. However, under the CCWLP, the primary sensitive variables regarding the MCEI contain rock mass restraint state, column irregularity degree, joint mechanical property. Additionally, the corresponding fitting models provide theoretical underpinnings for guiding the monitoring, enhancement, functioning, and servicing of CJB project undertakings.</p></div>","PeriodicalId":500,"journal":{"name":"Bulletin of Engineering Geology and the Environment","volume":"84 6","pages":""},"PeriodicalIF":3.7,"publicationDate":"2025-05-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143938644","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"The longitudinal loading mode evaluation of fault-crossing tunnel by a new load equivalence method based on a Pasternak model","authors":"Tianqiang Wang,&nbsp;Ping Geng,&nbsp;Chenyang Xiang,&nbsp;Guoguo Liu,&nbsp;Qi Wang,&nbsp;Lin Deng,&nbsp;Shiqiang Ma","doi":"10.1007/s10064-025-04273-4","DOIUrl":"10.1007/s10064-025-04273-4","url":null,"abstract":"<div><p>The deformation and failure of fault-crossing tunnel attributes to its weak mechanical behavior near the fault core under a certain longitudinal loading mode when a tunnel is passively subjected to a normal faulting. A reasonable analytical solution to deal with the complex dislocation issues for fault-crossing tunnel by an elastic beam model is convenient to reveal the dislocation response of tunnel. In view of this, the longitudinal loading mode of tunnel was discussed by a new load equivalence method based on a hypothetically longitudinal deformation pattern of tunnel along the transition zone. Subsequently, a numerical model was built to reveal the mechanism response of tunnel emerged in the strata containing a fault core to verify the correctness of the proposed analytical evaluation method. And then, the proposed load equivalence method was applied to the scenarios involving fault fracture zone and its adaptability was confirmed numerically and experimentally. The results showed that the presumptive loading mode under the proposed load equivalence method of tunnel was available relying on the deformation and force indices. Under the deformation profile mode at transition zone, the maximum bending moment appears at the ends of the deformed profiles, and the sections near fault core present the negative shear behavior. The analytical prediction at loading mode can be widely applied at the fault-crossing tunnel emerged in the field with fault fracture zone, and the length of the transition zone of tunnel is about 1.2 times of the width of fault fracture zone appropriately.</p></div>","PeriodicalId":500,"journal":{"name":"Bulletin of Engineering Geology and the Environment","volume":"84 6","pages":""},"PeriodicalIF":3.7,"publicationDate":"2025-05-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143938645","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Shear behavior of gravel-block soil of the Qinghai-Tibet Plateau based on large-scale direct shear test and numerical simulation","authors":"Tao Wen,&nbsp;Wenjun Jia,&nbsp;Zhi Quan,&nbsp;Wei Guo,&nbsp;Yankun Wang,&nbsp;Ningsheng Chen","doi":"10.1007/s10064-025-04346-4","DOIUrl":"10.1007/s10064-025-04346-4","url":null,"abstract":"<div><p>The shear behavior of gravel-block soil (GBS) is unique and significant for evaluation the failure mechanism of GBS landslide on the Qinghai-Tibet Plateau. This study focuses on interpreting the shear behavior observed in the GBS during large-scale direct shear tests conducted on a landslide in Jiacha County, Tibet, China. The tests considered coarse particle content (CPC), dry density, and moisture conditions. Additionally, a discrete element numerical model, scaled to match the laboratory testing dimensions, was developed to simulate the large-scale direct shear tests on GBS. Results indicated that an increase in CPC improves the strength of the GBS, as it enhances the framework strength through interlocking between gravel blocks and between gravel blocks and the soil mass. The critical CPC for shear failure of the GBS exhibits a decreasing trend as the dry density increases. Furthermore, particle crushing rate (PCR) of the GBS is positively correlated with CPC, vertical pressure, and dry density. The simulation results show good agreement with the test results, providing insights into the damage-shear fracture mechanism of typical GBS under large-scale direct shear tests. The research outcomes provide a theoretical basis for the prevention and control of geological hazards in the Qinghai-Tibet Plateau.</p></div>","PeriodicalId":500,"journal":{"name":"Bulletin of Engineering Geology and the Environment","volume":"84 6","pages":""},"PeriodicalIF":3.7,"publicationDate":"2025-05-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143938556","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Frictional characteristics of calcareous sand subjected to oscillatory normal stress","authors":"Wengang Dang,&nbsp;Zhenchao Bai,&nbsp;Linchong Huang,&nbsp;Kang Tao,&nbsp;Xiang Li,&nbsp;Jianjun Ma,&nbsp;Yu Liang","doi":"10.1007/s10064-025-04334-8","DOIUrl":"10.1007/s10064-025-04334-8","url":null,"abstract":"<div><p>Calcareous sand is commonly utilized in coastal and marine constructions, and dynamic stress disturbances caused by the natural environment or human activities make calcareous sand to exhibit different shear properties. Hence, it is important to understand calcareous sand’s frictional behavior under variable normal stress. A sequence of direct shear tests on calcareous sand sourced from South China Sea were conducted using a dynamic shear box device. Both constant and oscillatory normal stresses were considered. The frictional characteristics under oscillatory normal stress were investigated by applying a constant shear rate (1 ~ 20 mm/min) and a sinusoidal normal load with different normalized amplitudes (0% ~ 80%) and frequencies (0.1 ~ 2 Hz). Experimental findings indicate that the shear stress, apparent friction coefficient, and normal displacement varied periodically in response to oscillatory normal stress, with phase lags observed among the normal stress, shear stress, and apparent friction coefficient. Unlike constant normal stress, oscillatory stress can either increase or decrease the shear strength, and the frictional strengthening or weakening is dictated by critical rate–amplitude–frequency values. Additionally, the shear stress results are characterized using the rate and state friction law. These findings shed light on the frictional behavior of calcareous sand exposed to complex stress states and beneficial for the mitigation of relative geo-risks.</p></div>","PeriodicalId":500,"journal":{"name":"Bulletin of Engineering Geology and the Environment","volume":"84 6","pages":""},"PeriodicalIF":3.7,"publicationDate":"2025-05-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143930138","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Hybrid prediction for reservoir landslide deformation based on multi-source InSAR and deep learning","authors":"Qiuyu Ruan,&nbsp;Fasheng Miao,&nbsp;Yiping Wu,&nbsp;Beibei Yang,&nbsp;Fancheng Zhao,&nbsp;Weiwei Zhan","doi":"10.1007/s10064-025-04345-5","DOIUrl":"10.1007/s10064-025-04345-5","url":null,"abstract":"<div><p>Time series Interferometric Synthetic Aperture Radar (InSAR) technology has been proven to be an effective tool for measuring landslide movements. However, previous research has primarily focused on the innovation and application of InSAR technology, its exploration in the analysis and prediction of slope displacement remains to be explored. Analyzing the coupling relationship between InSAR derived displacement and triggering factors, and applying these into landslide displacement prediction, can provide valuable insights for landslide disaster prevention and mitigation early warning systems. In this study, multi-source InSAR data were collected to obtain the displacement of the Shuping landslide in the Three Gorges Reservoir area. We characterized the temporal and spatial displacement of the Shuping landslide and discussed the response mechanism between landslide movement and triggering factors. Subsequently, the landslide displacement was decomposed into trend and periodic term by the wavelet analysis (WA) algorithm. Long short-term memory (LSTM) and Bidirectional-LSTM (Bi-LSTM) algorithm were employed for time series prediction modeling, and parameter optimization was conducted using the grey wolf optimization (GWO) algorithm. Finally, combining InSAR data with displacement prediction models, we conducted InSAR-assisted displacement prediction research and confirmed its effectiveness in improving prediction accuracy. The findings demonstrate the feasibility of applying InSAR technology in landslide displacement prediction, offering a reference for the prediction and prevention of reservoir-induced landslides.</p></div>","PeriodicalId":500,"journal":{"name":"Bulletin of Engineering Geology and the Environment","volume":"84 6","pages":""},"PeriodicalIF":3.7,"publicationDate":"2025-05-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143930139","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Paleoliquefaction mapping for Bangladesh from the records of historical earthquakes","authors":"Tasnim Jabin Jui,&nbsp;Md. Shakhawat Hossain,&nbsp;A. S. M. Maksud Kamal,&nbsp;Md. Zillur Rahman,&nbsp;Masaba Tasneem Sahebi,&nbsp;Abrar Hossain","doi":"10.1007/s10064-025-04316-w","DOIUrl":"10.1007/s10064-025-04316-w","url":null,"abstract":"<div><p>Paleoliquefaction investigations are crucial for assessing seismic hazard potential and identifying regions susceptible to liquefaction, which is essential for seismic risk-sensitive land-use planning. This research aimed to identify paleoliquefaction sites by reviewing documented descriptions of the damages and ground deformations in Bangladesh during three significant historical earthquakes: the Bengal Earthquake (1885), the Great Assam Earthquake (1897), and the Srimangal Earthquake (1918). A paleoliquefaction map for Bangladesh was generated, locating the paleoliquefaction sites during these three major historical earthquakes. In addition, Standard Penetration Test (SPT) blow count and Down-hole Seismic Tests (DST) were conducted at selected locations to assess the Liquefaction Potential Index (LPI) by using deterministic (simplified) and probabilistic procedures. The results confirmed a high likelihood of liquefaction during future large-magnitude earthquakes. The research outcome will help to distinguish and characterize Bangladesh’s susceptible regions to soil liquefaction during potential earthquakes in the future and is recommended for consideration in large-scale construction or development plans.</p></div>","PeriodicalId":500,"journal":{"name":"Bulletin of Engineering Geology and the Environment","volume":"84 6","pages":""},"PeriodicalIF":3.7,"publicationDate":"2025-05-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143930140","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Stress–strain response and deformation control method of a roof roadway in transversely isotropic mudstone-clay strata","authors":"Huaidong Liu,&nbsp;Changyou Liu,&nbsp;Jingxuan Yang,&nbsp;Jun Wang","doi":"10.1007/s10064-025-04287-y","DOIUrl":"10.1007/s10064-025-04287-y","url":null,"abstract":"<div><p>The mudstone-clay composite roof roadway exhibits distinct transversely isotropic characteristics, rendering the prediction of uncoordinated deformation in surrounding rock complex and challenging. Based on transversely isotropic theory, the deformation parameters of mudstone-clay composite in different directions are calculated by true triaxial experiment, and the elastic modulus is determined as the key parameter affecting the uncoordinated deformation of composite roof roadway. The stress and strain expression of roadway surrounding rock is theoretically deduced, and the utilization of ultra-high strength bolts is proposed to control the uncoordinated roadway deformation. The results indicate significant variations in the elastic modulus of the mudstone-clay assemblage in both horizontal and vertical directions, with Poisson’s ratio showing a narrow range of variation. The composite with saturated clay exhibits reduced deformation resistance and more pronounced transverse isotropy compared to the composite with dry clay. The stress concentration is highest near the inflection point of the roadway. The roadway ribs experience vertical stress increase and horizontal stress decrease, while the roof and floor strata mainly undergo vertical stress decrease and horizontal stress increase. The strain in roadway surrounding rock mainly shows vertical strain, especially with the roadway roof exhibiting the highest vertical strain peak and the largest influence range. The vertical strain of roadway surrounding rock can be significantly reduced by increasing the value of its vertical elastic modulus <i>E</i>₂, if it is less than 0.15GPa. However, value higher than this has little effect on the strain. According to the field observation, the utilization of ultra-high strength bolt support (<i>E</i>₂ &gt; 0.15GPa) in comparison to Q₂₃₅ threaded steel resin bolt support (<i>E</i>₂ &lt; 0.15GPa) demonstrates a significant reduction in roadway uncoordinated deformation, thereby validating the accuracy of theoretical research.</p></div>","PeriodicalId":500,"journal":{"name":"Bulletin of Engineering Geology and the Environment","volume":"84 6","pages":""},"PeriodicalIF":3.7,"publicationDate":"2025-05-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143929916","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Analytical solution for calculating the water inflow in composite lined tunnels with waterproof and drainage systems","authors":"Jun Yu,&nbsp;Chi Zhang,&nbsp;Dongkai Li","doi":"10.1007/s10064-025-04347-3","DOIUrl":"10.1007/s10064-025-04347-3","url":null,"abstract":"<div><p>In the design of waterproof and drainage systems for composite lined tunnels, water inflow is an important parameter. As the main method for predicting tunnel inflow at present, the existing analytical solutions tend to ignore the influence of seepage in the surrounding rock around the tunnel and the two-dimensional seepage characteristics inside the tunnel structure. This study addresses these two issues and derives an analytical solution for the water inflow of a composite lined tunnel with waterproof drainage system using the conformal transformation and the sub-regional separation of variables method. Comparison with numerical simulations and existing solutions shows that the analytical solution has good accuracy under any geotextile permeability coefficient<i> k</i>_<i>g</i>, soil permeability coefficient and primary lining permeability coefficient, which is obviously superior to the existing solutions. Finally, parameter analyses based on actual projects are carried out to investigate the relationship between geotextile parameters, drainage pipe spacing, and tunnel groundwater discharge. It is found that as <i>k</i>_<i>g</i> decreases, the level of groundwater inflow into the tunnel decreases. For a fixed <i>k</i>_<i>g</i>, smaller drainage pipe spacing and larger geotextile thickness increase groundwater inflow. In addition, when the surrounding rock and lining parameters are fixed, it is feasible to adjust the geotextile parameters and drainage pipe spacing to regulate the tunnel groundwater discharge, and this study can provide a reference for the design of the waterproof and drainage systems of composite lined tunnels.</p></div>","PeriodicalId":500,"journal":{"name":"Bulletin of Engineering Geology and the Environment","volume":"84 6","pages":""},"PeriodicalIF":3.7,"publicationDate":"2025-05-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143929917","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}