To fully exploit the technical advantages of the large-depth and high-precision artificial source electromagnetic method in the complex structure area of southern Sichuan and compensate for the shortcomings of the con...To fully exploit the technical advantages of the large-depth and high-precision artificial source electromagnetic method in the complex structure area of southern Sichuan and compensate for the shortcomings of the conventional electromagnetic method in exploration depth,precision,and accuracy,the large-depth and high-precision wide field electromagnetic method is applied to the complex structure test area of the Luochang syncline and Yuhe nose anticline in the southern Sichuan.The advantages of the wide field electromagnetic method in detecting deep,low-resistivity thin layers are demonstrated.First,on the basis of the analysis of physical property data,a geological–geoelectric model is established in the test area,and the wide field electromagnetic method is numerically simulated to analyze and evaluate the response characteristics of deep thin shale gas layers on wide field electromagnetic curves.Second,a wide field electromagnetic test is conducted in the complex structure area of southern Sichuan.After data processing and inversion imaging,apparent resistivity logging data are used for calibration to develop an apparent resistivity interpretation model suitable for the test area.On the basis of the results,the characteristics of the electrical structure change in the shallow longitudinal formation of 6 km are implemented,and the transverse electrical distribution characteristics of the deep shale gas layer are delineated.In the prediction area near the well,the subsequent data verification shows that the apparent resistivity obtained using the inversion of the wide field electromagnetic method is consistent with the trend of apparent resistivity revealed by logging,which proves that this method can effectively identify the weak response characteristics of deep shale gas formations in complex structural areas.This experiment,it is shown shows that the wide field electromagnetic method with a large depth and high precision can effectively characterize the electrical characteristics of deep,low-resistivity thin layers in complex structural areas,and a new set of low-cost evaluation technologies for shale gas target layers based on the wide field electromagnetic method is explored.展开更多
We review the latest aeromagnetic geological data of continental China. We discuss the latest achievements in geological mapping and the newly detected features based on aeromagnetic data. Using aeromagnetic data coll...We review the latest aeromagnetic geological data of continental China. We discuss the latest achievements in geological mapping and the newly detected features based on aeromagnetic data. Using aeromagnetic data collected for more than 50 years, a series of 1:5000000 and 1:1000000 aeromagnetic maps of continental China were compiled using state-of-the-art digital technology, and data processing and transformation. Guided by plate tectonics and continental dynamics, rock physical properties, and magnetic anomalies, we compiled maps of the depth of the magnetic basement of continental China and the major geotectonic units, and presented newly detected geological structures based on the aeromagnefic data.展开更多
An electromagnetic field is generated through the accelerating movement of two equal but opposite charges of a single dipole. An electromagnetic field can also be generated by a time-varying infinitesimal point charge...An electromagnetic field is generated through the accelerating movement of two equal but opposite charges of a single dipole. An electromagnetic field can also be generated by a time-varying infinitesimal point charge. In this study, a comparison between the electromagnetic fields of an infinitesimal point charge and a dipole has been presented. First, the time-domain potential function of a point source in a 3D conductive medium is derived. Then the electric and magnetic fields in a 3D homogeneous lossless space are derived via the relation between the potential and field. The field differences between the infinitesimal point charge and the dipole in the step-off time, far-source, and near-source zones are analyzed, and the accuracy of the solutions from these sources is investigated. It is also shown that the field of the infinitesimal point charge in the near-source zone is different from that of the dipole, whereas the far-source zone fields of these two sources are identical. The comparison of real and simulated data shows that the infinitesimal point charge represents the real source better than the divole source.展开更多
Traditionally, airborne time-domain electromagnetic (ATEM) data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and ...Traditionally, airborne time-domain electromagnetic (ATEM) data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and over-determined problems in the inversion. The correlation complicates the mapping relation between the ATEM data and the earth parameters and thus increases the inversion complexity. To obviate this, we adopt principal component analysis to transform ATEM data into orthogonal principal components (PCs) to reduce the correlations and the data dimensionality and simultaneously suppress the unrelated noise. In this paper, we use an artificial neural network (ANN) to approach the PCs mapping relation with the earth model parameters, avoiding the calculation of Jacobian derivatives. The PC-based ANN algorithm is applied to synthetic data for layered models compared with data-based ANN for airborne time-domain electromagnetic inversion. The results demonstrate the PC-based ANN advantages of simpler network structure, less training steps, and better inversion results over data-based ANN, especially for contaminated data. Furthermore, the PC-based ANN algorithm effectiveness is examined by the inversion of the pseudo 2D model and comparison with data-based ANN and Zhody's methods. The results indicate that PC-based ANN inversion can achieve a better agreement with the true model and also proved that PC-based ANN is feasible to invert large ATEM datasets.展开更多
The airborne electromagnetic (AEM) method has a high sampling rate and survey flexibility. However, traditional numerical modeling approaches must use high-resolution physical grids to guarantee modeling accuracy, e...The airborne electromagnetic (AEM) method has a high sampling rate and survey flexibility. However, traditional numerical modeling approaches must use high-resolution physical grids to guarantee modeling accuracy, especially for complex geological structures such as anisotropic earth. This can lead to huge computational costs. To solve this problem, we propose a spectral-element (SE) method for 3D AEM anisotropic modeling, which combines the advantages of spectral and finite-element methods. Thus, the SE method has accuracy as high as that of the spectral method and the ability to model complex geology inherited from the finite-element method. The SE method can improve the modeling accuracy within discrete grids and reduce the dependence of modeling results on the grids. This helps achieve high-accuracy anisotropic AEM modeling. We first introduced a rotating tensor of anisotropic conductivity to Maxwell's equations and described the electrical field via SE basis functions based on GLL interpolation polynomials. We used the Galerkin weighted residual method to establish the linear equation system for the SE method, and we took a vertical magnetic dipole as the transmission source for our AEM modeling. We then applied fourth-order SE calculations with coarse physical grids to check the accuracy of our modeling results against a 1D semi-analytical solution for an anisotropic half-space model and verified the high accuracy of the SE. Moreover, we conducted AEM modeling for different anisotropic 3D abnormal bodies using two physical grid scales and three orders of SE to obtain the convergence conditions for different anisotropic abnormal bodies. Finally, we studied the identification of anisotropy for single anisotropic abnormal bodies, anisotropic surrounding rock, and single anisotropic abnormal body embedded in an anisotropic surrounding rock. This approach will play a key role in the inversion and interpretation of AEM data collected in regions with anisotropic geology.展开更多
To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional (2D) time- domain transient electromagnetic secondary field of t...To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional (2D) time- domain transient electromagnetic secondary field of the line source based on the DuFort- Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force (Emf), and determined the forward time step. By extending upward the earth-air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the fiat layered model with high precision--the maximum relative error is less than 0.01% between our method and the analytical method--and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.展开更多
Nuclear magnetic resonance measurements in hydrate-bearing sandstone samples from the Shenhu area, South China Sea were used to study the effect of gas hydrates on the sandstone permeability. The hydrate-bearing sampl...Nuclear magnetic resonance measurements in hydrate-bearing sandstone samples from the Shenhu area, South China Sea were used to study the effect of gas hydrates on the sandstone permeability. The hydrate-bearing samples contain pore-filling hydrates. The data show that the pore-filling hydrates greatly affect the formation permeability while depending on many factors that also bear on permeability; furthermore, with increasing hydrate saturation, the formation permeability decreases. We used the Masuda model and an exponent N = 7.9718 to formulate the empirical equation that describes the relation between relative permeability and hydrate saturation for the Shenhu area samples.展开更多
In an effort to reduce the shale gas exploration risks and costs, we applied the wide-field electromagnetic method (WFEM), because of its strong anti-interference capability, high resolution, ability to conduct expl...In an effort to reduce the shale gas exploration risks and costs, we applied the wide-field electromagnetic method (WFEM), because of its strong anti-interference capability, high resolution, ability to conduct exploration at large depths, and high efficiency, to the Bayan Syncline in the South Huayuan block, Hunan Province. We collected rock samples and analyzed their resistivity and induced polarization (IP) and built A series of two-dimensional models for geological conditions to investigate the applicability of WFEM to different geological structures. We also analyzed the correlation between TOC of shale and the resistivity and IP ratio to determine the threshold for identifying target formations. We used WFEM to identify the underground structures and determine the distribution, depth, and thickness of the target strata. Resistivity, IP, and total organic carbon were used to evaluate the shale gas prospects and select favorable areas (sweet spots) for exploration and development. Subsequently, drilling in these areas proved the applicability of WFEM in shale gas exploration.展开更多
We study the entanglement in anisotropie (1/2,1) mixed-spin Heisenberg XY model under the presence of an external magnetic field at thermal equilibrium. By adjusting the anisotropic parameter and the magnetic field,...We study the entanglement in anisotropie (1/2,1) mixed-spin Heisenberg XY model under the presence of an external magnetic field at thermal equilibrium. By adjusting the anisotropic parameter and the magnetic field, one is able to obtain entanglement at higher temperature. We find the evidence of the quantum phase transition in the model and observe that the quantum phase transition point at low temperature moves toward weak magnetic field with the increase of the anisotropic parameter.展开更多
Electrically anisotropic strata are abundant in nature, so their study can help our data interpretation and our understanding of the processes of geodynamics. However, current data processing generally assumes isotrop...Electrically anisotropic strata are abundant in nature, so their study can help our data interpretation and our understanding of the processes of geodynamics. However, current data processing generally assumes isotropic conditions when surveying anisotropic structures, which may cause discrepancies between reality and electromagnetic data interpretation. Moreover, the anisotropic interpretation of the time-domain airborne electromagnetic (TDAEM) method is still confined to one dimensional (1D) cases, and the corresponding three-dimensional (3D) numerical simulations are still in development. In this study, we expanded the 3D TDAEM modeling of arbitrarily anisotropic media. First, through coordinate rotation of isotropic conductivity, we obtained the conductivity tensor of an arbitrary anisotropic rock. Next, we incorporated this into Maxwell's equations, using a regular hexahedral grid of vector finite elements to subdivide the solution area. A direct solver software package provided the solution for the sparse linear equations that resulted. Analytical solutions were used to verify the accuracy and feasibility of the algorithm. The proven model was then applied to analyze the effects of arbitrary anisotropy in 3D TDAEM via the distribution of responses and amplitude changes, which revealed that different anisotropy situations strongly affected the responses of TDAEM.展开更多
Due to the strong electromagnetic interferences and human interference,traditional electromagnetic methods cannot obtain high quality resistivity data of mineral deposits in Chinese mines.The wide field electromagneti...Due to the strong electromagnetic interferences and human interference,traditional electromagnetic methods cannot obtain high quality resistivity data of mineral deposits in Chinese mines.The wide field electromagnetic method(WFEM),in which the pseudo-random signal is taken as the transmitter source,can extract high quality resistivity data in areas with sever interference by only measuring the electric field component.We use the WFEM to extract the resistivity information of the Dongguashan mine in southeast China.Compared with the audio magnetotelluric(AMT)method,and the controlled source audio-frequency magnetotelluric(CSAMT) method,the WFEM can obtain data with higher quality and simpler operations.The inversion results indicate that the WFEM can accurately identify the location of the main ore-body,which can be used for deep mine exploration in areas with strong interference.展开更多
Frequency-domain airborne electromagnetics is a proven geophysical exploration method.Presently,the interpretation is mainly based on resistivity-depth imaging and onedimensional layered inversion;nevertheless,it is d...Frequency-domain airborne electromagnetics is a proven geophysical exploration method.Presently,the interpretation is mainly based on resistivity-depth imaging and onedimensional layered inversion;nevertheless,it is difficult to obtain satisfactory results for two- or three-dimensional complex earth structures using 1D methods.3D forward modeling and inversion can be used but are hampered by computational limitations because of the large number of data.Thus,we developed a 2.5D frequency-domain airborne electromagnetic forward modeling and inversion algorithm.To eliminate the source singularities in the numerical simulations,we split the fields into primary and secondary fields.The primary fields are calculated using homogeneous or layered models with analytical solutions,and the secondary(scattered) fields are solved by the finite-element method.The linear system of equations is solved by using the large-scale sparse matrix parallel direct solver,which greatly improves the computational efficiency.The inversion algorithm was based on damping leastsquares and singular value decomposition and combined the pseudo forward modeling and reciprocity principle to compute the Jacobian matrix.Synthetic and field data were used to test the effectiveness of the proposed method.展开更多
The vector transformation and pole reduction from the total-field anomaly are signifi cant for the interpretation.We examined these industry-standard processing procedures in the Fourier domain.We propose a novel iter...The vector transformation and pole reduction from the total-field anomaly are signifi cant for the interpretation.We examined these industry-standard processing procedures in the Fourier domain.We propose a novel iteration algorithm for regional magnetic anomalies transformations to derive the vertical-component data from the total-field measurements with the variation in the core-fi eld direction over the region.Additionally,we use the same algorithm to convert the calculated vertical-component data into the corresponding data at the pole and realize the processing of diff erential reduction to the pole(DRTP).Unlike Arkani-Hamed’s DRTP method,the two types of iterative algorithms have the same forms,and DRTP is realized by implementing this algorithm twice.The synthetic model’s calculation results show that the method has high accuracy,and the fi eld data processing confi rms its practicality.展开更多
To improve the inversion accuracy of time-domain airborne electromagnetic data, we propose a parallel 3D inversion algorithm for airborne EM data based on the direct Gauss-Newton optimization. Forward modeling is perf...To improve the inversion accuracy of time-domain airborne electromagnetic data, we propose a parallel 3D inversion algorithm for airborne EM data based on the direct Gauss-Newton optimization. Forward modeling is performed in the frequency domain based on the scattered secondary electrical field. Then, the inverse Fourier transform and convolution of the transmitting waveform are used to calculate the EM responses and the sensitivity matrix in the time domain for arbitrary transmitting waves. To optimize the computational time and memory requirements, we use the EM "footprint" concept to reduce the model size and obtain the sparse sensitivity matrix. To improve the 3D inversion, we use the OpenMP library and parallel computing. We test the proposed 3D parallel inversion code using two synthetic datasets and a field dataset. The time-domain airborne EM inversion results suggest that the proposed algorithm is effective, efficient, and practical.展开更多
In order to study the distribution of shale gas reservoir in the Babaoshan Basin of Eastern Kunlun,the wide-field electromagnetic(WFEM)survey was carried out to obtain the spatial distribution characteristics of the u...In order to study the distribution of shale gas reservoir in the Babaoshan Basin of Eastern Kunlun,the wide-field electromagnetic(WFEM)survey was carried out to obtain the spatial distribution characteristics of the underground electrical volume resistivity based on the delineation of the scope of the Babaoshan Basin by regional gravity data.The basic characteristics of the basement,basin framework,and extension,vertical change,burial depth of dark mud shale in this area were identified,and the electrical distribution of the Babaoshan mud shale horizon was revealed,which has been proved to be a good geological effect by drilling.The exploration results show that the WFEM has significant effects on the exploration of shale gas occurrence strata,which meets the needs of investigation and evaluation of multi-layered and large-scale shale gas,and plays a good demonstration role in the follow-up shale gas exploration.展开更多
To make three-dimensional electromagnetic exploration achievable,the distributed wide field electromagnetic method(WFEM)based on the high-order 2^(n) sequence pseudo-random signal is proposed and realized.In this meth...To make three-dimensional electromagnetic exploration achievable,the distributed wide field electromagnetic method(WFEM)based on the high-order 2^(n) sequence pseudo-random signal is proposed and realized.In this method,only one set of high-order pseudo-random waveforms,which contains all target frequencies,is needed.Based on high-order sequence pseudo-random signal construction algorithm,the waveform can be customized according to different exploration tasks.And the receivers are independent with each other and dynamically adjust the acquisition parameters according to different requirements.A field test in the deep iron ore of Qihe−Yucheng showed that the distributed WFEM based on high-order pseudo-random signal realizes the high-efficiency acquisition of massive electromagnetic data in quite a short time.Compared with traditional controlled-source electromagnetic methods,the distributed WFEM is much more efficient.Distributed WFEM can be applied to the large scale and high-resolution exploration for deep resources and minerals.展开更多
A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been deve...A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.展开更多
The wide-field electromagnetic method is widely used in hydrocarbon exploration,mineral deposit detection,and geological disaster prediction.However,apparent resistivity and normalized field amplitude exceeding 2048 H...The wide-field electromagnetic method is widely used in hydrocarbon exploration,mineral deposit detection,and geological disaster prediction.However,apparent resistivity and normalized field amplitude exceeding 2048 Hz often exhibit upward warping in data,making geophysical inversion and interpretation challenging.The cumulative error of the crystal oscillator in signal transmission and acquisition contributes to an upturned apparent resistivity curve.To address this,a high-frequency information extraction method is proposed based on time-domain signal reconstruction,which helps to record a complete current data sequence;moreover,it helps estimate the crystal oscillator error for the transmitted signal.Considering the recorded error,a received signal was corrected using a set of reconstruction algorithms.After processing,the high-frequency component of the wide-field electromagnetic data was not upturned,while accurate high-frequency information was extracted from the signal.Therefore,the proposed method helped effectively extract high-frequency components of all wide-field electromagnetic data.展开更多
Cylindrically symmetric inhomogeneous cosmological model for bulk viscous fluid distribution with electro- magnetic field is obtained. The source of the magnetic field is due to an electric current produced along the ...Cylindrically symmetric inhomogeneous cosmological model for bulk viscous fluid distribution with electro- magnetic field is obtained. The source of the magnetic field is due to an electric current produced along the z-axis. F12 is the non-vanishing component of electromagnetic field tensor. To get the deterministic solution, it has been assumed that the expansion 0 in the model is proportional to the shear σ. The values of cosmological constant for these models are found to be small and positive at late time, which are consistent with the results from recent supernovae Ia observations. Physical and geometric aspects of the models are also discussed in presence and absence of magnetic field.展开更多
文摘To fully exploit the technical advantages of the large-depth and high-precision artificial source electromagnetic method in the complex structure area of southern Sichuan and compensate for the shortcomings of the conventional electromagnetic method in exploration depth,precision,and accuracy,the large-depth and high-precision wide field electromagnetic method is applied to the complex structure test area of the Luochang syncline and Yuhe nose anticline in the southern Sichuan.The advantages of the wide field electromagnetic method in detecting deep,low-resistivity thin layers are demonstrated.First,on the basis of the analysis of physical property data,a geological–geoelectric model is established in the test area,and the wide field electromagnetic method is numerically simulated to analyze and evaluate the response characteristics of deep thin shale gas layers on wide field electromagnetic curves.Second,a wide field electromagnetic test is conducted in the complex structure area of southern Sichuan.After data processing and inversion imaging,apparent resistivity logging data are used for calibration to develop an apparent resistivity interpretation model suitable for the test area.On the basis of the results,the characteristics of the electrical structure change in the shallow longitudinal formation of 6 km are implemented,and the transverse electrical distribution characteristics of the deep shale gas layer are delineated.In the prediction area near the well,the subsequent data verification shows that the apparent resistivity obtained using the inversion of the wide field electromagnetic method is consistent with the trend of apparent resistivity revealed by logging,which proves that this method can effectively identify the weak response characteristics of deep shale gas formations in complex structural areas.This experiment,it is shown shows that the wide field electromagnetic method with a large depth and high precision can effectively characterize the electrical characteristics of deep,low-resistivity thin layers in complex structural areas,and a new set of low-cost evaluation technologies for shale gas target layers based on the wide field electromagnetic method is explored.
基金supported by the China Land Aeromagnetic Characteristics and Tectonic Structures Research(No.1212011087009)part of the national geological and mineral resources investigation projects,and the Comprehensive Exploration of Aero Geophysical&Remote Sensing Survey and Interpretation System Research(No.2013AA063905)part of the planning for national high technology research and development
文摘We review the latest aeromagnetic geological data of continental China. We discuss the latest achievements in geological mapping and the newly detected features based on aeromagnetic data. Using aeromagnetic data collected for more than 50 years, a series of 1:5000000 and 1:1000000 aeromagnetic maps of continental China were compiled using state-of-the-art digital technology, and data processing and transformation. Guided by plate tectonics and continental dynamics, rock physical properties, and magnetic anomalies, we compiled maps of the depth of the magnetic basement of continental China and the major geotectonic units, and presented newly detected geological structures based on the aeromagnefic data.
基金supported by Chinese National Programs for Fundamental Research and Development(No.2012CB416605)the National Natural Science Foundation of China(No.41174090)Development Project of National Key Scientific Equipment(No.ZDYZ2012-1-05-04)
文摘An electromagnetic field is generated through the accelerating movement of two equal but opposite charges of a single dipole. An electromagnetic field can also be generated by a time-varying infinitesimal point charge. In this study, a comparison between the electromagnetic fields of an infinitesimal point charge and a dipole has been presented. First, the time-domain potential function of a point source in a 3D conductive medium is derived. Then the electric and magnetic fields in a 3D homogeneous lossless space are derived via the relation between the potential and field. The field differences between the infinitesimal point charge and the dipole in the step-off time, far-source, and near-source zones are analyzed, and the accuracy of the solutions from these sources is investigated. It is also shown that the field of the infinitesimal point charge in the near-source zone is different from that of the dipole, whereas the far-source zone fields of these two sources are identical. The comparison of real and simulated data shows that the infinitesimal point charge represents the real source better than the divole source.
基金supported by the National Natural Science Foundation of China (Grant No. 40974039)High-Tech Research and Development Program of China (Grant No.2006AA06205)Leading Strategic Project of Science and Technology, Chinese Academy of Sciences (XDA08020500)
文摘Traditionally, airborne time-domain electromagnetic (ATEM) data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and over-determined problems in the inversion. The correlation complicates the mapping relation between the ATEM data and the earth parameters and thus increases the inversion complexity. To obviate this, we adopt principal component analysis to transform ATEM data into orthogonal principal components (PCs) to reduce the correlations and the data dimensionality and simultaneously suppress the unrelated noise. In this paper, we use an artificial neural network (ANN) to approach the PCs mapping relation with the earth model parameters, avoiding the calculation of Jacobian derivatives. The PC-based ANN algorithm is applied to synthetic data for layered models compared with data-based ANN for airborne time-domain electromagnetic inversion. The results demonstrate the PC-based ANN advantages of simpler network structure, less training steps, and better inversion results over data-based ANN, especially for contaminated data. Furthermore, the PC-based ANN algorithm effectiveness is examined by the inversion of the pseudo 2D model and comparison with data-based ANN and Zhody's methods. The results indicate that PC-based ANN inversion can achieve a better agreement with the true model and also proved that PC-based ANN is feasible to invert large ATEM datasets.
基金financially supported by the Key Program of National Natural Science Foundation of China(No.41530320)China Natural Science Foundation for Young Scientists(No.41404093)+1 种基金Key National Research Project of China(Nos2016YFC0303100 and 2017YFC0601900)China Natural Science Foundation(No.41774125)
文摘The airborne electromagnetic (AEM) method has a high sampling rate and survey flexibility. However, traditional numerical modeling approaches must use high-resolution physical grids to guarantee modeling accuracy, especially for complex geological structures such as anisotropic earth. This can lead to huge computational costs. To solve this problem, we propose a spectral-element (SE) method for 3D AEM anisotropic modeling, which combines the advantages of spectral and finite-element methods. Thus, the SE method has accuracy as high as that of the spectral method and the ability to model complex geology inherited from the finite-element method. The SE method can improve the modeling accuracy within discrete grids and reduce the dependence of modeling results on the grids. This helps achieve high-accuracy anisotropic AEM modeling. We first introduced a rotating tensor of anisotropic conductivity to Maxwell's equations and described the electrical field via SE basis functions based on GLL interpolation polynomials. We used the Galerkin weighted residual method to establish the linear equation system for the SE method, and we took a vertical magnetic dipole as the transmission source for our AEM modeling. We then applied fourth-order SE calculations with coarse physical grids to check the accuracy of our modeling results against a 1D semi-analytical solution for an anisotropic half-space model and verified the high accuracy of the SE. Moreover, we conducted AEM modeling for different anisotropic 3D abnormal bodies using two physical grid scales and three orders of SE to obtain the convergence conditions for different anisotropic abnormal bodies. Finally, we studied the identification of anisotropy for single anisotropic abnormal bodies, anisotropic surrounding rock, and single anisotropic abnormal body embedded in an anisotropic surrounding rock. This approach will play a key role in the inversion and interpretation of AEM data collected in regions with anisotropic geology.
基金supported by the National High Technology Research and Development Program (863 Program)(2009AA06Z108)
文摘To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional (2D) time- domain transient electromagnetic secondary field of the line source based on the DuFort- Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force (Emf), and determined the forward time step. By extending upward the earth-air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the fiat layered model with high precision--the maximum relative error is less than 0.01% between our method and the analytical method--and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.
基金supported by the Key Program for International S&T Cooperation Projects of China(No.2010DFA21630)the National Natural Science Foundation of China(No.41306050)
文摘Nuclear magnetic resonance measurements in hydrate-bearing sandstone samples from the Shenhu area, South China Sea were used to study the effect of gas hydrates on the sandstone permeability. The hydrate-bearing samples contain pore-filling hydrates. The data show that the pore-filling hydrates greatly affect the formation permeability while depending on many factors that also bear on permeability; furthermore, with increasing hydrate saturation, the formation permeability decreases. We used the Masuda model and an exponent N = 7.9718 to formulate the empirical equation that describes the relation between relative permeability and hydrate saturation for the Shenhu area samples.
基金financially supported by the Thirteenth Five-Year-Plan Major Project "Marine Shale Gas Exploration and Evaluation over Laifengxianfeng and Hefeng Block"(No.2016ZX05034004-004)China Huadian Engineering Co.,LTD(No.CHEC-KJ-2014-Z10)
文摘In an effort to reduce the shale gas exploration risks and costs, we applied the wide-field electromagnetic method (WFEM), because of its strong anti-interference capability, high resolution, ability to conduct exploration at large depths, and high efficiency, to the Bayan Syncline in the South Huayuan block, Hunan Province. We collected rock samples and analyzed their resistivity and induced polarization (IP) and built A series of two-dimensional models for geological conditions to investigate the applicability of WFEM to different geological structures. We also analyzed the correlation between TOC of shale and the resistivity and IP ratio to determine the threshold for identifying target formations. We used WFEM to identify the underground structures and determine the distribution, depth, and thickness of the target strata. Resistivity, IP, and total organic carbon were used to evaluate the shale gas prospects and select favorable areas (sweet spots) for exploration and development. Subsequently, drilling in these areas proved the applicability of WFEM in shale gas exploration.
基金supported by National Natural Science Foundation of China under Grant No.10674181
文摘We study the entanglement in anisotropie (1/2,1) mixed-spin Heisenberg XY model under the presence of an external magnetic field at thermal equilibrium. By adjusting the anisotropic parameter and the magnetic field, one is able to obtain entanglement at higher temperature. We find the evidence of the quantum phase transition in the model and observe that the quantum phase transition point at low temperature moves toward weak magnetic field with the increase of the anisotropic parameter.
基金financially supported by National Nonprofit institute Research Grant of IGGE(Nos.AS2017J06,AS2017Y04,and AS2016J10)Survey on coastal area for airborne magnetic method of UNV in Jiangsu(No.DD20160151-03)+3 种基金Key National Research Project of China(No.2017YFC0601900)Key Program of National Natural Science Foundation of China(No.41530320)Natural Science Foundation(No.41274121)China Natural Science Foundation for Young Scientists(No.41404093)
文摘Electrically anisotropic strata are abundant in nature, so their study can help our data interpretation and our understanding of the processes of geodynamics. However, current data processing generally assumes isotropic conditions when surveying anisotropic structures, which may cause discrepancies between reality and electromagnetic data interpretation. Moreover, the anisotropic interpretation of the time-domain airborne electromagnetic (TDAEM) method is still confined to one dimensional (1D) cases, and the corresponding three-dimensional (3D) numerical simulations are still in development. In this study, we expanded the 3D TDAEM modeling of arbitrarily anisotropic media. First, through coordinate rotation of isotropic conductivity, we obtained the conductivity tensor of an arbitrary anisotropic rock. Next, we incorporated this into Maxwell's equations, using a regular hexahedral grid of vector finite elements to subdivide the solution area. A direct solver software package provided the solution for the sparse linear equations that resulted. Analytical solutions were used to verify the accuracy and feasibility of the algorithm. The proven model was then applied to analyze the effects of arbitrary anisotropy in 3D TDAEM via the distribution of responses and amplitude changes, which revealed that different anisotropy situations strongly affected the responses of TDAEM.
基金Project(2018YFC0807802)supported by the National Key R&D Program of ChinaProject(41874081)supported by the National Natural Science Foundation of China
文摘Due to the strong electromagnetic interferences and human interference,traditional electromagnetic methods cannot obtain high quality resistivity data of mineral deposits in Chinese mines.The wide field electromagnetic method(WFEM),in which the pseudo-random signal is taken as the transmitter source,can extract high quality resistivity data in areas with sever interference by only measuring the electric field component.We use the WFEM to extract the resistivity information of the Dongguashan mine in southeast China.Compared with the audio magnetotelluric(AMT)method,and the controlled source audio-frequency magnetotelluric(CSAMT) method,the WFEM can obtain data with higher quality and simpler operations.The inversion results indicate that the WFEM can accurately identify the location of the main ore-body,which can be used for deep mine exploration in areas with strong interference.
基金supported by the Doctoral Fund Project of the Ministry of Education(No.20130061110060 class tutors)the National Natural Science Foundation of China(No.41504083)National Basic Research Program of China(973Program)(No.2013CB429805)
文摘Frequency-domain airborne electromagnetics is a proven geophysical exploration method.Presently,the interpretation is mainly based on resistivity-depth imaging and onedimensional layered inversion;nevertheless,it is difficult to obtain satisfactory results for two- or three-dimensional complex earth structures using 1D methods.3D forward modeling and inversion can be used but are hampered by computational limitations because of the large number of data.Thus,we developed a 2.5D frequency-domain airborne electromagnetic forward modeling and inversion algorithm.To eliminate the source singularities in the numerical simulations,we split the fields into primary and secondary fields.The primary fields are calculated using homogeneous or layered models with analytical solutions,and the secondary(scattered) fields are solved by the finite-element method.The linear system of equations is solved by using the large-scale sparse matrix parallel direct solver,which greatly improves the computational efficiency.The inversion algorithm was based on damping leastsquares and singular value decomposition and combined the pseudo forward modeling and reciprocity principle to compute the Jacobian matrix.Synthetic and field data were used to test the effectiveness of the proposed method.
基金supported by the National Key R&D Program of China (No. 2017YFC0602000)the China Geological Survey Project (Nos. DD20191001 and DD20189410)。
文摘The vector transformation and pole reduction from the total-field anomaly are signifi cant for the interpretation.We examined these industry-standard processing procedures in the Fourier domain.We propose a novel iteration algorithm for regional magnetic anomalies transformations to derive the vertical-component data from the total-field measurements with the variation in the core-fi eld direction over the region.Additionally,we use the same algorithm to convert the calculated vertical-component data into the corresponding data at the pole and realize the processing of diff erential reduction to the pole(DRTP).Unlike Arkani-Hamed’s DRTP method,the two types of iterative algorithms have the same forms,and DRTP is realized by implementing this algorithm twice.The synthetic model’s calculation results show that the method has high accuracy,and the fi eld data processing confi rms its practicality.
基金supported by the Key Natural Science Foundation(No.41530320)Natural Science Foundation(No.41274121)+1 种基金Natural Science Foundation for young scientist(No.41404093)the Projects on the Development of the Key Equipment of Chinese Academy of Science(No.ZDYZ2012-1-03)
文摘To improve the inversion accuracy of time-domain airborne electromagnetic data, we propose a parallel 3D inversion algorithm for airborne EM data based on the direct Gauss-Newton optimization. Forward modeling is performed in the frequency domain based on the scattered secondary electrical field. Then, the inverse Fourier transform and convolution of the transmitting waveform are used to calculate the EM responses and the sensitivity matrix in the time domain for arbitrary transmitting waves. To optimize the computational time and memory requirements, we use the EM "footprint" concept to reduce the model size and obtain the sparse sensitivity matrix. To improve the 3D inversion, we use the OpenMP library and parallel computing. We test the proposed 3D parallel inversion code using two synthetic datasets and a field dataset. The time-domain airborne EM inversion results suggest that the proposed algorithm is effective, efficient, and practical.
基金Project(2019-SF-141)supported by Science and Technology Program of Qinghai Province,ChinaProjects(2017042105kc055,2017042014ky014)supported by Geological Exploration Foundation of Qinghai Province,China。
文摘In order to study the distribution of shale gas reservoir in the Babaoshan Basin of Eastern Kunlun,the wide-field electromagnetic(WFEM)survey was carried out to obtain the spatial distribution characteristics of the underground electrical volume resistivity based on the delineation of the scope of the Babaoshan Basin by regional gravity data.The basic characteristics of the basement,basin framework,and extension,vertical change,burial depth of dark mud shale in this area were identified,and the electrical distribution of the Babaoshan mud shale horizon was revealed,which has been proved to be a good geological effect by drilling.The exploration results show that the WFEM has significant effects on the exploration of shale gas occurrence strata,which meets the needs of investigation and evaluation of multi-layered and large-scale shale gas,and plays a good demonstration role in the follow-up shale gas exploration.
基金funded by the National Natural Science Foundation of China(No.42004056)the Natural Science Foundation of Shangdong Province,China(No.ZR2020QD052)China Postdoctoral Science Foundation(No.2019M652386)。
文摘To make three-dimensional electromagnetic exploration achievable,the distributed wide field electromagnetic method(WFEM)based on the high-order 2^(n) sequence pseudo-random signal is proposed and realized.In this method,only one set of high-order pseudo-random waveforms,which contains all target frequencies,is needed.Based on high-order sequence pseudo-random signal construction algorithm,the waveform can be customized according to different exploration tasks.And the receivers are independent with each other and dynamically adjust the acquisition parameters according to different requirements.A field test in the deep iron ore of Qihe−Yucheng showed that the distributed WFEM based on high-order pseudo-random signal realizes the high-efficiency acquisition of massive electromagnetic data in quite a short time.Compared with traditional controlled-source electromagnetic methods,the distributed WFEM is much more efficient.Distributed WFEM can be applied to the large scale and high-resolution exploration for deep resources and minerals.
基金Projects(41674080,41674079)supported by the National Natural Science Foundation of China
文摘A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.
基金Project(42004056)supported by the National Natural Science Foundation of ChinaProject(ZR2020QD052)supported by the Natural Science Foundation of Shandong Province,ChinaProject(2019YFC0604902)supported by the National Key Research and Development Program of China。
文摘The wide-field electromagnetic method is widely used in hydrocarbon exploration,mineral deposit detection,and geological disaster prediction.However,apparent resistivity and normalized field amplitude exceeding 2048 Hz often exhibit upward warping in data,making geophysical inversion and interpretation challenging.The cumulative error of the crystal oscillator in signal transmission and acquisition contributes to an upturned apparent resistivity curve.To address this,a high-frequency information extraction method is proposed based on time-domain signal reconstruction,which helps to record a complete current data sequence;moreover,it helps estimate the crystal oscillator error for the transmitted signal.Considering the recorded error,a received signal was corrected using a set of reconstruction algorithms.After processing,the high-frequency component of the wide-field electromagnetic data was not upturned,while accurate high-frequency information was extracted from the signal.Therefore,the proposed method helped effectively extract high-frequency components of all wide-field electromagnetic data.
文摘Cylindrically symmetric inhomogeneous cosmological model for bulk viscous fluid distribution with electro- magnetic field is obtained. The source of the magnetic field is due to an electric current produced along the z-axis. F12 is the non-vanishing component of electromagnetic field tensor. To get the deterministic solution, it has been assumed that the expansion 0 in the model is proportional to the shear σ. The values of cosmological constant for these models are found to be small and positive at late time, which are consistent with the results from recent supernovae Ia observations. Physical and geometric aspects of the models are also discussed in presence and absence of magnetic field.