镜质体成烃反应动力学模型的标定及其在热史恢复中的应用

镜质体反射率（Ｒｏ）与其生烃率有关，因此，由Ｒｏ应可求取其生烃率；同时，在求取了镜质体成烃反应动力学参数之后，若已知其受热史，则可计算出镜质体的生烃率；反过来，若已知镜质体的生烃率，应可反推其受热史。基于这一基本原理，本文在建立和标定镜质体成烃反应动力学模型的基础上推导了可将其应用于存在着抬升、剥蚀和古地温非线性变化等复杂地质条件下热史恢复的原理模型，结合海拉尔盆地海参７井的实测Ｒｏ值，有效地恢复了其古热史     

镜质体成烃反应动力学模型的标定及其在热史恢复中的应用

镜质体反射率（Ｒｏ）与其生烃率有关，因此，由Ｒｏ应可求取其生烃率；同时，在求取了镜质体成烃反应动力学参数之后，若已知其受热史，则可计算出镜质体的生烃率；反过来，若已知镜质体的生烃率，应可反推其受热史。基于这一基本原理，本文在建立和标定镜质体成烃反应动力学模型的基础上推导了可将其应用于存在着抬升、剥蚀和古地温非线性变化等复杂地质条件下热史恢复的原理模型，结合海拉尔盆地海参７井的实测Ｒｏ值，有效地恢复了其古热史。     

软岩巷道弹塑性变形的理论分析

本文采用莫尔－库仑屈服准则和非关联塑性流动法则分析了轴对称软岩巷道的变形规律。从轴对称平面应变问题的基本方程出发，导出了软岩巷道弹塑性变形的位移理论解答，并给出了现有应力与位移理论解的适用条件。通过实例计算结果表明，本文得出的位移理论解与有限元数值解相当吻合。     

Three-dimensional heat,moisture and air transfer in unsaturated soils

A new three-dimensional numerical model of coupled heat, moisture and air transfer in unsaturated soil is presented. In particular, the model accommodates moisture transfer in the form of liquid and vapour flow and heat transfer arising from conduction, convection and latent heat of vaporization. The bulk flow of dry air and the movement of air in a dissolved state are also included. The theoretical basis of the model, the finite element solution of the spatial terms and finite difference solution of the temporal terms are briefly presented. Attention is focused on the verification of the new numerical solution. This is achieved via comparisons with independent solutions of heat, moisture and air transfer in an unsaturated soil. The physical problem considered includes the highly non-linear hydraulic properties of sand. Thermal conductivity is also included as a function of soil moisture content. Excellent correlation of results is shown thus providing confidence in the new model. The new model is also applied to a number of test cases which illustrate the need for the development of a model which can fully include three-dimensional behaviour. In particular, three applications are presented each increasing in complexity. The first application illustrates three-dimensional heat transfer. This particular application is verified against existing commercial finite element software. Subsequent applications serve to illustrate how the coupled processes of heat moisture and air transfer combine to yield three-dimensional problems even within a simple geometric domain. Visualization of three-dimensional results is also addressed. © 1998 by John Wiley & Sons, Ltd.     

土体一维冻结问题温度场半解析解

针对土体的一维冻结过程,采用考虑土体冻结状态下未冻水存在的等效热容模型,建立了其温度场计算的半解析方法, 并与相关文献中的数值解进行了对比,验证了方法的正确性。对等效热容模型与显热容模型进行了对比, 计算结果表明：在瞬态阶段,显热容模型计算的冻结锋面推进速度较等效热容模型慢;而在接近稳态阶段,显热容模型计算的冻土区厚度较等效热容模型厚。在冻结诱导水流较弱的情形下,土体冻结状态下的持水特性越好,采用显热容模型进行温度场计算产生的误差越显著,并且其计算结果应用于人工冻结设计是偏于不安全的     

底部抽真空预压法砂井地基固结解析解

基于谢康和等应变条件砂井地基固结理论和Hansbo砂井固结理论,考虑底部抽真空预压法加固方法中真空作用面位于固结土层底部的实际边界条件,推导出忽略竖向渗流情况下的底部真空预压加固地基固结方程解析解答。根据超孔压固结方程形成过程以及其解析解表达式,分析其与一般负压径向固结解答的区别。通过室内模型试验实测数据与解析模型计算结果的对比表明,不同位置处孔压和固结度计算值与实测结果吻合较好,从而验证了该模型的合理性,同时运用该模型也可有效验证已有关于底部抽真空室内模拟及现场原位试验结果。底部抽真空轴对称固结解析解可为底部抽真空技术的实际工程应用提供基础性的理论支持,推动底部抽真空技术的大规模推广应用。     

Control Function Measures for Hydrodynamic Problems

This short paper shows that different choices made for control functions of information in hydrodynamic flow problems can have significant implications for interpretations of the system. Using as a simple illustration the case of a steady-state, one-dimensional flow model with no internal sources or sinks and with the hydraulic conductivity depending on a single parameter and the distance from the origin, it is shown that, even when a continuous, error-free head data set is provided, statements about the uniqueness or not of the inverse solution are conditioned on the choice of the control function. Care has to be exercised in obtaining physically meaningful results and, depending on the model assumptions and the data available, there may not be acceptable models. It is also shown that there may be more than one model behavior that is acceptable. The results have implications for the hydrodynamic upscaling problem for flow in permeable media, for ensemble averaging methods, and for parameter determination for deterministic models of permeable flow.     

四川盆地中生代地热特征与前陆盆地构造、沉积作用

沉积盆地热体制与油气资源勘探有着密切的联系.基于石油钻井的Ro数据,采用古地温梯度法恢复了四川盆地钻井的最高古地温.结果显示,四川盆地西部、川东南和川东北的钻井,顶部的剥蚀量在3 700 m左右,盆地中部的剥蚀量约3 100 m.钻井古地温梯度在13.2～24.5℃/km之间,古热流在31.2～60.1 mW/m2之间...     

早前寒武纪古地热状态的构造环境意义

区域变质作用是在造山运动中地热流和动力联合作用的体现。在早前寒武纪的造山带中,随着地壳的缩短,古地热流向构造边界聚敛,形成由高温到低温的区域变质相带序列,它既指示古地热流运移的方向,又是构造边界的指向。古地热流值的高低同构造环境的性质有关,如陆内裂陷带的区域变质仅达绿片岩相,而陆缘活动带的则可高达角闪岩相。古地热流从活动带的边界向中心运移,具有向心性,随着构造层的深埋,构成地热槽,它的片段则表现为地热斜。识别这些古地热状态,可为构造环境属性的鉴定提供有力的依据     

Effects of transient uplift/erosion on the surface heat flow and heat generation relationship in presence of small scale asthenospheric convection

In this paper the effect of transient uplift/erosion on the relationship between surface heat flow and heat generation for truncated exponential model of radiogenic heat source distribution and basal asthenospheric convection is investigated. Asthenospheric convection is described by a parameterized model, in the form of a nonlinear heat flux boundary condition involving basal temperature and mantle internal temperature. This boundary condition has been linearized and the analytical solution of the problem is obtained by the eigenvalue-eigenfunction expansion method. The analytical solution is used to derive the nature of surface heat flow and heat generation relationship. The results show that the linear relationship is maintained during the uplift/erosion and the estimates of the slope of the linear relationship are different from the depth scale of the exponential model and increase with the rate of uplift/erosion. The estimates of the reduced heat flow also increase with the rate of uplift/erosion. These results would find applications in the interpretation of linear surface heat flow and heat generation relationship which is observed in different tectonic environment.     

悬臂式刚性墙动力响应解析研究

从水平简谐振动作用下二维土?悬臂式刚性墙计算模型出发,基于波动力学理论,同时考虑土层的竖向应力和竖向位移,对二维场地中悬臂式刚性墙的动力响应特性进行了解析研究。首先对土层振动方程进行变换,得到关于体积应变? 的方程,通过分离变量法进行求解,再回代振动方程得到关于位移运动方程的非齐次方程,结合墙与土层的相互作用条件及远场边界条件得到振动方程定解,进而得到地下刚性墙墙上土压力、墙底剪力及弯矩的更为严格的解析解。将所得解与忽略竖向应力解、忽略竖向位移解进行了对比。研究表明,所得解能多反映出一个共振频率,且当土体泊松比大于0.45时,忽略竖向位移解失去意义。通过参数分析,表明激励频率与土体阻尼因子对墙体动力响应影响较大,考虑的振动模态阶数对墙体动力响应影响较小。     

A mixed solution for two-dimensional unsteady flow in fractured porous media

A mixed finite element–boundary element solution for the analysis of two-dimensional flow in porous media composed of rock blocks and discrete fractures is described. The rock blocks are modelled implicitly by using boundary elements whereas finite elements are adopted to model the discrete fractures. The computational procedure has been implemented in a hybrid code which has been validated first by comparing the numerical results with the closed-form solution for flow in a porous aquifer intercepted by a vertical fracture only. Then, a more complex problem has been solved where a pervious, homogeneous and isotropic matrix containing a net of fractures is considered. The results obtained are shown to describe satisfactorily the main features of the flow problem under study. © 1997 by John Wiley & Sons, Ltd.     

Simulation of chloride migration rates in Paleo Pennar Delta Region, Coastal Andhra Pradesh, India

 Quaternary alluvial aquifers in the paleo Pennar delta region of Andhra Pradesh (Long. 80°0′ and 80°12′; Lat. 14°40′ and 14°20′) constitute an important coastal strip, with potential fresh groundwater resources and several well fields in operation. The lineament patterns and traces of paleo channels provide basic information on the configuration and boundaries of paleo delta and the possible neotectonic movements in the region. The paleo delta region is essentially characterized by freshwater aquifer systems at the near-surface depths, (up to 50 m from ground level) with transmissivity values in the system ranging from 1200 to 2500 m²/day. This groundwater system has been subjected to heavy draft situation over the last two decades, parlicularly near the outer rim of the delta where an unlined brackish-water canal runs parallel to the coast imparting chloride contamination to the adjoining fresh aquifers. A two-dimensional solute transport model solution was applied to assess chloride migration rates inland under different hydraulic stresses, combining finite difference solution of flow equation and the method of characteristic solution of solute transport equation. Groundwater flow and chloride migration patterns/rates were obtained for different simulated stress events in the delta system, and measures required to protect the freshwater resources ara outlined. Received: 2 January 1997 · Accepted: 4 November 1997     

A novel analytical solution for constant‐head test in a patchy aquifer

A mathematical model describing the hydraulic head distribution for a constant‐head test performed in a well situated at the centre of a patchy aquifer is presented. The analytical solution for the mathematical model is derived by the Laplace transforms and the Bromwich integral method. The solution for the hydraulic head has been shown to satisfy the governing equations, related boundary conditions, and continuity requirements for the hydraulic head and flow rate at the interface of the patch and outer regions. An efficient numerical approach is proposed to evaluate the solution, which has an integral covering an integration range from zero to infinity and an integrand consisting the product and square of the Bessel functions. This solution can be used to produce the curves of dimensionless hydraulic head against dimensionless time for investigating the effect of the contrast of formation properties on the dimensionless hydraulic head distribution. Define the ratio of outer‐region transmissivity to patch‐region transmissivity as α. The dimensionless hydraulic head for α=0.1 case is about 2.72 times to that for α=10 case at dimensionless large time (e.g. τ?10⁶) when the dimensionless distance (ρ) equals 10. The results indicate that the hydraulic head distribution highly depends on the hydraulic properties of two‐zone formations. Copyright © 2006 John Wiley & Sons, Ltd.     

Thermoelastic stress due to a rectangular heat source in a semi-infinite medium. Presentation of an analytical solution

The thermoelastic response due to a time-dependent rectangular heat source in a semi-infinite medium is analyzed. The problem originates from studies of nuclear waste repositories in rock. Canisters containing heat-emitting nuclear waste are deposited over a large rectangular area deep below the ground surface. The solution for a time-dependent heat source is obtained from the corresponding instantaneous heat source by superposition. The thermoelastic problem for the instantaneous rectangular heat source in a infinite surrounding is solved exactly. An important step is the introduction of so-called quadrantal heat sources. The solution for the rectangle is obtained from four quadrantal solutions. The solution for the quadrantal heat source depends on the three dimelasionless coordinates only. Time occurs in the scale factors only. The condition of zero normal and shear stresses at the ground surface is fulfilled by using a mirror heat source and a boundary solution. The boundary solution accounts for the residual normal stress at the ground surface. Using a Hertzian potential, a surprisingly simple solution is obtained. The final analytical solution is quite tractable considering the complexity of the initial problem. The solution may be used to test numerical models for coupled thermoelastic processes. It may also be used in more detailed numerical simulations of the process near the heat sources as boundary conditions to account for the three-dimensional global process.     

敏感性分析在烃源岩成熟度史模拟中的应用——以川东北地区普光5井古生界海相烃源岩为例

通过分析输人模型的参数对输出结果的影响,可以确定影响烃源岩成熟度史模拟的敏感性参数.本文应用Easy% Ro化学动力学模型,以普光5井为例,对川东北地区各期构造运动剥蚀厚度、古地表温度和古地温梯度进行了相关的敏感性分析.分析结果表明:研究区下寒武统、下志留统、下二叠统和上二叠统烃源岩现今成熟度状态完全受控于燕山运动晚幕...     

Theoretical study of hydraulically fractured penny-shaped cracks in hot,dry rocks

A simple and analytical method of solving the problem of heat extraction from a penny-shaped crack having both inlet and outlet holes is developed by using a two-dimensional flow model when fluid is injected at a constant flow rate. Temperature distributions in both rock and fluid are obtained by taking into account the hydraulic and thermal growth of the crack. The outlet fluid temperature and energy extraction rate versus time are shown by some illustrative examples. Although the effect of thermal, contraction of rock upon the crack width is large, its influence upon flow rate and crack tip stress intensity factor is sufficiently small within the first several years so that the crack radius can be determined mainly by the mechanical deformation of, the rock.     

A Simple Method for Analysis of Point Anchored Rockbolts in Circular Tunnels in Elastic Ground

Summary. A simple analytical method for the analysis of point anchored rockbolts is presented in this paper. The solution has been derived for elastic ground and rockbolts, for plane strain conditions, and for tunnels with circular cross section. The method provides accurate results for the rockbolts’ loads and displacements and explicitly includes the connection of the rockbolts to the surrounding ground. The addition of such details to a Finite Element numerical model is critical; otherwise the solution obtained may be dependent on the discretization used and on the stiffness of rockbolts and ground. As an alternative to including details of the rockbolt head and anchor point in the numerical model, which could be computationally very expensive, an equivalent spring constant is proposed. The spring constant is obtained by matching numerical with analytical results for a simple case, but keeping the geometry, material properties, and discretization unchanged.     

A superposition scheme to obtain fundamental boundary element solutions in multi‐layered elastic media

A superposition scheme is proposed to obtain a fundamental solution for boundary elements in multi‐layered elastic media. A three‐layered elastic region is obtained by superposing two sets of bonded half‐planes and subtracting one infinite plane. Therefore, the solution for an element in a layered media can be expressed in terms of bonded half‐plane solutions and an infinite‐plane solution. The major advantages of this superposition scheme are: (1) it is unnecessary to introduce elements at the interface, (2) it can be extended to higher‐order element, and (3) it may be applicable to three dimensions easily. The accuracy and performance of the developed model is illustrated by two examples. For the problem of a pressurized two‐dimensional crack within a three‐layered system, the comparison with other numerical results shows the model is quite accurate and efficient. The model is also used for a study of a practical two‐dimensional mining problem in South Africa, i.e. stoping through a dyke with material properties different from the host rock. Copyright © 2000 John Wiley & Sons, Ltd.     

Heat extraction from aquifer geothermal systems

This study aims to model temperature distributions in an aquifer thermal extraction (ATE) system that contains a single extraction well in a thin confined aquifer. The aquifer is bounded by hot dry rocks with different thermomechanical properties and thicknesses. Based on the heat convection–conduction equation, a mathematical model is developed to describe the spatial and temporal temperature distributions of the ATE systems. The mechanisms of heat transfer in the model involve horizontal convection and thermal conduction in the aquifer, and vertical thermal conduction in both rocks. A semi‐analytical solution in dimensionless form is developed using the Laplace transform technique and its corresponding time‐domain result is computed by the modified Crump method. In addition, the steady‐state solution is obtained by applying the final value theorem. The simulation results from the semi‐analytical solution indicate that the aquifer temperature distributions are affected by aquifer thickness, the thermomechanical properties of the aquifer and rocks, geothermal gradient, outer boundary temperatures of the rocks, extraction rate, and operating time. The present solution can be used as a preliminary tool for assessing heat extraction efficiency in ATE systems. Copyright © 2010 John Wiley & Sons, Ltd.