In mineral exploration, the apparent resistivity and apparent frequency （or apparent polarizability） parameters of induced polarization method are commonly utilized to describe the induced polarization anomaly. When...In mineral exploration, the apparent resistivity and apparent frequency （or apparent polarizability） parameters of induced polarization method are commonly utilized to describe the induced polarization anomaly. When the target geology structure is significantly complicated, these parameters would fail to reflect the nature of the anomaly source, and wrong conclusions may be obtained. A wavelet approach and a metal factor method were used to comprehensively interpret the induced polarization anomaly of complex geologic bodies in the Adi Bladia mine. Db5 wavelet basis was used to conduct two-scale decomposition and reconstruction, which effectively suppress the noise interference of greenschist facies regional metamorphism and magma intrusion, making energy concentrated and boundary problem unobservable. On the basis of that, the ore-induced anomaly was effectively extracted by the metal factor method.展开更多
Traditional two-dimensional（2D） complex resistivity forward modeling is based on Poisson＇s equation but spectral induced polarization（SIP） data are the coproducts of the induced polarization（IP） and the electro...Traditional two-dimensional（2D） complex resistivity forward modeling is based on Poisson＇s equation but spectral induced polarization（SIP） data are the coproducts of the induced polarization（IP） and the electromagnetic induction（EMI） effects.This is especially true under high frequencies,where the EMI effect can exceed the IP effect.2D inversion that only considers the IP effect reduces the reliability of the inversion data.In this paper,we derive differential equations using Maxwell＇s equations.With the introduction of the Cole-Cole model,we use the finite-element method to conduct2 D SIP forward modeling that considers the EMI and IP effects simultaneously.The data-space Occam method,in which different constraints to the model smoothness and parametric boundaries are introduced,is then used to simultaneously obtain the four parameters of the Cole-Cole model using multi-array electric field data.This approach not only improves the stability of the inversion but also significantly reduces the solution ambiguity.To improve the computational efficiency,message passing interface programming was used to accelerate the 2D SIP forward modeling and inversion.Synthetic datasets were tested using both serial and parallel algorithms,and the tests suggest that the proposed parallel algorithm is robust and efficient.展开更多
The electromagnetically induced reflection(EIR)effect of graphene metamaterials has been investigated by finite difference time domain(FDTD)method.In this study,a metamaterial sandwich structure composed of silica(SiO...The electromagnetically induced reflection(EIR)effect of graphene metamaterials has been investigated by finite difference time domain(FDTD)method.In this study,a metamaterial sandwich structure composed of silica(SiO2),gold and graphene on terahertz band is designed.By changing the width of the two ribbons of graphene length and the incident angle of electromagnetic wave,the EIR effect of the structure is discussed,and it can be found that SiO2 is a kind of excellent dielectric material.The simulation results show that graphene metamaterial is not sensitive to polarized incident electromagnetic wave.Therefore,such EIR phenomena as insensitive polarization and large incident angle can be applied to optical communication filters and terahertz devices.展开更多
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.展开更多
The transient electromagnetic(TEM)method becomes more urgent than ever for marine ex-ploration due to abundant resource reserves and the increasing undersea engineering construction activities,especially in the offsho...The transient electromagnetic(TEM)method becomes more urgent than ever for marine ex-ploration due to abundant resource reserves and the increasing undersea engineering construction activities,especially in the offshore exploration of mineral deposits such as Sanshandao gold mine.However,the re-search and application of TEM method in marine environment are still challenged by many problems.Such contradiction motivates our study on the coincident-loop TEM in seafloor exploration.The TEM response of coincident loops is firstly derived in the integral form,based on the potential functions in Helmholtz equ-ations for a magnetic source locating in the whole-space layered model.The frequency-domain vertical magnetic field is described as the Hankel integral with double first-order Bessel functions of first kind.Se-condly,the time-domain induced voltage is obtained by transforming the frequency-domain response through the cosine transform and then taking the derivative of time.To simultaneously solve the Hankel transform and the cosine transform,a novel algorithm is introduced by adapting the fixed-point quadra-ture and extrapolation via the Shanks transformation.Finally,a typical conductivity model for marine po-lymetallic deposit is designed to investigate the characteristic of TEM response under various conditions.Numerical results demonstrate that existence of conductive seawater causes the TEM response to increase significantly and decay slower.The air-sea reflected electromagnetic waves lead to a significantly large fake negative response(NR)in shallower seawater with depth less than 300 m.Increase in the height of loops will weaken and delay the anomaly response and shorten the observation time-window.The height of configu-ration should be no more than 100 m for shallower targets and 50 m for deeper targets,respectively.The observation time-window should cover 10-1000 ms.Increase in the radius of loops only enhances the TEM response proportionally but hardly improves the relative anomaly.The vertical resolution on the low-resistivity target approximates 20 m for the configuration considered in the study.Decreases in D.C.resistivity and chargeability cause the positive response(PR)to increase significantly and decay more ra-pidly.Meanwhile,the NR is advanced and enlarged significantly and decays slower compared with the PR.The influence of time constant is not monotony and there exists an optimal value for producing the maxi-mum NR.As the frequency parameter increases,the PR is caused to decay more rapidly without magni-tude change and the NR is advanced and decays more rapidly with significant increase in magnitude.The influence of frequency parameter is more pronounced than that of time constant.展开更多
基金Project(41174103)supported by the National Natural Science Foundation of ChinaProject(2010-211)supported by the Foreign Mineral Resources Venture Exploration Special Fund of China
文摘In mineral exploration, the apparent resistivity and apparent frequency （or apparent polarizability） parameters of induced polarization method are commonly utilized to describe the induced polarization anomaly. When the target geology structure is significantly complicated, these parameters would fail to reflect the nature of the anomaly source, and wrong conclusions may be obtained. A wavelet approach and a metal factor method were used to comprehensively interpret the induced polarization anomaly of complex geologic bodies in the Adi Bladia mine. Db5 wavelet basis was used to conduct two-scale decomposition and reconstruction, which effectively suppress the noise interference of greenschist facies regional metamorphism and magma intrusion, making energy concentrated and boundary problem unobservable. On the basis of that, the ore-induced anomaly was effectively extracted by the metal factor method.
基金jointly sponsored by the National Natural Science Foundation of China(Grant No.41374078)the Geological Survey Projects of the Ministry of Land and Resources of China(Grant Nos.12120113086100 and 12120113101300)Beijing Higher Education Young Elite Teacher Project
文摘Traditional two-dimensional（2D） complex resistivity forward modeling is based on Poisson＇s equation but spectral induced polarization（SIP） data are the coproducts of the induced polarization（IP） and the electromagnetic induction（EMI） effects.This is especially true under high frequencies,where the EMI effect can exceed the IP effect.2D inversion that only considers the IP effect reduces the reliability of the inversion data.In this paper,we derive differential equations using Maxwell＇s equations.With the introduction of the Cole-Cole model,we use the finite-element method to conduct2 D SIP forward modeling that considers the EMI and IP effects simultaneously.The data-space Occam method,in which different constraints to the model smoothness and parametric boundaries are introduced,is then used to simultaneously obtain the four parameters of the Cole-Cole model using multi-array electric field data.This approach not only improves the stability of the inversion but also significantly reduces the solution ambiguity.To improve the computational efficiency,message passing interface programming was used to accelerate the 2D SIP forward modeling and inversion.Synthetic datasets were tested using both serial and parallel algorithms,and the tests suggest that the proposed parallel algorithm is robust and efficient.
基金Research Project of Anhui Province Education Department(No.KJ2020A0684)Innovation and Entrepreneurship Training Program for College Students(Nos.S201910375072,201910375050,201910375052,202010375030)。
文摘The electromagnetically induced reflection(EIR)effect of graphene metamaterials has been investigated by finite difference time domain(FDTD)method.In this study,a metamaterial sandwich structure composed of silica(SiO2),gold and graphene on terahertz band is designed.By changing the width of the two ribbons of graphene length and the incident angle of electromagnetic wave,the EIR effect of the structure is discussed,and it can be found that SiO2 is a kind of excellent dielectric material.The simulation results show that graphene metamaterial is not sensitive to polarized incident electromagnetic wave.Therefore,such EIR phenomena as insensitive polarization and large incident angle can be applied to optical communication filters and terahertz devices.
基金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.
基金This research was supported by the National Key R&D Program of China(No.2018YFC0603806).
文摘The transient electromagnetic(TEM)method becomes more urgent than ever for marine ex-ploration due to abundant resource reserves and the increasing undersea engineering construction activities,especially in the offshore exploration of mineral deposits such as Sanshandao gold mine.However,the re-search and application of TEM method in marine environment are still challenged by many problems.Such contradiction motivates our study on the coincident-loop TEM in seafloor exploration.The TEM response of coincident loops is firstly derived in the integral form,based on the potential functions in Helmholtz equ-ations for a magnetic source locating in the whole-space layered model.The frequency-domain vertical magnetic field is described as the Hankel integral with double first-order Bessel functions of first kind.Se-condly,the time-domain induced voltage is obtained by transforming the frequency-domain response through the cosine transform and then taking the derivative of time.To simultaneously solve the Hankel transform and the cosine transform,a novel algorithm is introduced by adapting the fixed-point quadra-ture and extrapolation via the Shanks transformation.Finally,a typical conductivity model for marine po-lymetallic deposit is designed to investigate the characteristic of TEM response under various conditions.Numerical results demonstrate that existence of conductive seawater causes the TEM response to increase significantly and decay slower.The air-sea reflected electromagnetic waves lead to a significantly large fake negative response(NR)in shallower seawater with depth less than 300 m.Increase in the height of loops will weaken and delay the anomaly response and shorten the observation time-window.The height of configu-ration should be no more than 100 m for shallower targets and 50 m for deeper targets,respectively.The observation time-window should cover 10-1000 ms.Increase in the radius of loops only enhances the TEM response proportionally but hardly improves the relative anomaly.The vertical resolution on the low-resistivity target approximates 20 m for the configuration considered in the study.Decreases in D.C.resistivity and chargeability cause the positive response(PR)to increase significantly and decay more ra-pidly.Meanwhile,the NR is advanced and enlarged significantly and decays slower compared with the PR.The influence of time constant is not monotony and there exists an optimal value for producing the maxi-mum NR.As the frequency parameter increases,the PR is caused to decay more rapidly without magni-tude change and the NR is advanced and decays more rapidly with significant increase in magnitude.The influence of frequency parameter is more pronounced than that of time constant.
