共查询到20条相似文献，搜索用时 218 毫秒
1.
Chongbin Zhao B. E. Hobbs A. Ord Shenglin Peng Liangming Liu H. B. Mühlhaus 《Mathematical Geology》2006,38(8):987-1010
Exact analytical solutions have been derived rigorously for the pore-fluid velocity, pore-fluid-flow focusing factor, stream
function and excess pore-fluid pressure around and within a buried inclined elliptic inclusion in pore-fluid-saturated porous
rocks. The geometric characteristics of the buried inclined elliptic inclusion are represented by the aspect ratio and dip
angle of the inclusion, while the hydrodynamic characteristic is represented by the permeability ratio of the elliptic inclusion
to its surrounding rock. Since an elliptic inclusion of any aspect ratio can be used to approximately represent geological
faults and cracks, the present analytical solutions can be used to investigate the pore-fluid-flow patterns around buried
faults and cracks within the crust of the Earth. Therefore, the present analytical solution not only provides a better understanding
of the physics behind the pore-fluid-flow focusing problem around and within buried faults and cracks, but also provides a
valuable benchmark solution for validating any numerical method in dealing with this kind of pore-fluid-flow focusing problem.
The pore-fluid-flow focusing factor of a buried elliptic inclusion is demonstrated to be dependent on the aspect ratio, the
permeability ratio and the dip angle. 相似文献
2.
Chongbin Zhao Bruce E. Hobbs Alison Ord Peter Hornby Hans Mühlhaus Shenglin Peng 《Computers and Geotechnics》2008,35(3):357-371
In this paper, theoretical and numerical methods are used to investigate pore-fluid flow focused temperature distribution patterns around geological faults and cracks of any length-scales in hydrothermal systems. If the far field inflow is uniform and the long axis of an elliptical fault of any length-scale is parallel to the far field inflow direction, a complete set of analytical solutions has been presented for the pore-fluid velocity, stream function and excess pore-fluid pressure around the elliptical fault embedded in fluid-saturated porous media. Because the analytical solutions are explicitly expressed in the conventional Cartesian coordinate system, not only can they be used to gain a theoretical insight into the pore-fluid flow patterns around geological faults and large cracks, but also they can be used as valuable benchmark solutions for validating any numerical methods. After a finite element computational model is validated by comparing the numerical solutions with the present analytical solutions, it is used to investigate pore-fluid flow focused heat transfer around geological faults in hydrothermal systems. Some interesting conclusions in relation to the effects of geological faults on pore-fluid flow focused heat transfer have been made through both the theoretical and the numerical analyses. 相似文献
3.
利用复变函数法和Green函数法给出了SH波对垂直界面附近椭圆形夹杂散射问题的解析解答。首先,将待求的半空间模型沿垂直边界分割为区域Ⅰ和区域Ⅱ两个直角域。通过保角映射的方法将区域Ⅰ内椭圆形夹杂的外域映射为单位圆外域,并利用镜像方法构造出两个区域内满足直角域边界条件的散射波场及适用的Green函数;其次,利用界面契合思想,通过在界面处添加附加力系的方法建立起满足界面处位移和应力连续条件的无穷代数方程组,并截断有限项求解;最后,给出了求解地表位移幅值的具体算例。结果表明,入射波数、入射角度、夹杂位置、垂直界面以及材料参数都对地表位移幅值的分布有影响。 相似文献
4.
Theoretical investigation of convective instability in inclined and fluid-saturated three-dimensional fault zones 总被引：2，自引：0，他引：2
Chongbin Zhao B.E. Hobbs A. Ord Shenglin Peng H.B. Mühlhaus Liangming Liu 《Tectonophysics》2004,387(1-4):47-64
The convective instability of pore-fluid flow in inclined and fluid-saturated three-dimensional fault zones has been theoretically investigated in this paper. Due to the consideration of the inclined three-dimensional fault zone with any values of the inclined angle, it is impossible to use the conventional linear stability analysis method for deriving the critical condition (i.e., the critical Rayleigh number) which can be used to investigate the convective instability of the pore-fluid flow in an inclined three-dimensional fault zone system. To overcome this mathematical difficulty, a combination of the variable separation method and the integration elimination method has been used to derive the characteristic equation, which depends on the Rayleigh number and the inclined angle of the inclined three-dimensional fault zone. Using this characteristic equation, the critical Rayleigh number of the system can be numerically found as a function of the inclined angle of the three-dimensional fault zone. For a vertically oriented three-dimensional fault zone system, the critical Rayleigh number of the system can be explicitly derived from the characteristic equation. Comparison of the resulting critical Rayleigh number of the system with that previously derived in a vertically oriented one has demonstrated that the characteristic equation of the Rayleigh number is correct and useful for investigating the convective instability of pore-fluid flow in the inclined three-dimensional fault zone system. The related numerical results from this investigation have indicated that: (1) the convective pore-fluid flow may take place in the inclined three-dimensional fault zone; (2) if the height of the fault zone is used as the characteristic length of the system, a decrease in the inclined angle of the inclined fault zone stabilizes the three-dimensional fundamental convective flow in the inclined three-dimensional fault zone system; (3) if the thickness of the stratum is used as the characteristic length of the system, a decrease in the inclined angle of the inclined fault zone destabilizes the three-dimensional fundamental convective flow in the inclined three-dimensional fault zone system; and that (4) the shape of the inclined three-dimensional fault zone may affect the convective instability of pore-fluid flow in the system. 相似文献
5.
The pressure path of solid inclusions in minerals: the retention of coesite inclusions during uplift 总被引：4，自引：0，他引：4
An inclusion model allows calculation of the stresses in and around minerals included in other minerals: during changes in pressure and temperature. The equations are applied to illustrate cooling and uplift histories of quarzo-feldspathic rocks from 500°C and various pressures to ambient conditions. Even in the absence of pore-fluid pressure, microfractures may open at external pressures of 200–400 MPa and temperatures of 200–400°C due to differential volume changes of the constituent minerals. Coesite included in garnet cannot have formed during progressive metamorphism from quartz at lithostatic pressures below the coesite stability field because of differential volume changes. Coesite inclusions are captured by their host minerals at ultra-high pressures and they persist to lower pressures because of the large volume increase occurring at the coesite to -quartz transition. The Pi-T path followed by the SiO2 inclusion traces the boundary between the stability fields of coesite and -quartz until radial fractures develop in the host at low external pressure. 相似文献
6.
We extend Riemannian wavefield extrapolation to pre- and post-stack migration using 2D logpolar coordinate system in order to enhance imaging of subsurface structures with imaging of dipping reflectors and turning waves. The logpolar coordinate is a coordinate system adequate for propagating both the source and receiver wavefields. The logpolar extrapolation wavenumber introduces an isotropic slowness model stretch to the single square root operator that enables the use of Cartesian finite-difference extrapolators for propagating wavefields in logpolar meshes. Wavefield extrapolation in this coordinate does not require any ray tracing at all and it can be readily extended to high-order finite-difference schemes since it is an analytic extrapolation operator. Two migration examples illustrate the advantages of the logpolar coordinate for imaging complex geological structures. We compared results of post- and pre-stack depth migration in logpolar and elliptic coordinate systems and we found that steeply dipping events can be better imaged in logpolar coordinate than in elliptic coordinate. 相似文献
7.
在反射波地震勘探中使用动校正确定反射点位置用到很多假设,如勘探深度要远远大于炮检距;假定倾斜反射界面的反射点与水平反射界面反射点都位于炮检距的中点;近似认为倾斜反射界面的动校正量等于水平反射界面的动校正量;假定反射面倾角固定且较小等,上述假设一定会造成勘探误差。由于反射点位置和反射面倾角未知,所以理论上无法唯一地确定反射点位置。如果反射波传播的介质的波速一定,从炮点发出的地震波,经反射点后,在接收点被接收,其可能的反射点是椭圆的一部分,但还不能唯一确定反射点;再取炮点和另外一个接收点,其可能的反射点是另外一个椭圆的一部分。如果假定反射面是平面,可以是水平面,也可以是有固定倾角的倾斜平面,该平面在地震波射线平面内是一条直线,该直线一定是两椭圆的公切线。把两椭圆方程和切线方程联立,就可以求解出公切点位置,公切点位置就是反射点位置,这就是用双椭圆确定反射点位置的方法。通过建立模型对勘探方法进行了检验,证实了用双椭圆方法确定反射点位置的有效性。双椭圆方法有一个重要的副产品,就是在确定反射面位置的同时计算出反射面的视倾角。 相似文献
8.
在反射波地震勘探中使用动校正确定反射点位置用到很多假设,如勘探深度要远远大于炮检距;假定倾斜反射界面的反射点与水平反射界面反射点都位于炮检距的中点;近似认为倾斜反射界面的动校正量等于水平反射界面的动校正量;假定反射面倾角固定且较小等,上述假设一定会造成勘探误差。由于反射点位置和反射面倾角未知,所以理论上无法唯一地确定反射点位置。如果反射波传播的介质的波速一定,从炮点发出的地震波,经反射点后,在接收点被接收,其可能的反射点是椭圆的一部分,但还不能唯一确定反射点;再取炮点和另外一个接收点,其可能的反射点是另外一个椭圆的一部分。如果假定反射面是平面,可以是水平面,也可以是有固定倾角的倾斜平面,该平面在地震波射线平面内是一条直线,该直线一定是两椭圆的公切线。把两椭圆方程和切线方程联立,就可以求解出公切点位置,公切点位置就是反射点位置,这就是用双椭圆确定反射点位置的方法。通过建立模型对勘探方法进行了检验,证实了用双椭圆方法确定反射点位置的有效性。双椭圆方法有一个重要的副产品,就是在确定反射面位置的同时计算出反射面的视倾角。 相似文献
9.
10.
层状可压缩岩基三维固结问题的状态空间解 总被引：1，自引：0，他引：1
采用状态空间法求解层状可压缩岩基的三维固结问题。首先从直角坐标系下考虑可压缩性的三维Biot固结问题的控制方程出发,通过Laplace-Fourier变换得到状态空间方程,解此方程并通过Cayley-Hamilton定理,得到单层可压缩岩基三维固结问题的传递矩阵;然后利用传递矩阵法,结合层间连续性条件和边界条件,得到了层状可压缩岩基三维固结问题在积分变换域内的解答;最后应用Laplace-Fourier逆变换技术,得到层状可压缩岩基三维固结问题在物理域内的理论解答。编制了相应的计算程序,进行了数值计算与分析,证明了压缩性对岩基固结问题的影响 相似文献
11.
ABSTRACT The decrepitation behaviour of fluid inclusions in quartz at one atmosphere confining pressure has been evaluated using pure H2O synthetic inclusions formed by healing fractures in natural quartz. Three different modes of non-elastic deformation, referred to as stretching, leakage or partial decrepitation, and total decrepitation have been observed. The internal pressure required to initiate non-elastic deformation is inversely related to inclusion size according to the equation: internal pressure (kbar) = 4.26 D-0.423 where D is the inclusion diameter in microns. Regularly shaped inclusions require a higher internal pressure to initiate non-elastic deformation than do irregularly shaped inclusions of similar size. Heating inclusions through the α/β quartz inversion results in mechanical instability in the quartz crystal and leads to mass decrepitation of inclusions owing to structural mismatches generated by pressure gradients in the quartz around each inclusion. Long-term heating experiments (∼2 years) suggest that the internal pressure required to initiate non-elastic deformation does not decrease significantly with time and indicates that short-lived thermal fluctuations in natural systems should not alter the inclusion density and homogenization temperature. Inclusions that do exhibit decreased density (higher homogenization temperature) are, however, always accompanied by a change in shape from irregular to that of a negative crystal. Observations of this study are consistent with elasticity theory related to fracture generation and propagation around inclusions in minerals. These results indicate that an inclusion will not be influenced by a neighbouring inclusion, or other defect in the host phase, as long as the distance between the two is >2–4 diameters of the larger of the two inclusions. 相似文献
12.
Residual pressure around mineral inclusions in diamond can provide useful information on the depth of diamond origin. Differential stress between an inclusion and host diamond arises from differences in thermal expansion and compressibility between host diamond and minerals. We determined residual pressure around mineral inclusions in a diamond from the Internationalnaya pipe, Yakutia, Russia, using the three-dimensional Raman mapping system developed recently by our group. The maximum residual pressures around the olivine and chromite inclusions were determined to be 0.69 GPa and 0.75 GPa, respectively. We proposed an advanced method for determining simultaneously pressure and temperature conditions where the mineral inclusions were trapped in the host diamond. The obtained values were 3.0 GPa and 447 °C, but these values are lower than typical P-T conditions in the mantle. Several technical possibilities for the discrepancy are discussed. 相似文献
13.
Fluid overpressure at the base of low-permeability strata reduces effective stress, allowing for gravitational sliding of the overlying cover. The force driving sliding is the slope-parallel component of the weight of the cover, whereas the resisting forces are the friction at the base of the cover and the buttressing resistance to shortening, which can be critically reduced by incision at the base of the slope. We developed an analytical model and undertook a series of analogue experiments to better understand the evolution of a sedimentary cover sliding above a low-permeability layer subjected to fluid overpressure. Where a downslope buttress was present, the sliding sheet length decreased with increasing pore-fluid pressure. In the absence of such buttress, the slide's length increased exponentially with increasing pore-fluid pressure. Another important difference dealt with geometry and kinematics. Buttressed slides consisted of one large slope-parallel mass rigidly translated and bounded by downslope thrusts and upslope normal faults. With increasing pore-fluid pressure, the contractional structures propagated upslope. By contrast, non-buttressed slides showed intense strain: deformation started with normal faults forming near the incision, then propagating upslope throughout the slide's evolution. 相似文献
14.
激光拉曼光谱法在单个流体包裹体研究中的应用进展 总被引：2，自引：2，他引：0
激光拉曼光谱法以其非破坏性、高灵敏度和高分辨率等特性,一直以来都是研究流体包裹体的重要方法之一。近年来,围绕激光拉曼光谱在流体包裹体中的应用而展开的研究主要集中在半定量-定量测试方面,即在准确定性的基础上,采用高斯-洛仑兹去卷积分峰、低温原位等定量方法获取流体包裹体的成分和含量,从而克服了激光拉曼光谱法应用于溶液中阳离子的定量分析的灵敏度、准确度较低的问题,不仅能获得流体包裹体中一些常温下不具拉曼活性的盐类物质拉曼特征峰信息,还能根据特征峰与浓度、内压之间的线性关系,进一步对盐度和压力等性质进行测定,从而拓展了激光拉曼光谱法在流体包裹体中的应用范围。本文对激光拉曼光谱法应用于分析流体包裹体成分、盐度、同位素和压力所取得的最新进展进行了评述,并认为随着分析测试技术的不断进步,激光拉曼光谱法未来的分析方向也将继续围绕多元复杂体系及其定量方面的研究展开。 相似文献
15.
This paper deals with analytical and numerical modelling of the internal stress generated in argillaceous rocks during humidification/desiccation processes, which is an essential issue for damage study. This local stress field arises from two mechanisms: (i) complex interactions between free swelling/shrinking clay matrix and non‐strained inclusions of carbonate and quartz and (ii) a self‐restraint effect induced by the moisture gradient during the transient moisture exchange process. The inclusion–matrix interaction is investigated in different cases. Firstly, the analytical solution of the stress around a cylindrical inclusion embedded in an infinite swelling matrix is derived: The inclusion would suffer tension (compression) under humidification (desiccation), and the resulting cracking patterns are discussed. Then, the problem of two inclusions with different distances in an infinite swelling matrix is considered, and it is shown that the local stress around an inclusion will be perturbed and amplified by neighbouring inclusions. Finally, an inclusion outcropping at the free surface of a swelling matrix is modelled as to investigate the effect of free surface: The inclusion–matrix interface undergoes shear stresses of which the maximum is found at the free surface. In addition to the inclusion–matrix interaction, the self‐restraint effect is investigated: The induced stress is maximal at the beginning of humidification/desiccation processes and vanishes gradually with time. The quantity of the self‐restraint stress is strongly controlled by the hydric loading rate. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
16.
了解单个石油包裹体详细的组分信息不仅有助于微观上定量认识油气充注过程中其成熟度的演化,还对精细刻画油气成藏过程中动力演化具有重要意义.根据石油包裹体热动力学原理,利用10个成熟度依次增大的已知组分的原油在匹配其饱和压力和改进气、液相摩尔体积计算前提下,建立了捕获不同组分原油的包裹体均一温度(Thoil)与其在室温下(20 ℃)气泡充填度(Fv)关系的标准图版,从而通过对比实测的Thoil和其室温下Fv与标准Thoil-Fv图版,反推出油包裹体近似等效流体组分.对东营凹陷北带民丰深层油包裹体进行系统荧光分析、显微测温和包裹体体积3D重构等一系列的测定.结果表明丰8井和丰深1井各发生过两幕油气充注,其中丰8井第一幕油气充注(甲烷摩尔含量位于31%~35%之间)对其油气成藏贡献最大;而丰深1井第二幕充注的油气(甲烷摩尔含量大于60%)控制丰深1井凝析油藏成藏.从而为进一步厘定北带深层油气充注PVTxt史奠定基础. 相似文献
17.
In present study, the non-linear variations of soil compressibility, hydraulic and electro-osmosis conductivities were analyzed through laboratory experiments, and incorporated in a one-dimensional model. The analytical solutions for excess pore water pressure and degree of consolidation were derived, and numerical simulations were performed to verify its effectiveness. The results indicated that both the non-linear variations of hydraulic and electro-osmosis conductivities showed remarkable impacts on the excess pore water pressure and degree of consolidation, especially for soils with relative high compressibility. A further comparison with previous analytical solutions indicated that more accurate predictions could be obtained with the proposed analytical solutions. 相似文献
18.
Petroleum inclusions in sedimentary basins: systematics, analytical methods and applications 总被引：21，自引：0，他引：21
Fluid dynamics in sedimentary basins is of tremendous interest, both from a scientific and an economic point of view. Integration between fluid inclusion and present-day fluid data provides the time aspect necessary for reconstruction of fluid flow paths, and can be used for mapping fluid dynamics both on a regional basin scale or on the more local scale of petroleum reservoirs. This paper presents a review of analytical and modelling methods for petroleum in fluid inclusions. Essentially, four types of data for petroleum inclusions can be established through analysis and modelling: (1) textures (which give indirect information of the time aspect), (2) fluid composition, (3) fluid properties and (4) pressure–temperature of trapping. During the last decade, development of analytical methods for determination of inclusion fluid compositions has taken place. Traditional correlation studies for characterisation of maturity and source facies may now therefore include fluid inclusion data. The development has also been directed towards an improved understanding of the physical properties of petroleum inclusions. Although these methods are in their early stages concerning precision and accuracy, data necessary for fluid flow modelling, such as fluid densities, viscosities, pressure and temperature, can be estimated. 相似文献
19.
《Journal of Structural Geology》2003,25(10):1713-1733
This paper aims at reconstructing paleostress history and deciphering the pore-fluid pressure conditions during the reactivation of a regional-scale fault. Paleostress analyses of the mesoscopic structures suggest that three successive events of reactivation on the Great Boundary Fault occurred in thrust-type, strike-slip type and thrust-type tectonic-settings, respectively. Whereas the pore-fluid pressure was supralithostatic during the first and third events of reactivation, it was sublithostatic during the second event. Each event of reactivation induced a fluid pressure gradient, which resulted in the focused and enhanced flow of syntectonic fluids within the high permeability locales, primarily fractures and faults.Fluid inclusion study on strike-slip veins reveals that the syntectonic fluids were highly dense, Na–Ca–Cl brines of formational water origin. Stratigraphic evidence in favour of a 2-km-thick column of overburden above Kaimur sandstone beds implies that the strike-slip faulting occurred at 160–202 °C temperature and 53 MPa pressure. Variation in homogenization temperature reflects fluctuation in pore-fluid pressure during entrapment of syntectonic fluids and points to seismic pumping as a possible mechanism of fluid flow during faulting. High paleogeothermal gradient, obtained by fluid inclusion data, is ascribed to the high heat flow due to crustal stretching during the Proterozoic rifting, the basal and intermittent volcanism in the basin, and occurrence of Berach granite as the basement. 相似文献
20.
Shear‐lag analysis is used to obtain closed‐form solutions for the problem of a stiff inclusion embedded in a poroelastic soil matrix. The following assumptions are made: the soil matrix and the inclusion are elastic; plane strain conditions apply; and shear stresses at the soil‐inclusion interface follow Coulomb's friction law. Two solutions are obtained, the first one for drained conditions where no excess pore pressures are generated, and the second one for undrained conditions where excess pore pressures are produced and the soil does not change volume during pullout. The solutions are verified by comparing analytical predictions with numerical results obtained using a finite element method. Predictions from the analytical solutions are also compared with results from experiments conducted in a large‐scale pullout box. Both comparisons show good agreement. The analytical solution shows that the pullout capacity in drained and undrained conditions is overall independent of the relative stiffness of the soil and the inclusion. The most important factor controlling the pullout capacity is the coefficient of friction between the soil and the inclusion. Both drained and undrained pullout capacities increase with the coefficient of friction; although the drained capacity shows a proportional increase, it is not so for the undrained capacity. The ratio of undrained to drained pullout capacity is about 0.9 for friction coefficients smaller than 0.2, but can be as small as 0.6 for a coefficient of friction of 1.0. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献