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.展开更多
Temperature gradient and cooling rate have an obvious effect on formation of methane hydrate. The process for formation of methane hydrate in coarse sand is monitored to tmderstand the relationship between temperature...Temperature gradient and cooling rate have an obvious effect on formation of methane hydrate. The process for formation of methane hydrate in coarse sand is monitored to tmderstand the relationship between temperature gradient and cooling rate and nucleation, growth and distribution of methane hydrate by using the electrical resistivity method. The results show that the change of resistivity can better reflect the nucleation and growth and distribution of methane hydrate. Temperature gradient promotes the nucleation, formation, and formation rate of methane hydrate. At a temperature gradient of 0.11℃/cm, the rate of methane hydrate formation and saturation reaches a maximum. Cooling rate has little effect on the methane hydrate formation process. Judging from the outcome of final spatial distribution of methane hydrate, the cooling rate has an obvious but irregular effect in coarse sand. The effect of tempera^re gradient on distribution of methane hydrate in coarse sand is less than that of cooling rate. At a temperature gradient of 0.07℃/cm, methane hydrate is distributed uniformly in the sample. If the temperature gradient is higher or lower than this value, the hydrate is enriched in the upper layer of sample.展开更多
Multi-electrodes Resistivity Imaging Survey(MRIS)is an array method of electrical survey.In practice how to choose a reasonable array is the key to get reliable survey results.Based on four methods of MRIS such as Wen...Multi-electrodes Resistivity Imaging Survey(MRIS)is an array method of electrical survey.In practice how to choose a reasonable array is the key to get reliable survey results.Based on four methods of MRIS such as Wenner,Schlumberger,Pole-pole and Dipole-dipole the authors established the model,by studying the result of the forward numerical simulation modeling and inverse modeling,and analyzed the differences among the different forms of detection devices.展开更多
Taking into account the quantum size effects and considering three types of scattering from bulk impurities,rough surface and rough interfaces, we use quantum-statistical Green's function approach and Kubo theory ...Taking into account the quantum size effects and considering three types of scattering from bulk impurities,rough surface and rough interfaces, we use quantum-statistical Green's function approach and Kubo theory to calculate the electronic conductivity and the giant magnetoresistance in magnetic multilayered cylindrical systems. It is found that in the limit of weakly scattering from impurities surface and interfaces, the total conductivity is given by a sum of conductivities of all the subbands and two spin-channels. For each subband and each spin-channel the scattering rate due to the impurities, surface and interfaces is added up.展开更多
Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid...Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.展开更多
基金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.
基金supported by the Chinese Academy of Sciences Action-plan for Western Project(No.KZCX2-XB3-03)the National Natural Science Foundation of China(No.41001038,51266005)the National Natural Science Foundation of China(No.41101070,1106ZBB007)
文摘Temperature gradient and cooling rate have an obvious effect on formation of methane hydrate. The process for formation of methane hydrate in coarse sand is monitored to tmderstand the relationship between temperature gradient and cooling rate and nucleation, growth and distribution of methane hydrate by using the electrical resistivity method. The results show that the change of resistivity can better reflect the nucleation and growth and distribution of methane hydrate. Temperature gradient promotes the nucleation, formation, and formation rate of methane hydrate. At a temperature gradient of 0.11℃/cm, the rate of methane hydrate formation and saturation reaches a maximum. Cooling rate has little effect on the methane hydrate formation process. Judging from the outcome of final spatial distribution of methane hydrate, the cooling rate has an obvious but irregular effect in coarse sand. The effect of tempera^re gradient on distribution of methane hydrate in coarse sand is less than that of cooling rate. At a temperature gradient of 0.07℃/cm, methane hydrate is distributed uniformly in the sample. If the temperature gradient is higher or lower than this value, the hydrate is enriched in the upper layer of sample.
基金Supported by Project of the National High Technology Research and Development Program of China(No.2007AA06Z215)
文摘Multi-electrodes Resistivity Imaging Survey(MRIS)is an array method of electrical survey.In practice how to choose a reasonable array is the key to get reliable survey results.Based on four methods of MRIS such as Wenner,Schlumberger,Pole-pole and Dipole-dipole the authors established the model,by studying the result of the forward numerical simulation modeling and inverse modeling,and analyzed the differences among the different forms of detection devices.
文摘Taking into account the quantum size effects and considering three types of scattering from bulk impurities,rough surface and rough interfaces, we use quantum-statistical Green's function approach and Kubo theory to calculate the electronic conductivity and the giant magnetoresistance in magnetic multilayered cylindrical systems. It is found that in the limit of weakly scattering from impurities surface and interfaces, the total conductivity is given by a sum of conductivities of all the subbands and two spin-channels. For each subband and each spin-channel the scattering rate due to the impurities, surface and interfaces is added up.
基金Projects(2006AA06Z105, 2007AA06Z134) supported by the National High-Tech Research and Development Program of ChinaProjects(2007, 2008) supported by China Scholarship Council (CSC)
文摘Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.