The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the tim...The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.展开更多
A geophysical survey was conducted in the Kelle-Bidjocka village, Messondo subdivision, in the Centre Region, Cameroon. The data acquisition was made by combining Schlumberger profiling and electrical soundings along ...A geophysical survey was conducted in the Kelle-Bidjocka village, Messondo subdivision, in the Centre Region, Cameroon. The data acquisition was made by combining Schlumberger profiling and electrical soundings along six (06) profiles of 1500 m in length for a total of 64 geoelectrical stations’ survey conducted through a variable mesh 100 m × 200 m, or 100 m × 300 m. The equipment used is the DC resistivimeter Syscal Junior 48 (Iris Instrument). Processing and modelling of field data are made by using the Res2Dinv, Qwseln and Surfer software. The investigation methods used are electrical resistivity (DC) and induced polarization (IP) methods. The analyses and interpretations have helped to highlight areas of weakness or conductive discontinuities (fractures, faults, shear zones, etc.) in Precambrian gneiss formations, sometimes undergoing weathering processes. They identify the weathering or mineralogical accumulation horizons, the most promising is a mineralization channel identified in the NE-SW direction. The highlighted mineralization is characterized by strong gradients of chargeability or polarization. Samples and other geological evidences observed in the area are used to associate the most polarizable structures with ferriferous formations. Weakly polarizable and particularly conductive backgrounds identified by the inverse pseudo-sections are thought to be sulphate minerals or groundwater targets for future hydrogeological studies.展开更多
Earth medium is not completely elastic, with its viscosity resulting in attenuation and dispersion of seismic waves. Most viscoelastic numerical simulations are based on the finite-difference and finite-element method...Earth medium is not completely elastic, with its viscosity resulting in attenuation and dispersion of seismic waves. Most viscoelastic numerical simulations are based on the finite-difference and finite-element methods. Targeted at viscoelastic numerical modeling for multilayered media, the constant-Q acoustic wave equation is transformed into the corresponding wave integral representation with its Green's function accounting for viscoelastic coefficients. An efficient alternative for full-waveform solution to the integral equation is proposed in this article by extending conventional frequency-domain boundary element methods to viscoelastic media. The viscoelastic boundary element method enjoys a distinct characteristic of the explicit use of boundary continuity conditions of displacement and traction, leading to a semi-analytical solution with sufficient accuracy for simulating the viscoelastic effect across irregular interfaces. Numerical experiments to study the viscoelastic absorption of different Q values demonstrate the accuracy and applicability of the method.展开更多
In this paper the authors propose the natural field induced polarization (IP) method and present the way to pick up IP effect. The relations between the object and anomaly are studied by taking field experiments as ex...In this paper the authors propose the natural field induced polarization (IP) method and present the way to pick up IP effect. The relations between the object and anomaly are studied by taking field experiments as examples. The effectiveness and usability of the method are testified.展开更多
To achieve the current study, geoelectrical surveys along six (06) profiles of 4 km long in a 100 m × 200 m grid defined according to the triangulation principle in the Ngoura area (Tindikala-Boutou villages) hav...To achieve the current study, geoelectrical surveys along six (06) profiles of 4 km long in a 100 m × 200 m grid defined according to the triangulation principle in the Ngoura area (Tindikala-Boutou villages) have been made through electrical sounding and profiling following Schlumberger array. The instrument is the resistivimeter Syscal Junior 48 (IRIS Instrument) which uses the electrical current. The data have been processed and modelled with Res2Dinv and Winsev softwares then interpolated with Surfer software. Electrical methods used are the Direct current (DC) and the Induced Polarization (IP). Interpretation and analyses of results from each investigation method highlight weak zones or conductive discontinuities and mineralized zones. Conductive zones have been identified as shear zones within granitic structures of the Precambrian basement, according to the geologic and tectonic background of the area. The structural trend of these shear zones is E-W. The mineralization within it is N-S and characterized by high values of chargeability, essentially in the eastern part of the area under study. This mineralization proves the presence of metalliferous or sulphide heaps disseminated in weathered quartz veins which cross shear zones. Also, the poor mineralization and conductive structures in shear zones characterize the groundwater zones. The intense activities of gold washers encountered in the area enable to link that mineralization to gold within quartz veins. The near surface gold mineralization is eluvial or alluvial, and in depth this mineralization is primary.展开更多
A geophysical investigation of subsurface structures using the Syscal Junior 48 resistivity-meter was conducted in Ngoura subdivision (East Cameroon) following a combined geoelectrical direct current (DC) approach inv...A geophysical investigation of subsurface structures using the Syscal Junior 48 resistivity-meter was conducted in Ngoura subdivision (East Cameroon) following a combined geoelectrical direct current (DC) approach involving Resistivity and IP methods. This investigation was allowed to collect data on forty-five (45) profiling lines at three acquisition levels (AB = 100 m, MN = 10 m;AB = 200 m, MN = 20 m and AB = 500 m, MN = 50 m) and two electric panels L1 and L4, using the Schlumberger array. Processing, modeling and interpretation of data collected using the Winsev, Res2Dinv and Surfer software helped in highlighting a conductive and strongly mineralized discontinuity in granite formations, which lined up with the NE-SW Kadei tectonic line. It extends beyond 100 m depth over an average width of 600 m. The mineralization associated with this discontinuity is identified by a high concentration of disseminated metalliferous minerals in quartz or pegmatite veins. The mining reconnaissance works in the study area and those of several authors have characterized this anomaly to a lode gold quartz or large pegmatite. The results of this study correlate with geological and tectonic data for the region marked by NE-SW Kadei tectonic line. Therefore, they confirm the reliability of a geoelectrical DC investigation method combining Resistivity and IP to the identification of ore bodies.展开更多
在网络世界中,网际互连协议(Internet Protocol,IP)承载网扮演着信息流动的骨干角色,其运行表现和服务水准至关重要。文章详细探讨了IP承载网中服务质量(Quality of Service,QoS)策略的设计与优化方法,重点涵盖分类服务与标记、流量控...在网络世界中,网际互连协议(Internet Protocol,IP)承载网扮演着信息流动的骨干角色,其运行表现和服务水准至关重要。文章详细探讨了IP承载网中服务质量(Quality of Service,QoS)策略的设计与优化方法,重点涵盖分类服务与标记、流量控制与优先级队列管理以及链路和节点策略的具体实现。通过流量分析与预测、智能资源分配与调度、动态策略调整以及故障检测与容错机制,提出一系列提升网络性能和用户体验的有效措施。展开更多
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.展开更多
During transient electromagnetic method (TEM) exploration of a copper mine, we detected the late-channel TEM signal reversal phenomenon (a voltage change from positive to negative) caused by the influence of the i...During transient electromagnetic method (TEM) exploration of a copper mine, we detected the late-channel TEM signal reversal phenomenon (a voltage change from positive to negative) caused by the influence of the induced polarization (IP) effect, which affects the depth and precision of the TEM detection. The conventional inversion method is inefficient because it is difficult to process the data. In this paper, the Cole-Cole model is adopted to analyze the effect of Dc resistivity, chargeability, time constant, and frequency exponent on the TEM response in an homogeneous half space model. Singular Value Decomposition (SVD) is used to invert the measured TEM data, and the Dc resistivity, chargeability, time constant and frequency exponent were extracted from the measured TEM data in the mine area. The extracted parameters are used for interpreting the detection result as a supplement. This reveals why the TEM data acquired in the area has a low resolution. It was found that the DC resistivity and time constant do not significantly change the results, however, the chargeability and frequency exponent have a significant effect. Because of these influences, the SVD method is more accurate than the conventional method in the apparent resistivity profile. The area of the copper mine is confined accurately based on the SVD inverted data. The conclusion has been verified by drill and is identical to the practical geological situation.展开更多
基金supported by the Joint Fund of Seismological Science(Grant No.U1839206)the National R&D Program on Monitoring,Early Warning and Prevention of Major Natural Disaster(Grant No.2017YFC1500301)+2 种基金supported by IGGCAS Research Start-up Funds(Grant No.E0515402)National Natural Science Foundation of China(Grant No.E1115401)supported by National Natural Science Foundation of China(Grant No.11971258).
文摘The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.
文摘A geophysical survey was conducted in the Kelle-Bidjocka village, Messondo subdivision, in the Centre Region, Cameroon. The data acquisition was made by combining Schlumberger profiling and electrical soundings along six (06) profiles of 1500 m in length for a total of 64 geoelectrical stations’ survey conducted through a variable mesh 100 m × 200 m, or 100 m × 300 m. The equipment used is the DC resistivimeter Syscal Junior 48 (Iris Instrument). Processing and modelling of field data are made by using the Res2Dinv, Qwseln and Surfer software. The investigation methods used are electrical resistivity (DC) and induced polarization (IP) methods. The analyses and interpretations have helped to highlight areas of weakness or conductive discontinuities (fractures, faults, shear zones, etc.) in Precambrian gneiss formations, sometimes undergoing weathering processes. They identify the weathering or mineralogical accumulation horizons, the most promising is a mineralization channel identified in the NE-SW direction. The highlighted mineralization is characterized by strong gradients of chargeability or polarization. Samples and other geological evidences observed in the area are used to associate the most polarizable structures with ferriferous formations. Weakly polarizable and particularly conductive backgrounds identified by the inverse pseudo-sections are thought to be sulphate minerals or groundwater targets for future hydrogeological studies.
基金supported by the National Natural Science Foundation of China (No. 41130418)the Strategic Leading Science and Technology Programme (Class B) of the Chinese Academy of Sciences (No. XDB10010400)
文摘Earth medium is not completely elastic, with its viscosity resulting in attenuation and dispersion of seismic waves. Most viscoelastic numerical simulations are based on the finite-difference and finite-element methods. Targeted at viscoelastic numerical modeling for multilayered media, the constant-Q acoustic wave equation is transformed into the corresponding wave integral representation with its Green's function accounting for viscoelastic coefficients. An efficient alternative for full-waveform solution to the integral equation is proposed in this article by extending conventional frequency-domain boundary element methods to viscoelastic media. The viscoelastic boundary element method enjoys a distinct characteristic of the explicit use of boundary continuity conditions of displacement and traction, leading to a semi-analytical solution with sufficient accuracy for simulating the viscoelastic effect across irregular interfaces. Numerical experiments to study the viscoelastic absorption of different Q values demonstrate the accuracy and applicability of the method.
基金The study is supported by the National Natural Science Foundation of China
文摘In this paper the authors propose the natural field induced polarization (IP) method and present the way to pick up IP effect. The relations between the object and anomaly are studied by taking field experiments as examples. The effectiveness and usability of the method are testified.
文摘To achieve the current study, geoelectrical surveys along six (06) profiles of 4 km long in a 100 m × 200 m grid defined according to the triangulation principle in the Ngoura area (Tindikala-Boutou villages) have been made through electrical sounding and profiling following Schlumberger array. The instrument is the resistivimeter Syscal Junior 48 (IRIS Instrument) which uses the electrical current. The data have been processed and modelled with Res2Dinv and Winsev softwares then interpolated with Surfer software. Electrical methods used are the Direct current (DC) and the Induced Polarization (IP). Interpretation and analyses of results from each investigation method highlight weak zones or conductive discontinuities and mineralized zones. Conductive zones have been identified as shear zones within granitic structures of the Precambrian basement, according to the geologic and tectonic background of the area. The structural trend of these shear zones is E-W. The mineralization within it is N-S and characterized by high values of chargeability, essentially in the eastern part of the area under study. This mineralization proves the presence of metalliferous or sulphide heaps disseminated in weathered quartz veins which cross shear zones. Also, the poor mineralization and conductive structures in shear zones characterize the groundwater zones. The intense activities of gold washers encountered in the area enable to link that mineralization to gold within quartz veins. The near surface gold mineralization is eluvial or alluvial, and in depth this mineralization is primary.
文摘A geophysical investigation of subsurface structures using the Syscal Junior 48 resistivity-meter was conducted in Ngoura subdivision (East Cameroon) following a combined geoelectrical direct current (DC) approach involving Resistivity and IP methods. This investigation was allowed to collect data on forty-five (45) profiling lines at three acquisition levels (AB = 100 m, MN = 10 m;AB = 200 m, MN = 20 m and AB = 500 m, MN = 50 m) and two electric panels L1 and L4, using the Schlumberger array. Processing, modeling and interpretation of data collected using the Winsev, Res2Dinv and Surfer software helped in highlighting a conductive and strongly mineralized discontinuity in granite formations, which lined up with the NE-SW Kadei tectonic line. It extends beyond 100 m depth over an average width of 600 m. The mineralization associated with this discontinuity is identified by a high concentration of disseminated metalliferous minerals in quartz or pegmatite veins. The mining reconnaissance works in the study area and those of several authors have characterized this anomaly to a lode gold quartz or large pegmatite. The results of this study correlate with geological and tectonic data for the region marked by NE-SW Kadei tectonic line. Therefore, they confirm the reliability of a geoelectrical DC investigation method combining Resistivity and IP to the identification of ore bodies.
文摘在网络世界中,网际互连协议(Internet Protocol,IP)承载网扮演着信息流动的骨干角色,其运行表现和服务水准至关重要。文章详细探讨了IP承载网中服务质量(Quality of Service,QoS)策略的设计与优化方法,重点涵盖分类服务与标记、流量控制与优先级队列管理以及链路和节点策略的具体实现。通过流量分析与预测、智能资源分配与调度、动态策略调整以及故障检测与容错机制,提出一系列提升网络性能和用户体验的有效措施。
基金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 Technology R&D Program in the 11th Five year Plan of China(No.2007BAQ00168-1-1)the National Natural Science Foundation of China(No. 41103052/D0309)the Shanxi Province Excellent Graduate Innovation Program(No. 20113038)
文摘During transient electromagnetic method (TEM) exploration of a copper mine, we detected the late-channel TEM signal reversal phenomenon (a voltage change from positive to negative) caused by the influence of the induced polarization (IP) effect, which affects the depth and precision of the TEM detection. The conventional inversion method is inefficient because it is difficult to process the data. In this paper, the Cole-Cole model is adopted to analyze the effect of Dc resistivity, chargeability, time constant, and frequency exponent on the TEM response in an homogeneous half space model. Singular Value Decomposition (SVD) is used to invert the measured TEM data, and the Dc resistivity, chargeability, time constant and frequency exponent were extracted from the measured TEM data in the mine area. The extracted parameters are used for interpreting the detection result as a supplement. This reveals why the TEM data acquired in the area has a low resolution. It was found that the DC resistivity and time constant do not significantly change the results, however, the chargeability and frequency exponent have a significant effect. Because of these influences, the SVD method is more accurate than the conventional method in the apparent resistivity profile. The area of the copper mine is confined accurately based on the SVD inverted data. The conclusion has been verified by drill and is identical to the practical geological situation.