The disaster of seawater intrusion seriously affects people's lives and restricts economic development,so the detection and treatment of seawater intrusion is a long-term task.On the basis of field investigation a...The disaster of seawater intrusion seriously affects people's lives and restricts economic development,so the detection and treatment of seawater intrusion is a long-term task.On the basis of field investigation and water quality analysis,according to the change characteristics of apparent resistivity of groundwater after Cl-reaches 250 mg/L,the theoretical basis for the application of high-density resistivity method was determined,and the characteristic values of apparent resistivity for seawater intrusion interfaces in different geological characteristic regions in Laizhou Bay area were determined by typical profile tests.Combined with water quality investigation and other means,profiles for the high-density resistivity method were arranged,and the interfaces between saline and fresh water were accurately divided.展开更多
The authors employ the high-density resistivity method to image the subsurface structure of a mountain in Erdaojiang District,Tonghua City,Jilin Province,China,to evaluate the potential risk of slope failure on surrou...The authors employ the high-density resistivity method to image the subsurface structure of a mountain in Erdaojiang District,Tonghua City,Jilin Province,China,to evaluate the potential risk of slope failure on surrounding residential areas and infrastructure,and identify a shallow fault that extends across the center of the mountain and is perpendicular to the mountain slope and accurately locate the spatial position and depth of another fault on the southern side of the mountain.The results provide an important basis for evaluating mountain slope stability.This study also demonstrates that the high-density resistivity method is effective for detecting mountain faults.展开更多
By determining the distribution and extent of geological structures surrounding the Mingyan Tunnel,Xicheng Town,Helong City,Jilin Province,we can evaluate the stability of the rock mass and assess potential hazards du...By determining the distribution and extent of geological structures surrounding the Mingyan Tunnel,Xicheng Town,Helong City,Jilin Province,we can evaluate the stability of the rock mass and assess potential hazards during tunnel construction.We use the high-density resistivity method to analyze the subsurface structure of the study area.Conductive anomalies are likely to represent joint and fissure systems within strongly weathered host rocks,and the bedrock surrounding the tunnel is relatively stable and does not contain well-developed faults.High-density resistivity analysis can provide valuable information in the context of tunnel engineering and safety.展开更多
Some unfavorable geological conditions can affect the construction of tunnels.In order to evaluate the damage degree of tunnel construction and determine the surrounding rock grade and stability of the tunnel,the auth...Some unfavorable geological conditions can affect the construction of tunnels.In order to evaluate the damage degree of tunnel construction and determine the surrounding rock grade and stability of the tunnel,the authors used high-density resistivity method to detect the surrounding rocks of Shimodong tunnel in Xicheng Town of Helong City.The underground resistivity structures of the entrance,exit and middle parts of the tunnel are obtained.Through analysis,it is found that there are no bedrock faults near the tunnel,although some joints and fissures are developed in some locations,which are characterized by low-resistivity anomalies.The tunnel structures are stable overall,favorable for safe and efficient construction.The study also proves the good application effect of the high-density resistivity method in tunnel safety detection.展开更多
The authors employ the high-density resistivity method during an archaeological investigation of Sumicheng site,an ancient city of the Tang Dynasty,to find evidence of human activities and locate a favorable target fo...The authors employ the high-density resistivity method during an archaeological investigation of Sumicheng site,an ancient city of the Tang Dynasty,to find evidence of human activities and locate a favorable target for archaeological excavation in the southern part of the outer city.There are two obvious high-resistivity structures,the south wall of the inner city and an ancient building near the south gate along the outer city wall,of which the resistivities are indicative of rammed soil foundations.The south wall of the inner city is continuous but is cut off abruptly to the east,which we suggest it is due to either wall damage or destruction.The resistivity signature of the target area is verified by archaeological excavation,proving the feasibility and effectiveness of implementing the high-density resistivity method for archaeological exploration.展开更多
High-density resistivity method is a new, efficient electrical prospecting method, which can complete a two-dimensional (vertical and horizontal) prospecting process, possesses certain imaging functions for the geo-...High-density resistivity method is a new, efficient electrical prospecting method, which can complete a two-dimensional (vertical and horizontal) prospecting process, possesses certain imaging functions for the geo-electric structure, and integrates electric profiling method with electric sounding method together. In this paper, the basic principle, data processing, and result explanation and inference of high- density resistivity method are introduced by taking the application of high-density resistivity method to the prospecting project in the slope of Gongchangling Open Pit, Liaoyang. The result of the prospecting result map analysis showed that the prospecting result was basically in line with the actual situation and proved the great significance of high-density resistivity method to the evaluation on the slope stability of Gongchangling open pit.展开更多
The measurement of resistivity in a compressed material within a diamond anvil cell presents significant challenges.The high-pressure exper-imental setup makes it difficult to directly measure the size changes induced...The measurement of resistivity in a compressed material within a diamond anvil cell presents significant challenges.The high-pressure exper-imental setup makes it difficult to directly measure the size changes induced by pressure in the three crystallographic directions of the sample.In this study,we introduce a novel and effective method that addresses these technical challenges.This method is anticipated to offer a valuable foundation for high-pressure investigations on quantum materials,particularly those with anisotropic layered structures.展开更多
The electrical resistivity method is a geophysical tool used to characterize the subsoil and can provide an important information for precision agriculture. The lack of knowledge about agronomic properties of the soil...The electrical resistivity method is a geophysical tool used to characterize the subsoil and can provide an important information for precision agriculture. The lack of knowledge about agronomic properties of the soil tends to affect the agricultural coffee production system. Therefore, research related to geoelectrical properties of soil such as resistivity for characterization the region of the study for coffee cultivation purposes can improve and optimize the production. This resistivity method allows to investigate the subsurface through different techniques: 1D vertical electrical sounding and electrical imaging. The acquisition of data using these techniques permitted the creation of 2D resistivity cross section from the study area. The geoelectrical data was acquired by using a resistivity meter equipment and was processed in different softwares. The results of the geoelectrical characterization from 1D resistivity model and 2D resistivity electrical sections show that in the study area of Kabiri, there are 8 varieties of geoelectrical layers with different resistivity or conductivity. Near survey in the study area, the lowest resistivity is around 0.322 Ω·m, while the highest is about 92.1 Ω·m. These values illustrated where is possible to plant coffee for suggestion of specific fertilization plan for some area to improve the cultivation.展开更多
The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy ...The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.展开更多
In this paper,we construct a high-order discontinuous Galerkin(DG)method which can preserve the positivity of the density and the pressure for the viscous and resistive magnetohydrodynamics(VRMHD).To control the diver...In this paper,we construct a high-order discontinuous Galerkin(DG)method which can preserve the positivity of the density and the pressure for the viscous and resistive magnetohydrodynamics(VRMHD).To control the divergence error in the magnetic field,both the local divergence-free basis and the Godunov source term would be employed for the multi-dimensional VRMHD.Rigorous theoretical analyses are presented for one-dimensional and multi-dimensional DG schemes,respectively,showing that the scheme can maintain the positivity-preserving(PP)property under some CFL conditions when combined with the strong-stability-preserving time discretization.Then,general frameworks are established to construct the PP limiter for arbitrary order of accuracy DG schemes.Numerical tests demonstrate the effectiveness of the proposed schemes.展开更多
To speed up three-dimensional (3D) DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method (AGMG). We first discretize the differential equation of the secondar...To speed up three-dimensional (3D) DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method (AGMG). We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner (AGMG-CG) to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling (3DDCXH). In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.展开更多
Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell's equations for an anisotropic mo...Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell's equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic (EM) fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.展开更多
High-density electrical method has been proved to be an effective method for probing shallow sedimentary layers.It is principally used to identify the boundary between the Quaternary soil layer and bedrock according t...High-density electrical method has been proved to be an effective method for probing shallow sedimentary layers.It is principally used to identify the boundary between the Quaternary soil layer and bedrock according to the vertical change of apparent resistivity.However,the artificial filling layer has the characteristics of heterogeneity and high porosity,which makes it challenging to detect the artificial filling layer by high-density electrical method.The key to solve this problem is to detect the difference of conductivity between the filling layer and the underlying bedrock.This paper takes the land in Chengjiangshan area of Huaibei City,Anhui Province as the detection target.On the basis of fully analyzing the physical properties of the artificial filling layer,two-dimensional high-density electrical survey and inversion are used to define the thickness of the artificial filling layer.The research shows that the highdensity resistivity method has obvious advantages in delineating the distribution of bedrock and the thickness of the filling layer,and the reliability of the high-density electrical method in the detection of the artificial filling layer,and delineates the scope of the filling layer is verified by the borehole data.展开更多
The research about subsurface characteristics by using transient electromagnetic method(TEM) and high density resistivity method(HDRM) were already conducted in Ordos. The objective of this research is to detect c...The research about subsurface characteristics by using transient electromagnetic method(TEM) and high density resistivity method(HDRM) were already conducted in Ordos. The objective of this research is to detect coalmine goaf areas based on rock resistivity. The data processing using wavelet transform, three point smoothing, RES2 DINV and Maxwell processing software to obtain 2D resistivity structure. The results showed that the layers with maximum resistivity values(30e33 U m on Line 1, 30e31 U m on Line 2, 32e40 U m on Line3) are founded at station 1e7, and 14e20 on Line 1,13e18 on Line 2, and 8e13 and 16e20 on Line 3 which is predicted as goaf layer, and the minimum resistivity values(20e26 U m of TEM, 45e75 U m of HDRM) at the other layers. This resistivity difference was caused by the geology and characteristics of the study area which is located close by the cleugh with rich coal, so the goaf area distinguishable with aquifer layer and coal seam. The results were also significant accidents and serious destruction of ecological environment.展开更多
Electrical measurement was employed to investigate the early hydration characteristics of cement pastes with different dosages of superplasticizer in the same W/C ratio. The hyperbolic method was applied to analyze th...Electrical measurement was employed to investigate the early hydration characteristics of cement pastes with different dosages of superplasticizer in the same W/C ratio. The hyperbolic method was applied to analyze the electrical resistivity development. The peak point (Ph) on the hyperbolic curve could be easily read. The time (th) to reach the point Ph had strong relations with the setting time. th was delayed with the increment of the dosage of superplasticizer. The time th was used to plot the relationship between the initial setting time and final setting time. The hyperbolic equation was established to predict the ultimate resistivity. The retardation effect of the superplasticizer was confirmed in the same W/C ratio by setting time and isothermal heat evolution.展开更多
The objective of this work is to develop a novel methodology for determining real resistivity of conductive asphalt concrete based on two-electrode method.Due to an influence of contact resistance,the measured resisti...The objective of this work is to develop a novel methodology for determining real resistivity of conductive asphalt concrete based on two-electrode method.Due to an influence of contact resistance,the measured resistivity is always not equal to the real resistivity.To determine the real resistivity,a linear relationship of the measured resistivity,contact resistance and the real resistivity was established.Then experiments for six specimens with varying graphite contents were designed and performed to validate the formulation.Results of experiments demonstrate that the slope of the line represents contact resistance,and the intercept indicates the real resistivity.The effects of graphite content on contact resistance and real resistivity are also revealed.Finally,results show that the influence of contact resistance on accuracy of resisitvity measurement becomes more serious if graphite content is beyond 3%.Hence,it is the time to choose this novel methodology to determine the real resistivity of asphalt concrete by taking account of contact resistance.展开更多
The detecting examples using the high density resistivity method, about the evaluation of containing water characteristic from the floor rock and the height of overburden failure, were given. It expresses that the hig...The detecting examples using the high density resistivity method, about the evaluation of containing water characteristic from the floor rock and the height of overburden failure, were given. It expresses that the high density resistivity method has good effect with strong maneuverability and continuous observing section during the prevention and cure for mine water disaster. At the same time, the article pointed out that the study of space data inversion and dynamic watching technology about the high density resistivity method must be enhanced in the future because of special condition of data collecting in mine.展开更多
The compressive strength development of Portland cement pastes was investigated by the electrical resistivity method and the maturity method.The experiments were carried out on the cement pastes with different water-c...The compressive strength development of Portland cement pastes was investigated by the electrical resistivity method and the maturity method.The experiments were carried out on the cement pastes with different water-cement ratios at different curing temperatures.The results show that the application of the maturity method has limitation to obtain the strength.It is found that both of the compressive strength and the electrical resistivity follow hyperbolic trend for all the mixes.The hyperbolic equation of each mix is obtained to estimate the ultimate resistivity value which can probably be reached.The relationship between electrical resistivity and compressive strength of the cement pastes is established based on the test results and interpreted by the empirical Archie equation and a strength-porosity equation.The relationship between the electrical resistivity after temperature correction and the compressive strength was linear and independent of curing temperature and water-cement ratio.展开更多
In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position inf...In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position information of measuring point.Utilizing very fast simulated annealing(VFSA) to fit the theoretical electromagnetic force(EMF) and measured EMF could obtain the all-time apparent resistivity of the measuring points in rectangular transmitting loop.The selective cope of initial model of VFSA could be confirmed by taking the late time apparent resistivity of transient electromagnetic method as the prior information.For verifying the correctness,the all-time apparent resistivities of the geoelectric models were calculated by VFSA and dichotomy,respectively.The results indicate that the relative differences of apparent resistivities calculated by these two methods are within 3%.The change of measuring point position has little influence on the tracing pattern of all-time apparent resistivity.The first branch of the curve of all-time apparent resistivity is close to the resistivity of the first layer medium and the last branch is close to the resistivity of the last layer medium,which proves the correctness of the arithmetics proposed.展开更多
The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and st...The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and stored in two parts separately. One part is associated with the volume integral and the other is associated with the subsurface boundary integral. The equivalent multiple linear systems with closer right-hand sides than the original systems were constructed. A recycling Krylov subspace technique was employed to solve the multiple linear systems. The solution of the seed system was used as an initial guess for the subsequent systems. The results of two numerical experiments show that the improved algorithm reduces the iterations and CPU time by almost 50%, compared with the classical preconditioned conjugate gradient method.展开更多
文摘The disaster of seawater intrusion seriously affects people's lives and restricts economic development,so the detection and treatment of seawater intrusion is a long-term task.On the basis of field investigation and water quality analysis,according to the change characteristics of apparent resistivity of groundwater after Cl-reaches 250 mg/L,the theoretical basis for the application of high-density resistivity method was determined,and the characteristic values of apparent resistivity for seawater intrusion interfaces in different geological characteristic regions in Laizhou Bay area were determined by typical profile tests.Combined with water quality investigation and other means,profiles for the high-density resistivity method were arranged,and the interfaces between saline and fresh water were accurately divided.
基金Supported by National Key R&D Program of China and Fundamental Research Funds for the Central Universities(2017YFC0601305)。
文摘The authors employ the high-density resistivity method to image the subsurface structure of a mountain in Erdaojiang District,Tonghua City,Jilin Province,China,to evaluate the potential risk of slope failure on surrounding residential areas and infrastructure,and identify a shallow fault that extends across the center of the mountain and is perpendicular to the mountain slope and accurately locate the spatial position and depth of another fault on the southern side of the mountain.The results provide an important basis for evaluating mountain slope stability.This study also demonstrates that the high-density resistivity method is effective for detecting mountain faults.
基金Supported by The National Natural Science Foundation of China(41504076)Jilin Science and Technological Development Program(20180101093JC)。
文摘By determining the distribution and extent of geological structures surrounding the Mingyan Tunnel,Xicheng Town,Helong City,Jilin Province,we can evaluate the stability of the rock mass and assess potential hazards during tunnel construction.We use the high-density resistivity method to analyze the subsurface structure of the study area.Conductive anomalies are likely to represent joint and fissure systems within strongly weathered host rocks,and the bedrock surrounding the tunnel is relatively stable and does not contain well-developed faults.High-density resistivity analysis can provide valuable information in the context of tunnel engineering and safety.
基金National Key R&D Program of China(2017YFC0601305)Fundamental Research Funds for the Central Universities.
文摘Some unfavorable geological conditions can affect the construction of tunnels.In order to evaluate the damage degree of tunnel construction and determine the surrounding rock grade and stability of the tunnel,the authors used high-density resistivity method to detect the surrounding rocks of Shimodong tunnel in Xicheng Town of Helong City.The underground resistivity structures of the entrance,exit and middle parts of the tunnel are obtained.Through analysis,it is found that there are no bedrock faults near the tunnel,although some joints and fissures are developed in some locations,which are characterized by low-resistivity anomalies.The tunnel structures are stable overall,favorable for safe and efficient construction.The study also proves the good application effect of the high-density resistivity method in tunnel safety detection.
基金projects of Jilin Provincial Institute of Cultural Relics and Archaeology(No.3S318B564423)National Key R&D Program of China(No.2017YFC0601305)Fundamental Research Funds for the Central Universities.
文摘The authors employ the high-density resistivity method during an archaeological investigation of Sumicheng site,an ancient city of the Tang Dynasty,to find evidence of human activities and locate a favorable target for archaeological excavation in the southern part of the outer city.There are two obvious high-resistivity structures,the south wall of the inner city and an ancient building near the south gate along the outer city wall,of which the resistivities are indicative of rammed soil foundations.The south wall of the inner city is continuous but is cut off abruptly to the east,which we suggest it is due to either wall damage or destruction.The resistivity signature of the target area is verified by archaeological excavation,proving the feasibility and effectiveness of implementing the high-density resistivity method for archaeological exploration.
文摘High-density resistivity method is a new, efficient electrical prospecting method, which can complete a two-dimensional (vertical and horizontal) prospecting process, possesses certain imaging functions for the geo-electric structure, and integrates electric profiling method with electric sounding method together. In this paper, the basic principle, data processing, and result explanation and inference of high- density resistivity method are introduced by taking the application of high-density resistivity method to the prospecting project in the slope of Gongchangling Open Pit, Liaoyang. The result of the prospecting result map analysis showed that the prospecting result was basically in line with the actual situation and proved the great significance of high-density resistivity method to the evaluation on the slope stability of Gongchangling open pit.
基金supported by the National Key Research and Development Program of China(Grant Nos.2022YFA1403900 and 2021YFA1401800)the NSF of China(Grant Nos.U2032214 and 12104487).
文摘The measurement of resistivity in a compressed material within a diamond anvil cell presents significant challenges.The high-pressure exper-imental setup makes it difficult to directly measure the size changes induced by pressure in the three crystallographic directions of the sample.In this study,we introduce a novel and effective method that addresses these technical challenges.This method is anticipated to offer a valuable foundation for high-pressure investigations on quantum materials,particularly those with anisotropic layered structures.
文摘The electrical resistivity method is a geophysical tool used to characterize the subsoil and can provide an important information for precision agriculture. The lack of knowledge about agronomic properties of the soil tends to affect the agricultural coffee production system. Therefore, research related to geoelectrical properties of soil such as resistivity for characterization the region of the study for coffee cultivation purposes can improve and optimize the production. This resistivity method allows to investigate the subsurface through different techniques: 1D vertical electrical sounding and electrical imaging. The acquisition of data using these techniques permitted the creation of 2D resistivity cross section from the study area. The geoelectrical data was acquired by using a resistivity meter equipment and was processed in different softwares. The results of the geoelectrical characterization from 1D resistivity model and 2D resistivity electrical sections show that in the study area of Kabiri, there are 8 varieties of geoelectrical layers with different resistivity or conductivity. Near survey in the study area, the lowest resistivity is around 0.322 Ω·m, while the highest is about 92.1 Ω·m. These values illustrated where is possible to plant coffee for suggestion of specific fertilization plan for some area to improve the cultivation.
基金Projects(41074085,41374118)supported by the National Natural Science Foundation of ChinaProject(20120162110015)supported by Doctoral Fund of Ministry of Education of ChinaProject(NCET-12-0551)supported by Program for New Century Excellent Talents in University,China
文摘The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.
基金supported by the NSFC Grant 11901555,12271499the Cyrus Tang Foundationsupported by the NSFC Grant 11871448 and 12126604.
文摘In this paper,we construct a high-order discontinuous Galerkin(DG)method which can preserve the positivity of the density and the pressure for the viscous and resistive magnetohydrodynamics(VRMHD).To control the divergence error in the magnetic field,both the local divergence-free basis and the Godunov source term would be employed for the multi-dimensional VRMHD.Rigorous theoretical analyses are presented for one-dimensional and multi-dimensional DG schemes,respectively,showing that the scheme can maintain the positivity-preserving(PP)property under some CFL conditions when combined with the strong-stability-preserving time discretization.Then,general frameworks are established to construct the PP limiter for arbitrary order of accuracy DG schemes.Numerical tests demonstrate the effectiveness of the proposed schemes.
基金supported by the Natural Science Foundation of China(Nos.41404057,41674077 and 411640034)the Nuclear Energy Development Project of China,and the‘555’Project of Gan Po Excellent People
文摘To speed up three-dimensional (3D) DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method (AGMG). We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner (AGMG-CG) to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling (3DDCXH). In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.
基金financially supported by the National Hi-tech Research and Development Program of China(863 Program)(No.2012AA09A20103)
文摘Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell's equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic (EM) fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.
文摘High-density electrical method has been proved to be an effective method for probing shallow sedimentary layers.It is principally used to identify the boundary between the Quaternary soil layer and bedrock according to the vertical change of apparent resistivity.However,the artificial filling layer has the characteristics of heterogeneity and high porosity,which makes it challenging to detect the artificial filling layer by high-density electrical method.The key to solve this problem is to detect the difference of conductivity between the filling layer and the underlying bedrock.This paper takes the land in Chengjiangshan area of Huaibei City,Anhui Province as the detection target.On the basis of fully analyzing the physical properties of the artificial filling layer,two-dimensional high-density electrical survey and inversion are used to define the thickness of the artificial filling layer.The research shows that the highdensity resistivity method has obvious advantages in delineating the distribution of bedrock and the thickness of the filling layer,and the reliability of the high-density electrical method in the detection of the artificial filling layer,and delineates the scope of the filling layer is verified by the borehole data.
基金supported by the Institute of Seismology Foundation, China Earthquake Administration (201326126)
文摘The research about subsurface characteristics by using transient electromagnetic method(TEM) and high density resistivity method(HDRM) were already conducted in Ordos. The objective of this research is to detect coalmine goaf areas based on rock resistivity. The data processing using wavelet transform, three point smoothing, RES2 DINV and Maxwell processing software to obtain 2D resistivity structure. The results showed that the layers with maximum resistivity values(30e33 U m on Line 1, 30e31 U m on Line 2, 32e40 U m on Line3) are founded at station 1e7, and 14e20 on Line 1,13e18 on Line 2, and 8e13 and 16e20 on Line 3 which is predicted as goaf layer, and the minimum resistivity values(20e26 U m of TEM, 45e75 U m of HDRM) at the other layers. This resistivity difference was caused by the geology and characteristics of the study area which is located close by the cleugh with rich coal, so the goaf area distinguishable with aquifer layer and coal seam. The results were also significant accidents and serious destruction of ecological environment.
基金the National Natural Science Foundation of China(No.50778078)
文摘Electrical measurement was employed to investigate the early hydration characteristics of cement pastes with different dosages of superplasticizer in the same W/C ratio. The hyperbolic method was applied to analyze the electrical resistivity development. The peak point (Ph) on the hyperbolic curve could be easily read. The time (th) to reach the point Ph had strong relations with the setting time. th was delayed with the increment of the dosage of superplasticizer. The time th was used to plot the relationship between the initial setting time and final setting time. The hyperbolic equation was established to predict the ultimate resistivity. The retardation effect of the superplasticizer was confirmed in the same W/C ratio by setting time and isothermal heat evolution.
基金Project(51178348)supported by the National Natural Science Foundation of China
文摘The objective of this work is to develop a novel methodology for determining real resistivity of conductive asphalt concrete based on two-electrode method.Due to an influence of contact resistance,the measured resistivity is always not equal to the real resistivity.To determine the real resistivity,a linear relationship of the measured resistivity,contact resistance and the real resistivity was established.Then experiments for six specimens with varying graphite contents were designed and performed to validate the formulation.Results of experiments demonstrate that the slope of the line represents contact resistance,and the intercept indicates the real resistivity.The effects of graphite content on contact resistance and real resistivity are also revealed.Finally,results show that the influence of contact resistance on accuracy of resisitvity measurement becomes more serious if graphite content is beyond 3%.Hence,it is the time to choose this novel methodology to determine the real resistivity of asphalt concrete by taking account of contact resistance.
文摘The detecting examples using the high density resistivity method, about the evaluation of containing water characteristic from the floor rock and the height of overburden failure, were given. It expresses that the high density resistivity method has good effect with strong maneuverability and continuous observing section during the prevention and cure for mine water disaster. At the same time, the article pointed out that the study of space data inversion and dynamic watching technology about the high density resistivity method must be enhanced in the future because of special condition of data collecting in mine.
基金Funding by the National Natural Science Foundation of China (Nos.50778078 and 51178202)the Doctoral Research Fund from Wuhan Institute of Technology
文摘The compressive strength development of Portland cement pastes was investigated by the electrical resistivity method and the maturity method.The experiments were carried out on the cement pastes with different water-cement ratios at different curing temperatures.The results show that the application of the maturity method has limitation to obtain the strength.It is found that both of the compressive strength and the electrical resistivity follow hyperbolic trend for all the mixes.The hyperbolic equation of each mix is obtained to estimate the ultimate resistivity value which can probably be reached.The relationship between electrical resistivity and compressive strength of the cement pastes is established based on the test results and interpreted by the empirical Archie equation and a strength-porosity equation.The relationship between the electrical resistivity after temperature correction and the compressive strength was linear and independent of curing temperature and water-cement ratio.
基金Projects(40804027,41074085) supported by the National Natural Science Foundation of ChinaProject(09JJ3048) supported by the Natural Science Foundation of Hunan Province,ChinaProject(200805331082) supported by the Research Fund for the Doctoral Program of Higher Education,China
文摘In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position information of measuring point.Utilizing very fast simulated annealing(VFSA) to fit the theoretical electromagnetic force(EMF) and measured EMF could obtain the all-time apparent resistivity of the measuring points in rectangular transmitting loop.The selective cope of initial model of VFSA could be confirmed by taking the late time apparent resistivity of transient electromagnetic method as the prior information.For verifying the correctness,the all-time apparent resistivities of the geoelectric models were calculated by VFSA and dichotomy,respectively.The results indicate that the relative differences of apparent resistivities calculated by these two methods are within 3%.The change of measuring point position has little influence on the tracing pattern of all-time apparent resistivity.The first branch of the curve of all-time apparent resistivity is close to the resistivity of the first layer medium and the last branch is close to the resistivity of the last layer medium,which proves the correctness of the arithmetics proposed.
基金Projects(40974077,41164004)supported by the National Natural Science Foundation of ChinaProject(2007AA06Z134)supported by the National High Technology Research and Development Program of China+2 种基金Projects(2011GXNSFA018003,0832263)supported by the Natural Science Foundation of Guangxi Province,ChinaProject supported by Program for Excellent Talents in Guangxi Higher Education Institution,ChinaProject supported by the Foundation of Guilin University of Technology,China
文摘The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and stored in two parts separately. One part is associated with the volume integral and the other is associated with the subsurface boundary integral. The equivalent multiple linear systems with closer right-hand sides than the original systems were constructed. A recycling Krylov subspace technique was employed to solve the multiple linear systems. The solution of the seed system was used as an initial guess for the subsequent systems. The results of two numerical experiments show that the improved algorithm reduces the iterations and CPU time by almost 50%, compared with the classical preconditioned conjugate gradient method.