A new scheme for numerical modelling of flow and transport processes in 3D fractured porous media

The objective of this work is to develop a new numerical approach for the three-dimensional modelling of flow and transient solute transport in fractured porous media which would provide an accurate and efficient treatment of 3D complex geometries and inhomogeneities. For this reason, and in order to eliminate as much as possible the number of degrees of freedom, the fracture network, fractures and their intersections, are solved with a coupled 2D–1D model while the porous matrix is solved independently with a 3D model. The interaction between both models is accounted for by a coupling iterative technique. In this way it is possible to improve efficiency and reduce CPU usage by avoiding 3D mesh refinements of the fractures. The approach is based on the discrete-fracture model in which the exact geometry and location of each fracture in the network must be provided as an input. The formulation is based on a multidimensional coupling of the boundary element method-multidomain (BEM-MD) scheme for the flow and boundary element dual reciprocity method-multidomain (BE-DRM-MD) scheme for the transport. Accurate results and high efficiency have been obtained and are reported in this paper.     

Convective stability analysis of the long-term storage of carbon dioxide in deep saline aquifers

Deep saline aquifers are one of the most suitable geologic formations for carbon sequestration. The linear and global stability analysis of the time-dependent density-driven convection in deep saline aquifers is presented for long-term storage of carbon dioxide (CO₂). The convective mixing that can greatly accelerate the CO₂ dissolution into saline aquifers arises because the density of brine increases upon the dissolution of CO₂ and such a density difference may induce instability. The effects of anisotropic permeability on the stability criteria, such as the critical time for the appearance of convective phenomena and the critical wavelength of the most unstable perturbation, are investigated with linear and global stability analysis. The linear stability analysis provides a sufficient condition for instability while the global stability analysis yields a sufficient condition for stability. The results obtained from these two approaches are not exactly the same but show a consistent trend, both indicating that the anisotropic system becomes more unstable when either the vertical or horizontal permeability increases.     

Mantle process beneath Philippine Sea back-arc spreading ridges: A synthesis of peridotite petrology and tectonics

Abstract In order to obtain a general view of the mantle process beneath a back‐arc basin spreading ridge, the diversity of peridotite petrology and tectonic occurrences in two back‐arc basin spreading ridges from the Philippine Sea were examined: the Parece Vela Rift and the Mariana Trough. The Parece Vela Basin spreading ridge (Parece Vela Rift) was a physically fast/intermediate‐spreading ridge, although many tectono‐magmatic features resemble those of slow‐ to ultraslow‐spreading ridges. Two unusual features of the Parece Vela Rift further demonstrate the uniqueness of the ridge: full‐axial development of oceanic core complexes and exposure of mantle peridotite at segment midpoints. The Parece Vela Rift yields a lithological assemblage of residual but still fertile lherzolite/harzburgite, plagioclase‐bearing harzburgite and dunite; similar assemblages are reported from the equatorial Mid‐Atlantic Ridge at the Romanche Fracture Zone and the ultraslow‐spreading ridges from the Indian and Arctic Oceans. The tectono‐magmatic characteristics of the Parece Vela Rift suggest that diffuse porous melt flow and pervasive melt–mantle interaction were the important mantle processes there. Globally, this ‘porous melt flow‐type’ mantle process is likely to occur beneath a segment midpoint of the ridge having a thick lithosphere, typically an ultraslow‐spreading ridge. In contrast, the Mariana Trough is a typical slow‐spreading ridge, exposing mantle peridotite at segment ends. The Mariana Trough yields a lithological assemblage of residual harzburgite and veined harzburgite, a common assemblage among the global abyssal peridotite suite. The tectono‐magmatic characteristics of the Mariana Trough suggest that channeled melt/fluid flow and limited melt–mantle interaction are the important mantle processes there, because of the colder wall‐rock peridotite in the segment end. This ‘channeled melt flow‐type’ mantle process is likely to occur in the shallow lithospheric mantle at the segment ends of any spreading ridges.     

非均匀介质中交错网格高阶有限差分数值模拟

地震波场的数值模拟一直是地球物理学的一个重要的研究领域,而在数值正演模拟方法的研究中,计算精度和计算效率是评价该方法有效性及优越性的二个关键问题。这里从一阶速度—应力弹性波动方程出发,着重介绍如何构造离散化模型的网格,如何求解空间导数,如何选取边界条件等内容,从而更有效地提高数值计算的精度与计算效率。文中构造了不同类型的介质模型,并在交错网格中,利用高阶有限差分模拟非均匀介质的波场传播。模拟结果表明,该方法实现简单,具有很好地稳定性和较高的精度,能够直观、高效地反映出介质中波场的传播规律。     

深层地热水开采与地面沉降的关系研究

本文根据天津塘沽区深层地沉观测标观测到的地面沉降数据,分析了深层地下热水开采对地面沉降的影响比例.在分析地热水开采引起地面沉降机理的基础上,通过统计回归方法建立了深层孔隙型地热水开采量与地面沉降量的相关模型;提出了进行地热尾水回灌和建立深层沉降监测网以加强综合研究的防治对策.     

多孔介质中水敏性应用的研究综述

水敏现象在石油、岩土和环境工程等领域广泛存在。有效地控制水敏现象的发生，不仅具有科学意义，还有一定的经济价值。本文在广泛查阅国内外研究成果的基础上，总结了水敏性在石油、岩土和环境工程方面应用的研究进展，并在此基础上提出了水敏性研究的发展方向。     

Errors due to the anisotropic-common-ray approximation of the coupling ray theory

The common-ray approximation eliminates problems with ray tracing through S-wave singularities and also considerably simplifies the numerical algorithm of the coupling ray theory for S waves, but may introduce errors in travel times due to the perturbation from the common reference ray. These travel-time errors can deteriorate the coupling-ray-theory solution at high frequencies. It is thus of principal importance for numerical applications to estimate the errors due to the common-ray approximation applied. The anisotropic-common-ray approximation of the coupling ray theory is more accurate than the isotropic-common-ray approximation. We derive the equations for estimating the travel-time errors due to the anisotropic-common-ray (and also isotropic-common-ray) approximation of the coupling ray theory. The errors of the common-ray approximations are calculated along the anisotropic common rays in smooth velocity models without interfaces. The derivation is based on the general equations for the second-order perturbations of travel time.     

多孔介质中非均匀流动特性的染色示踪试验研究

通过土壤染色剂进行的4组试验,对不同介质结构条件下的水流和溶质非均匀运动规律,非均匀流动变异信息分布特征关系以及全局性非均匀流动示踪方法进行了研究。结果表明,即使在相对比较均匀的介质条件下,流动也表现出明显的非均匀特性;对数正态分布能够较好的反映水流的运动分布模式,相比水流运动,溶质的运动和分布规律明显不同,表现出更多的不确定性和变异性。     

线性粘弹介质中地震波场数值模拟

本文将以往两种粘弹介质中地震波模拟方法的优点结合起来，以模型理论和积分本构方程为基础，从理论上分析了模型对地震波场的影响；采用交错网格有限差分法对粘弹介质中的地震波进行数值模拟.数值计算结果表明该方法不仅便于计算，同时也便于从力学的角度来分析地震波的传播.数值计算结果与理论分析一致，说明这种方法可以更为有效地模拟粘弹介质中地震波的传播.     

TTI介质的相速度研究

推导了二维TTI介质的相速度表达式,并且依据推导出来的相速度表达式,模拟并分析了二维TTI介质相速度的传播快照以及TI介质相速度的传播快照;对比并分析了TTI介质和TI介质模型的相速度理论计算值的X分量特征的差异。TTI介质的相速度研究具有较高的理论研究价值和实际应用价值．