Geophysical Prospecting最新文献

{"title":"Seismic crosshole full-waveform inversion of high-frequency SV-waves for glacial sediment characterization","authors":"Sarah Beraus,&nbsp;Daniel Köhn,&nbsp;Thomas Bohlen,&nbsp;Thomas Burschil,&nbsp;Bennet Schuster,&nbsp;Hermann Buness,&nbsp;Gerald Gabriel","doi":"10.1111/1365-2478.70024","DOIUrl":"https://doi.org/10.1111/1365-2478.70024","url":null,"abstract":"<p>We present a two-dimensional, high-resolution full-waveform inversion approach of crosshole seismic data to image fine glacial sediments in an overdeepened Alpine valley, the Tannwald Basin (ICDP site 5068_1) north of Lake Constance. A vertically polarizing shear-wave source is employed to excite elastic waves, which are recorded by an eight-station three-component geophone string at 105–134 m depth. Based on traveltime tomography models derived separately from the P- and S-wave first arrival picks, we reconstruct the vertically polarized shear-wave velocity distributions from the vertical component data using isotropic, elastic full-waveform inversion. In order to mitigate the effects of source and receiver coupling, we use the global correlation norm as an objective function. The resolution is enhanced by applying an anisotropic gradient filter. The final model shows small-scale structures that look similar to a structural image obtained by seismic imaging at frequencies above 200 Hz. It provides a satisfying data fit and resolves layers of at least 1 m thickness that correlate well with the lithology derived from the core information within the errors of the source and receiver positions.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 5","pages":"1587-1605"},"PeriodicalIF":1.8,"publicationDate":"2025-04-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.70024","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143949731","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Seismic full waveform inversion for fracture parameters in anisotropic media","authors":"Ujjwal Shekhar,&nbsp;Morten Jakobsen,&nbsp;Ivan Pšenčík,&nbsp;Kui Xiang","doi":"10.1111/1365-2478.70025","DOIUrl":"https://doi.org/10.1111/1365-2478.70025","url":null,"abstract":"<div>\u0000 \u0000 <p>Vertical fractures are often reported in sedimentary rocks. The detection of these inherent fractures is important before carrying out the carbon dioxide sequestration in these rocks. The detection of the fractures is also crucial for an accurate estimation of the moment tensor from microseismic waveform data. In this study, we use the distorted Born iterative method to perform seismic full waveform inversion for the parameters of vertical fractures in sedimentary formations. The distorted Born iterative method is based on transforming a nonlinear inverse scattering problem into a series of linear inverse problems by using the distorted Born approximation. We work in the frequency domain and use a volume integral equation method to solve the direct scattering problem. A heterogeneous, generally anisotropic medium, is iteratively updated using the matrix-free formulation of Fréchet derivatives and their adjoint. In the distorted Born iterative method, the heterogeneous medium Green's function is also updated after each iteration, which is not done in the classical Born iterative method. In our implementation, we assume that the fractures are thin, vertical and parallel to each other. The background, in which fractures are embedded, is transversely isotropic with a vertical axis of symmetry. The vertically transversely isotropic background can be inhomogeneous. In an isotropic background, it is common to invert for a single tangential fracture weakness along with a normal fracture weakness. However, in a vertically transversely isotropic background, the horizontal-tangential and the vertical-tangential fracture weaknesses vary, and therefore we invert for three fracture weaknesses. In numerical experiments, we employ a cross-hole seismic configuration and invert synthetic waveform data for fracture weaknesses. The radiation pattern analysis is performed to investigate the cross-talk among different fracture weaknesses. It is found that the horizontal-tangential fracture weakness is better resolved than the other fracture weaknesses, which is confirmed through the numerical results.</p></div>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 5","pages":"1606-1634"},"PeriodicalIF":1.8,"publicationDate":"2025-04-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143950233","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Horizontal loop electromagnetic data inversion through novel metaheuristics and its applications to ore deposit explorations","authors":"Hanbing Ai,&nbsp;Arkoprovo Biswas,&nbsp;Yunus Levent Ekinci,&nbsp;Çağlayan Balkaya,&nbsp;Ahmad Alvandi,&nbsp;Yue Shen,&nbsp;Zhuofan Liu,&nbsp;Arka Roy,&nbsp;Kejia Su","doi":"10.1111/1365-2478.70027","DOIUrl":"https://doi.org/10.1111/1365-2478.70027","url":null,"abstract":"<p>The success of horizontal loop electromagnetic method in determining ore body model parameters is intricately tied to the inversion process. However, geophysical inverse problems often grapple with ill-posed and non-unique mathematical difficulties. To handle this situation, we presented two new metaheuristics: the Barnacles Mating Optimizer and the Hunger Games Search algorithm. This pioneering approach is the first in the literature to employ both optimizers simultaneously to invert electromagnetic anomalies. To test the core functionality of these inversion tools, we simplified the horizontal loop electromagnetic problem by using a single frequency. Before carrying out the data inversions, a synthetic anomaly caused by a simple multi-source model was used to perform some modal analyses. These analyses revealed the composite modality character of the inversion process and highlighted the necessity of an algorithm that effectively balances global exploration (thoroughly searching the model space while avoiding local minima) and local exploitation, which focuses on approximating the global optimum closely. In addition, the findings indicated the importance of conducting post-inversion analyses to evaluate potential uncertainties in the estimations. Initial experiments with nine challenging benchmark functions revealed that Hunger Games Search outperforms in tackling optimization problems. The performances of both algorithms were then examined using synthetic horizontal loop electromagnetic anomalies and two field data sets from Australia. To achieve a balanced comparison, we fine-tuned the control parameters of both algorithms, optimizing their performance under comparable conditions. These assessments provided valuable insights into the effectiveness of each algorithm when properly configured. Uncertainty appraisal analyses mostly confirmed the consistency and mathematical consistency of the obtained solutions. In the real data inversion cases, the outputs of Hunger Games Search showed strong agreement with the drilling data. This study illustrates the stronger performance of the Hunger Games Search algorithm relative to Barnacles Mating Optimizer in the presented scenarios and demonstrates its effectiveness in estimating source parameters using horizontal loop electromagnetic anomalies.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 5","pages":"1566-1586"},"PeriodicalIF":1.8,"publicationDate":"2025-04-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143950265","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Two-dimensional transient-extracting transform for characterization of complex tight channel reservoirs","authors":"Zhang Jia,&nbsp;Chen Hui,&nbsp;Hu Ying,&nbsp;Chen Xuping,&nbsp;Xie Yutao","doi":"10.1111/1365-2478.70023","DOIUrl":"https://doi.org/10.1111/1365-2478.70023","url":null,"abstract":"<div>\u0000 \u0000 <p>As one of the most promising types of tight oil and gas reservoirs for continental sedimentation in China, complex tight channel reservoirs have important research significance. However, this type of reservoir is characterized by thin sand bodies, vertical multistage development, complex lateral changes and tight lithology, making their exploration and development extremely difficult. We propose a new time–frequency analysis method termed the two-dimensional transient-extracting transform, which introduces a space–time window to capture the lateral information of seismic signals. This paper defines a two-dimensional transient-extracting operator, which provides highly concentrated time–frequency representations for transient signals with time–frequency curves almost parallel to the time and space axes. Then, on the basis of the relationship between the frequency and wavenumber, the wavenumber parameter is determined through two-dimensional transient-extracting transform energy maximization. Finally, a laterally continuous and energy-concentrated time–frequency spectrum is obtained. Compared with the other methods, the proposed method results in better lateral continuity, more concentrated energy and the ability to clearly differentiate thin stacked seismic reflections, enabling precise characterization of the channels. The wedge model and actual data show that the proposed method can describe the morphology of complex tight channel sand bodies more clearly and reliably than other methods.</p>\u0000 </div>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 5","pages":"1552-1565"},"PeriodicalIF":1.8,"publicationDate":"2025-04-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143950122","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Improving generalization performance of deep learning–based seismic data interpolation","authors":"Jiho Park,&nbsp;Sooyoon Kim,&nbsp;Soon Jee Seol,&nbsp;Joongmoo Byun","doi":"10.1111/1365-2478.70020","DOIUrl":"https://doi.org/10.1111/1365-2478.70020","url":null,"abstract":"<p>Seismic data interpolation techniques are vital for preprocessing, as spatial undersampling in seismic data presents processing challenges. Recently, multiple deep learning–based interpolation techniques have emerged, each catering to distinct missing data scenarios, including regular, irregular or large gaps. However, this standardized approach can induce a creeping overfitting issue in terms of various missing types, notably undermining the generalization capability of trained deep learning models. It is worthy of serious consideration for performance generalization of deep learning–based trace interpolation in terms of various missing patterns. This study introduces an innovative approach, redefining deep learning–based seismic data interpolation to focus on enhancing generalized performance be treating unseen data. We highlight how data biases in the training dataset substantially impair interpolation performance on target data with varying features. Then we offer some guidelines to counter these biases during training dataset construction. Furthermore, we propose a versatile, single deep learning model applicable to any case of missing data in real-field situations, utilizing U-Net3+ as the backbone. Experiments using field data considering various missing scenarios reveal that our method excels in interpolating unseen target data; it does this by using an unbiased dataset, bolstering general interpolation performance. This study emphasizes the importance of a systematically designed training dataset to augment generalization in deep learning–based interpolation and indicates the need for more comprehensive research to create a universally applicable deep learning–based seismic data interpolation network for practical use.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 5","pages":"1534-1551"},"PeriodicalIF":1.8,"publicationDate":"2025-04-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.70020","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143950000","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Bayesian trans-dimensional soil behaviour type inference using conditional posterior proposals","authors":"Michael Conrad Koch,&nbsp;Kazunori Fujisawa,&nbsp;Anandaroop Ray","doi":"10.1111/1365-2478.70021","DOIUrl":"https://doi.org/10.1111/1365-2478.70021","url":null,"abstract":"<p>Identification of subsurface geological profiles is indispensable to geotechnical design and construction. Subsurface stratification through Bayesian inversion of soil behaviour type index data, obtained from cone penetration tests, is achieved through the development of a novel three-block Markov chain Monte Carlo algorithm. Working in a trans-dimensional context, where the number of layers, layer depths and soil random field parameters are unknown, the algorithm is able to estimate the range of non-unique solutions or the uncertainty of these parameters. A blocking strategy has been applied that allows for the development of a formulation that primarily involves computationally inexpensive tasks such as sampling from truncated normal and Inv-Gamma distributions and evaluation of general normal densities. Part of this strategy involves the design of a novel proposal density for jumping between parameter spaces of different dimensions in the reversible jump Markov chain Monte Carlo applied in the first block. Optimal sampling in trans-dimensional problems with a single reversible jump Markov chain using random walk Metropolis–Hastings proposals is often difficult and requires ad hoc concatenation of multiple independent chains or sophisticated methods like parallel tempering or delayed rejection. The formulation presented in this study renders the conditional posterior density over the mean of the random field representing the soil parameters to be analytical, thereby allowing the corresponding proposals to be made directly from the conditional posterior. Hence, unlike most other existing algorithms, we avoid random walks altogether by sampling from the conditional posterior distribution directly. The algorithm is validated using synthetic and real soil behaviour type index data from benchmark problems. A standard normality check of the decorrelated residuals is used as a measure to test algorithm performance. Results show that the algorithm is able to identify the soil stratification parameters and random field properties correctly and also identify their uncertainties.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 5","pages":"1510-1533"},"PeriodicalIF":1.8,"publicationDate":"2025-03-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.70021","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143950337","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Monitoring waterflooding front based on the three-dimensional direct-current method with steel casing","authors":"Jiewei Shu,&nbsp;Yunhe Liu,&nbsp;Changchun Yin,&nbsp;Xiuyan Ren,&nbsp;Bo Zhang,&nbsp;Yang Su,&nbsp;Zhihao Rong,&nbsp;Laonao Wei,&nbsp;Zhiyuan Ke","doi":"10.1111/1365-2478.70002","DOIUrl":"https://doi.org/10.1111/1365-2478.70002","url":null,"abstract":"<p>When using electrical prospecting methods to monitor the front edge of waterflooding in the residual oil development, the electromagnetic responses generated by the high conductivity of the steel casing in the observation area usually obscure the electrical abnormal signal of the formation caused by waterflooding, largely affecting the imaging accuracy. Considering the densely distributed casing wells in the working area, we propose a novel method that steel casing serves as paths not only for transmitting, but also transferring the waterflooding front signals to the observations that are carried out on the surface. As the electric signals going through the steel casing have nearly zero attenuation, we can effectively observe the weak electric signals of underground waterflooding fluids at the surface receivers. To achieve a high-precision simulation of steel casings, we used an electrical conductivity equivalent method to overcome the problem of excessively large comparative geometrical dimensions. Meanwhile, we utilized unstructured tetrahedral grids to accurately model the irregular shape of waterflooding fluids migration. Based on numerical experiments, we analysed the effects of source type, location, and waterflooding depth on the conduction effect of steel casing, confirming that our method can effectively monitor the migration of waterflooding and capture the boundary of the waterflooding front edge.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 5","pages":"1471-1487"},"PeriodicalIF":1.8,"publicationDate":"2025-03-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143949723","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Prediction of anisotropic coalbed methane reservoir brittleness index based on rock physics modelling","authors":"Zhaoji Zhang,&nbsp;Fei Gong,&nbsp;Guangui Zou,&nbsp;Guowei Zhu,&nbsp;Wei Liu,&nbsp;Hao Chen","doi":"10.1111/1365-2478.70022","DOIUrl":"https://doi.org/10.1111/1365-2478.70022","url":null,"abstract":"<p>The brittleness index is a significant indicator for forecasting the characteristics of coalbed methane reservoirs. However, the brittleness index of coalbed methane reservoirs is susceptible to a number of influencing factors, including the pore structure, coalbed methane content and pore fluid. Moreover, coalbed methane reservoirs are characterized by vertical transverse isotropy anisotropy, which presents a significant challenge to the accurate prediction of brittleness. Accordingly, this paper presents a new rock physics model for coalbed methane reservoirs, based on experimental test data obtained from coal samples. This model considers the vertical transverse isotropy anisotropy characteristics of coalbed methane reservoirs and the influence of pore structure, adsorbed gas and pore fluid in coal. The accuracy of the model is corroborated by the results of the experimental tests. The model predictions indicate that organic matter and clay content exert a greater influence than the stratification indicator factor. The impact of content is smaller than that of pore structure, and the influence of coalbed methane increases following water saturation. The model may be employed to forecast the elastic anisotropy and brittleness index of vertical transverse isotropy-type coalbed methane reservoirs. The inversion results of the logging data indicate that the brittleness index of coalbed methane reservoirs exhibits some variation in different directions and the brittleness index is greater along the bedding direction. The study assists in elucidating the interrelationship between the rock physical parameters of coalbed methane reservoirs and the brittleness index of coal seams, thereby furnishing a basis for anticipating the anisotropic sweet spot in coalbed methane reservoirs.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 5","pages":"1488-1509"},"PeriodicalIF":1.8,"publicationDate":"2025-03-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143949724","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Determining structural boundaries using blocky magnetotelluric data inversion","authors":"Kaijun Xu,&nbsp;Yaoguo Li","doi":"10.1111/1365-2478.70018","DOIUrl":"https://doi.org/10.1111/1365-2478.70018","url":null,"abstract":"<div>\u0000 \u0000 <p>The magnetotelluric method has large depths of investigation and can provide important structural information in many exploration problems. The one-dimensional magnetotelluric inversion also has been applied to extract boundary information and provide the constraints for the interpretation of complementary datasets. Traditional smooth inversion based on the <span></span><math>\u0000 <semantics>\u0000 <msub>\u0000 <mi>L</mi>\u0000 <mn>2</mn>\u0000 </msub>\u0000 <annotation>${L}_2$</annotation>\u0000 </semantics></math> norm only provides a single smooth model that it is difficult to detect the location of the geological boundary. Trans-dimensional inversion provides an effective means to determine the boundaries with uncertainty quantification but incurs significant computational costs. We present an efficient method to detect distinct interfaces from one-dimensional blocky magnetotelluric inversions using an Ekblom norm. The method leverages the Ekblom norm to assess the change in the recovered resistivity model with the threshold parameter as a means to delineate the significant boundaries in the subsurface. The threshold parameter specific to the Ekblom-norm inversion is then used to probe the variability of the inversion to obtain a more robust interface detection. Once the interfaces are detected, we calculate the average resistivity value between detected interfaces to form a final conductivity model. As a demonstration, we apply this method to a synthetic example and the field data from East Tennant in Australia. The results show that the method is effective in obtaining boundaries.</p></div>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 5","pages":"1456-1470"},"PeriodicalIF":1.8,"publicationDate":"2025-03-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143950280","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Stable Q-compensated viscoelastic reverse-time migration based on the modified fractional Laplacian wave equations","authors":"Haiqiang Tang,&nbsp;Jianping Huang,&nbsp;Qiang Mao,&nbsp;Xinru Mu,&nbsp;Fei Li,&nbsp;Juan Chen","doi":"10.1111/1365-2478.70012","DOIUrl":"https://doi.org/10.1111/1365-2478.70012","url":null,"abstract":"<p>The attenuation property of earth media can lead to amplitude loss and phase dispersion effects on multicomponent elastic data. Ignoring their impacts during imaging process will result in blurred and dislocated imaging profiles. To accurately characterize the attenuation effect in viscoelastic media, we first derive a new viscoelastic wave equation with decoupled fractional Laplacians. Numerical tests show that the proposed wave equation can accurately capture the propagation characteristics of seismic waves in viscoelastic media. Furthermore, our new wave equation can be modified to yield a decomposition equation, which enables the separated propagation of vector P- and S-wavefields. Building on the derived viscoelastic forward propagator, we develop a stable <i>Q</i>-compensated viscoelastic reverse-time migration approach. Usually, the inner product imaging condition is used to obtain imaging results. However, the result of inner product is affected by the angle between vectors, making the resulting images contaminated with the angle information. In this article, we introduce the magnitude- and sign-based imaging condition for PS imaging and conduct a cross-correlation imaging condition based on the scalar P-wavefield for PP imaging. In contrast to the inner product imaging condition, our imaging scheme is capable of overcoming the contamination by the angle information. In addition, high-frequency noise is amplified exponentially during the attenuation compensation process, affecting imaging precision. To address this problem, we derive the stabilized <i>Q</i>-compensation wave equations explicitly for vector- and scalar wavefields. Numerical examples demonstrate that the proposed <i>Q</i>-compensated viscoelastic reverse-time migration method can effectively correct the viscoelastic effects, yielding high-quality PP- and PS-imaging profiles.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 5","pages":"1431-1455"},"PeriodicalIF":1.8,"publicationDate":"2025-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143950324","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}