The conventional finite-element（FE） method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such...The conventional finite-element（FE） method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such as dipping interfaces and rough topography. We present an adaptive FE method for 2.5D forward modeling of induced polarization（IP）. In the presented method, an unstructured triangulation mesh that allows for local mesh refinement and flexible description of arbitrary model geometries is used. Furthermore, the mesh refinement process is guided by dual error estimate weighting to bias the refinement towards elements that affect the solution at the receiver locations. After the final mesh is generated, the Jacobian matrix is used to obtain the IP response on 2D structure models. We validate the adaptive FE algorithm using a vertical contact model. The validation shows that the elements near the receivers are highly refined and the average relative error of the potentials converges to 0.4 % and 1.2 % for the IP response. This suggests that the numerical solution of the adaptive FE algorithm converges to an accurate solution with the refined mesh. Finally, the accuracy and flexibility of the adaptive FE procedure are also validated using more complex models.展开更多
The study of induced polarization （IP） information extraction from magnetotelluric （MT） sounding data is of great and practical significance to the exploitation of deep mineral, oil and gas resources. The linear i...The study of induced polarization （IP） information extraction from magnetotelluric （MT） sounding data is of great and practical significance to the exploitation of deep mineral, oil and gas resources. The linear inversion method, which has been given priority in previous research on the IP information extraction method, has three main problems as follows： 1） dependency on the initial model, 2） easily falling into the local minimum, and 3） serious non-uniqueness of solutions. Taking the nonlinearity and nonconvexity of IP information extraction into consideration, a two-stage CO-PSO minimum structure inversion method using compute unified distributed architecture （CUDA） is proposed. On one hand, a novel Cauchy oscillation particle swarm optimization （CO-PSO） algorithm is applied to extract nonlinear IP information from MT sounding data, which is implemented as a parallel algorithm within CUDA computing architecture; on the other hand, the impact of the polarizability on the observation data is strengthened by introducing a second stage inversion process, and the regularization parameter is applied in the fitness function of PSO algorithm to solve the problem of multi-solution in inversion. The inversion simulation results of polarization layers in different strata of various geoelectric models show that the smooth models of resistivity and IP parameters can be obtained by the proposed algorithm, the results of which are relatively stable and accurate. The experiment results added with noise indicate that this method is robust to Gaussian white noise. Compared with the traditional PSO and GA algorithm, the proposed algorithm has more efficiency and better inversion results.展开更多
Determination of rock mechanical parameters is the most important step in rock mass quality evaluation and has significant impacts on geotechnical engineering practice.Rock mass integrity coefficient(KV)is one of the ...Determination of rock mechanical parameters is the most important step in rock mass quality evaluation and has significant impacts on geotechnical engineering practice.Rock mass integrity coefficient(KV)is one of the most efficient parameters,which is conventionally determined from boreholes.Such approaches,however,are time-consuming and expensive,offer low data coverage of point measurements,require heavy equipment,and are hardly conducted in steep topographic sites.Hence,borehole approaches cannot assess the subsurface thoroughly for rock mass quality evaluation.Alternatively,use of geophysical methods is non-invasive,rapid and economical.The proposed geophysical approach makes useful empirical correlation between geophysical and geotechnical parameters.We evaluated the rock mass quality via integration between KV measured from the limited boreholes and inverted resistivity obtained from electrical resistivity tomography(ERT).The borehole-ERT correlation provided KV along various geophysical profiles for more detailed 2D/3D(two-/three-dimensional)mapping of rock mass quality.The subsurface was thoroughly evaluated for rock masses with different engineering qualities,including highly weathered rock,semi-weathered rock,and fresh rock.Furthermore,ERT was integrated with induced polarization(IP)to resolve the uncertainty caused by water/clay content.Our results show that the proposed method,compared with the conventional approaches,can reduce the ambiguities caused by inadequate data,and give more accurate insights into the subsurface for rock mass quality evaluation.展开更多
The investigated study area locates about 72 km from Jeddah City, Makkah District, Kingdom of Saudi Arabia. The study mainly aimed to define the most significant zones of possible mineralization and outline their subs...The investigated study area locates about 72 km from Jeddah City, Makkah District, Kingdom of Saudi Arabia. The study mainly aimed to define the most significant zones of possible mineralization and outline their subsurface parameters (location and strike) in the southeast of Jabal Al-Qashah. Several geophysical methods have been conducted to carry out the goal, including ground magnetic, self-potential (SP), and induced polarization (IP) methods. Integrating these methods aims to help delineate the possible mineralization in the study area. The magnetic survey was conducted along 17 profiles where these profiles were chosen to be perpendicular to the strike of the quartz shear zone. Self-potential was applied along with five profiles covering the study area. At the same time, induced polarization was used along with one profile located at the western side of the study area corresponding to some magnetic and SP profiles. The most interesting zones of mineralization were successfully determined by comparing the results of residual magnetic profile (3), SP profile (1), and IP profile, where geological structures control some mineralization.展开更多
基金financially supported by the National Natural Science Foundation of China(No.41204055,41164003,and 41104074)Opening Project(No.SMIL-2014-06) of Hubei Subsurface Multi-scale Imaging Lab(SMIL),China University of Geosciences(Wuhan)
文摘The conventional finite-element（FE） method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such as dipping interfaces and rough topography. We present an adaptive FE method for 2.5D forward modeling of induced polarization（IP）. In the presented method, an unstructured triangulation mesh that allows for local mesh refinement and flexible description of arbitrary model geometries is used. Furthermore, the mesh refinement process is guided by dual error estimate weighting to bias the refinement towards elements that affect the solution at the receiver locations. After the final mesh is generated, the Jacobian matrix is used to obtain the IP response on 2D structure models. We validate the adaptive FE algorithm using a vertical contact model. The validation shows that the elements near the receivers are highly refined and the average relative error of the potentials converges to 0.4 % and 1.2 % for the IP response. This suggests that the numerical solution of the adaptive FE algorithm converges to an accurate solution with the refined mesh. Finally, the accuracy and flexibility of the adaptive FE procedure are also validated using more complex models.
基金Projects(41604117,41204054)supported by the National Natural Science Foundation of ChinaProjects(20110490149,2015M580700)supported by the Research Fund for the Doctoral Program of Higher Education,China+1 种基金Project(2015zzts064)supported by the Fundamental Research Funds for the Central Universities,ChinaProject(16B147)supported by the Scientific Research Fund of Hunan Provincial Education Department,China
文摘The study of induced polarization （IP） information extraction from magnetotelluric （MT） sounding data is of great and practical significance to the exploitation of deep mineral, oil and gas resources. The linear inversion method, which has been given priority in previous research on the IP information extraction method, has three main problems as follows： 1） dependency on the initial model, 2） easily falling into the local minimum, and 3） serious non-uniqueness of solutions. Taking the nonlinearity and nonconvexity of IP information extraction into consideration, a two-stage CO-PSO minimum structure inversion method using compute unified distributed architecture （CUDA） is proposed. On one hand, a novel Cauchy oscillation particle swarm optimization （CO-PSO） algorithm is applied to extract nonlinear IP information from MT sounding data, which is implemented as a parallel algorithm within CUDA computing architecture; on the other hand, the impact of the polarizability on the observation data is strengthened by introducing a second stage inversion process, and the regularization parameter is applied in the fitness function of PSO algorithm to solve the problem of multi-solution in inversion. The inversion simulation results of polarization layers in different strata of various geoelectric models show that the smooth models of resistivity and IP parameters can be obtained by the proposed algorithm, the results of which are relatively stable and accurate. The experiment results added with noise indicate that this method is robust to Gaussian white noise. Compared with the traditional PSO and GA algorithm, the proposed algorithm has more efficiency and better inversion results.
基金supported by Xinjiang Key Laboratory of Geohazards Prevention(Grant No.XKLGP2022K07)Key R&D Program of Xinjiang Uygur Autonomous Region(Grant No.2022B03001-2)the Third Xinjiang Scientific Expedition Program(Grant No.2022xjkk1305).
文摘Determination of rock mechanical parameters is the most important step in rock mass quality evaluation and has significant impacts on geotechnical engineering practice.Rock mass integrity coefficient(KV)is one of the most efficient parameters,which is conventionally determined from boreholes.Such approaches,however,are time-consuming and expensive,offer low data coverage of point measurements,require heavy equipment,and are hardly conducted in steep topographic sites.Hence,borehole approaches cannot assess the subsurface thoroughly for rock mass quality evaluation.Alternatively,use of geophysical methods is non-invasive,rapid and economical.The proposed geophysical approach makes useful empirical correlation between geophysical and geotechnical parameters.We evaluated the rock mass quality via integration between KV measured from the limited boreholes and inverted resistivity obtained from electrical resistivity tomography(ERT).The borehole-ERT correlation provided KV along various geophysical profiles for more detailed 2D/3D(two-/three-dimensional)mapping of rock mass quality.The subsurface was thoroughly evaluated for rock masses with different engineering qualities,including highly weathered rock,semi-weathered rock,and fresh rock.Furthermore,ERT was integrated with induced polarization(IP)to resolve the uncertainty caused by water/clay content.Our results show that the proposed method,compared with the conventional approaches,can reduce the ambiguities caused by inadequate data,and give more accurate insights into the subsurface for rock mass quality evaluation.
文摘The investigated study area locates about 72 km from Jeddah City, Makkah District, Kingdom of Saudi Arabia. The study mainly aimed to define the most significant zones of possible mineralization and outline their subsurface parameters (location and strike) in the southeast of Jabal Al-Qashah. Several geophysical methods have been conducted to carry out the goal, including ground magnetic, self-potential (SP), and induced polarization (IP) methods. Integrating these methods aims to help delineate the possible mineralization in the study area. The magnetic survey was conducted along 17 profiles where these profiles were chosen to be perpendicular to the strike of the quartz shear zone. Self-potential was applied along with five profiles covering the study area. At the same time, induced polarization was used along with one profile located at the western side of the study area corresponding to some magnetic and SP profiles. The most interesting zones of mineralization were successfully determined by comparing the results of residual magnetic profile (3), SP profile (1), and IP profile, where geological structures control some mineralization.
