地下水非稳定流的LTFAM算法

渗流域内应用拉普拉斯变换(LT)建立相应的有限分析(FAM)方程,顾及渗流域内地下水流的初始条件和边界条件,可在LT空间构成一个封闭的以水头像函数为变量的线性方程组。将此方程组所得的解,通过Stehfest数值反演公式,可归化为时间域的解(水头)。由于时间t被隐含在数值方程内,从而克服了传统数值法按时段(△t)逐步迭代的缺陷,提高了计算效率,也为用嵌入法建立地下水流管理模型提供了一条捷径。     

地铁施工对济南白泉泉群流量的影响

为研究地铁建设对济南白泉泉群的影响，在综合分析白泉泉域地质、水文地质条件的基础上，假定研究区岩溶强径流带位置及水力性质，利用FEFLOW软件建立地下水流数值模型。以规划地铁M1号线为研究对象，分析了济南东站、梁王站、梁王东站分别施工及3个站同时施工4种情景下，采用施工降水或施工降水+人工回灌两种施工方式对白泉泉群流量的影响。结果表明：单独采用施工降水的施工方式使得白泉泉群流量衰减，其中3个站同时施工对泉流量的影响最大，泉流量最大衰减达5.48%；各站分别施工时，济南东站对泉流量影响最大，泉流量较未施工时减少了0.043×10⁴ m³/d。采用施工降水+人工回灌的施工方式，能够有效缓解泉流量的衰减，各车站施工时的泉流量衰减由仅施工降水时的2.26%~5.48%降低至0.08%~1.21%。岩溶强径流带有利于地下水形成优势径流，促进白泉泉群补给，一定程度上缓解因地铁施工引起的泉流量衰减。     

漳州热田的对流热流和传导热流的研究

漳州地热系统属对流型地热系统.漳州热田是我国东南沿海地区目前所见热田中温度最高（121.5℃）的一个.地表热流值的研究表明,热田中心具有最大的实测热流值（359mw/m²）.本文根据热田内152个钻孔的测温资料和56块岩石样品的热导率数据,采用三种不同方法计算出漳州热田及其邻近地区的大地热流值,并讨论了热流值分布的特点.     

云南大地热流及地热地质问题

本文报道了云南20多个大地热流测试数据。据热流分布特点并参考前人有关资料,结合地球物理测深资料、岩石热物性与放射性生热率测试数据,对各区地壳深度范围的温度、热流配分组成作了粗略分析,划分出各有特征的若干地质-热流区。云南普遍地热偏高(滇东南小范围除外)的主要原因是新构造运动强烈;怒江以西近代岩浆活动亦强烈,具有典型的板块交汇带的基本地热特点     

Fracture analysis in the south-western Corinth rift (Greece) and implications on fault hydraulic behavior

This paper reviews the data concerning the fracture network and the hydraulic characteristics of faults in an active zone of the Gulf of Corinth. Pressure gap measured through fault planes shows that in this area the active normal faults (Aigion, Helike) act, at least temporarily and locally, as transversal seal. The analysis of the carbonate cements in the fractures on both the hangingwall and the footwall of the faults also suggests that they have acted as local seals during the whole fault zone evolution. However, the pressure and the characteristics of the water samples measured in the wells indicate that meteoric water circulates from the highest part of the relief to the coast, which means it goes through the fault zones. Field quantitative analysis and core studies from the AIG-10 well have been performed to define both regional and fault-related fracture networks. Then laboratory thin section observations have been done to recognize the different fault rocks characterizing the fault zone components. These two kinds of approach give information on the permeability characteristics of the fault zone. To synthesize the data, a schematic conceptual 3D fluid flow modeling has been performed taking into account fault zone permeability architecture, sedimentation, fluid flow, fault vertical offset and meteoric water influx, as well as compaction water flow. This modeling allows us to fit all the data with a model where the fault segments act as a seal whereas the relays between these segments allow for the regional flow from the Peloponnese topographic highs to the coast.     

Statistical evaluation of methods for the calculation of static formation temperatures in geothermal and oil wells using an extension of the error propagation theory

The feasibility of using the more sophisticated weighted least-squares (WLS) model, as opposed to the traditional ordinary least-squares (OLS), in linear regressions of BHT data to estimate the static formation temperatures (SFT) was investigated. The most commonly used analytical methods (line-source; spherical and radial heat flow; and cylindrical heat source) were applied. Error propagation equations were derived to calculate errors in the time function of each method. These errors were combined with the BHT measurement errors to compute weighting factors for applying the WLS. Intercept uncertainties were estimated for all regressions using sets of synthetic and actual borehole logs taken from geothermal and oil applications. SFT computed with the spherical and radial heat flow method were generally greater than those from the other two methods.     

Temperature–composition–time (T–X–t) data from authigenic K-feldspar: An integrated methodology for dating fluid flow events

Integration of fluid inclusion analysis with high spatial resolution Ar–Ar dating of K-feldspar cements has been used to resolve and reconstruct palaeo-fluid flow. Fluid inclusion analysis allows discrimination of distinct cement phases, thereby identifying discrete episodes of fluid flow. Ar–Ar dating of the same cements via high spatial resolution laserprobe establishes absolute age constraints on the framework previously constructed. Integration of these two datasets yields temperature–composition–time data.     

Reservoir description using dynamic parameterisation selection with a combined stochastic and gradient search

There is a correspondence between flow in a reservoir and large scale permeability trends. This correspondence can be derived by constraining reservoir models using observed production data. One of the challenges in deriving the permeability distribution of a field using production data involves determination of the scale of resolution of the permeability. The Adaptive Multiscale Estimation (AME) seeks to overcome the problems related to choosing the resolution of the permeability field by a dynamic parameterisation selection. The standard AME uses a gradient algorithm in solving several optimisation problems with increasing permeability resolution. This paper presents a hybrid algorithm which combines a gradient search and a stochastic algorithm to improve the robustness of the dynamic parameterisation selection. At low dimension, we use the stochastic algorithm to generate several optimised models. We use information from all these produced models to find new optimal refinements, and start out new optimisations with several unequally suggested parameterisations. At higher dimensions we change to a gradient-type optimiser, where the initial solution is chosen from the ensemble of models suggested by the stochastic algorithm. The selection is based on a predefined criterion. We demonstrate the robustness of the hybrid algorithm on sample synthetic cases, which most of them were considered insolvable using the standard AME algorithm.     

Modèles géostatistiques de concentrations ou de débits le long des cours d'eau

Estimating concentrations or flow rates along a stream network requires specific models. Two classes of models, recently proposed in the literature, are generalized, to the intrinsic case in particular. We present a global construction by ‘streams’, i.e. on the whole set of paths between sources and outlet. Combining stationary or intrinsic one-dimensional random functions leads to stationary or intrinsic models on segments, with discontinuities at the forks. A construction from outlet to sources, leads to stationary or intrinsic models on each stream, without any discontinuity at the forks. The linear variogram is found as a particular case. The extension to the linear model of coregionalization is immediate, allowing a multivariate modelling of concentrations. To cite this article: C. de Fouquet, C. Bernard-Michel, C. R. Geoscience 338 (2006).