砂土液化是常见的地震灾害,目前应用于研究砂土液化动力特性的室内试验以及模型试验还不能全面反映土体液化全过程。计算流体动力学(computational fluid dynamics,CFD)与离散元法(discrete element method,DEM)耦合模拟方法能够准确地...砂土液化是常见的地震灾害,目前应用于研究砂土液化动力特性的室内试验以及模型试验还不能全面反映土体液化全过程。计算流体动力学(computational fluid dynamics,CFD)与离散元法(discrete element method,DEM)耦合模拟方法能够准确地模拟各类水土耦合问题。通过二次开发的CFD-DEM流固耦合模块实现离散元软件PFC3D与计算流体力学软件OpenFOAM之间的力学信息交互,利用颗粒水下自由沉降验证该方法的可行性。利用PFC3D软件模拟室内循环三轴试验标定出具有真实饱和砂土动力特性的数值砂样。根据已有的参数信息以及耦合模拟方法建立了饱和砂土的场地液化模型。模拟结果表明,离散元法能够复现室内砂土液化试验,标定参数可应用于场地液化模拟;单颗粒沉降速度与理论解一致验证了CFD-DEM耦合方法的准确性;峰值加速度0.25g下不同深度处土体均会发生液化,液化时超孔压比无法达到1,超孔压累计值由浅层往深层递增;液化后土体强度自下而上逐渐恢复,再固结的场地土体结构呈现均匀化发展趋势。展开更多
Crustal stresses play an important role in both exploration and development in the oil and gas industry.However,it is difficult to simulate crustal stress distributions accurately,because of the incompatibilities that...Crustal stresses play an important role in both exploration and development in the oil and gas industry.However,it is difficult to simulate crustal stress distributions accurately,because of the incompatibilities that exist among different software.Here,a series of algorithms is developed and integrated in the Petrel2ANSYS to carry out two-way conversions between the 3D attribute models that employ corner-point grids used in Petrel and the 3D finite-element grids used in ANSYS.Furthermore,a modified method of simulating stress characteristics and analyzing stress fields using the finite-element method and multiple finely resolved 3D models is proposed.Compared to the traditional finite-element simulation-based approach,which involves describing the heterogeneous within a rock body or sedimentary facies in detail and simulating the stress distribution,the single grid cell-based approach focuses on a greater degree on combining the rock mechanics described by 3D corner-point grid models with the finely resolved material characteristics of 3D finite-element models.Different models that use structured and unstructured grids are verified in Petrel2ANSYS to assess the feasibility.In addition,with minor modifications,platforms based on the present algorithms can be extended to other models to convert corner-point grids to the finite-element grids constructed by other software.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
A fixed artificial source(greater than 200 kW) was used and the source location was selected at a high resistivity region to ensure high emission efficiency. Some publications used the "earth-ionosphere" mod...A fixed artificial source(greater than 200 kW) was used and the source location was selected at a high resistivity region to ensure high emission efficiency. Some publications used the "earth-ionosphere" mode in modeling the electromagnetic(EM) fields with the offset up to a thousand kilometer, and such EM fields still have a signal/noise ratio of 10-20 dB. This means that a new EM method with fixed source is feasible, but in their calculation, the displacement in air was neglected. In this work, some three-layer modeling results were presented to illustrate the basic EM fields' characteristics in the near, far and waveguide areas under "earth-ionosphere" mode, and a standard is given to distinguish the boundary of near, far and waveguide areas. Due to the influence of the ionosphere and displacement current in the air, the "earth-ionosphere" mode EM fields have an extra waveguide zone, where the fields' behavior is very different from that of the far field zone.展开更多
文摘砂土液化是常见的地震灾害,目前应用于研究砂土液化动力特性的室内试验以及模型试验还不能全面反映土体液化全过程。计算流体动力学(computational fluid dynamics,CFD)与离散元法(discrete element method,DEM)耦合模拟方法能够准确地模拟各类水土耦合问题。通过二次开发的CFD-DEM流固耦合模块实现离散元软件PFC3D与计算流体力学软件OpenFOAM之间的力学信息交互,利用颗粒水下自由沉降验证该方法的可行性。利用PFC3D软件模拟室内循环三轴试验标定出具有真实饱和砂土动力特性的数值砂样。根据已有的参数信息以及耦合模拟方法建立了饱和砂土的场地液化模型。模拟结果表明,离散元法能够复现室内砂土液化试验,标定参数可应用于场地液化模拟;单颗粒沉降速度与理论解一致验证了CFD-DEM耦合方法的准确性;峰值加速度0.25g下不同深度处土体均会发生液化,液化时超孔压比无法达到1,超孔压累计值由浅层往深层递增;液化后土体强度自下而上逐渐恢复,再固结的场地土体结构呈现均匀化发展趋势。
基金Project(2017ZX05013002-002)supported by Major National Science and Technology Projects of ChinaProject(RIPED-2016-JS-276)supported by Petro-China Research Institute of Petroleum Exploration and Development
文摘Crustal stresses play an important role in both exploration and development in the oil and gas industry.However,it is difficult to simulate crustal stress distributions accurately,because of the incompatibilities that exist among different software.Here,a series of algorithms is developed and integrated in the Petrel2ANSYS to carry out two-way conversions between the 3D attribute models that employ corner-point grids used in Petrel and the 3D finite-element grids used in ANSYS.Furthermore,a modified method of simulating stress characteristics and analyzing stress fields using the finite-element method and multiple finely resolved 3D models is proposed.Compared to the traditional finite-element simulation-based approach,which involves describing the heterogeneous within a rock body or sedimentary facies in detail and simulating the stress distribution,the single grid cell-based approach focuses on a greater degree on combining the rock mechanics described by 3D corner-point grid models with the finely resolved material characteristics of 3D finite-element models.Different models that use structured and unstructured grids are verified in Petrel2ANSYS to assess the feasibility.In addition,with minor modifications,platforms based on the present algorithms can be extended to other models to convert corner-point grids to the finite-element grids constructed by other software.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
基金Projects(41204054,41541036,41604111)supported by the National Natural Science Foundation of China
文摘A fixed artificial source(greater than 200 kW) was used and the source location was selected at a high resistivity region to ensure high emission efficiency. Some publications used the "earth-ionosphere" mode in modeling the electromagnetic(EM) fields with the offset up to a thousand kilometer, and such EM fields still have a signal/noise ratio of 10-20 dB. This means that a new EM method with fixed source is feasible, but in their calculation, the displacement in air was neglected. In this work, some three-layer modeling results were presented to illustrate the basic EM fields' characteristics in the near, far and waveguide areas under "earth-ionosphere" mode, and a standard is given to distinguish the boundary of near, far and waveguide areas. Due to the influence of the ionosphere and displacement current in the air, the "earth-ionosphere" mode EM fields have an extra waveguide zone, where the fields' behavior is very different from that of the far field zone.
