Comments on “Steady- and transient-state inversion in hydrogeology by successive flux estimation” by P. Pasquier and D. Marcotte

Pasquier and Marcotte [Pasquier P, Marcotte D. Steady- and transient-state inversion in hydrogeology by successive flux estimation. Adv Wat Res 2006;29:1934–52] propose some modifications to the Comparison Model Method (CMM), in order to apply it to transient 3D ground water flow data for conductivity identification. We present some remarks on that paper to improve the comprehension of the basic features of the CMM and of the real value of the novelties introduced by Pasquier and Marcotte.     

Ocean Model Open Boundary Conditions with Volume,Heat and Salinity Conservation Constraints简

Open boundary conditions （OBCs） for a regional ocean model that can be integrated stably over a long timeframe, as well as satisfy the volume, heat and salinity conservation constraints, were developed. First, the idea that the inward and outward flux information can be treated separately in the OBCs was adopted. Second, in order to maintain the property that the volume, heat and salinity remains conserved in the simulation domain, conservation constraints were added to the OBCs, and an inverse method utilized to solve the constraint equations. Ideal experiments were designed to investigate the conservation property, and the OBCs were found to work efficiently to maintain the volume, heat and salinity conservation. It was found that simulations were comparable to observations when the OBCs were applied to a regional ocean model.     

L1-L2范数联合约束稀疏脉冲反演的应用

稀疏脉冲反演实际上就是利用反褶积原理, 从带有噪声的地震道中计算出具有稀疏分布特征的反射系数的振幅和时间.稀疏脉冲反演是非线性优化问题, 通常都是把非线性优化问题转化为线性优化问题, 然后用线性优化算法求解.以范数约束为基础, 提出L1-L2范数联合约束求解的方法.该算法采用了目前国际流行的内点算法, 与传统的优化算法相比, 这种算法具有精度高, 速度快的优点.通过研究人工模型和南海某油田实际数据, 表明该算法的计算结果和测井记录吻合好, 提高了地震分辨率, 目的层段分辨率优于8m.利用反射系数剖面预测的储层厚度和开发井吻合很好, 大大地降低了资源量计算的风险和油田开发的不确定性.      

加权最小范数平差的伪观测解法

首先证明了加权最小范数平差一般解的一个重要性质，由此提出了利用附加的观测值来求解加权最小范数平差问题的方法，并证明了该法与利用一般解法等价性的充要条件。     

Assessing the impact of the hydraulic properties of a crusted soil on overland flow modelling at the field scale

Soil surface crusts are widely reported to favour Hortonian runoff, but are not explicitly represented in most rainfall‐runoff models. The aim of this paper is to assess the impact of soil surface crusts on infiltration and runoff modelling at two spatial scales, i.e. the local scale and the plot scale. At the local scale, two separate single ring infiltration experiments are undertaken. The first is performed on the undisturbed soil, whereas the second is done after removal of the soil surface crust. The HYDRUS 2D two‐dimensional vertical infiltration model is then used in an inverse modelling approach, first to estimate the soil hydraulic properties of the crust and the subsoil, and then the effective hydraulic properties of the soil represented as a single uniform layer. The results show that the crust hydraulic conductivity is 10 times lower than that of the subsoil, thus illustrating the limiting role the crust has on infiltration. Moving up to the plot scale, a rainfall‐runoff model coupling the Richards equation to a transfer function is used to simulate Hortonian overland flow hydrographs. The previously calculated hydraulic properties are used, and a comparison is undertaken between a single‐layer and a double‐layer representation of the crusted soil. The results of the rainfall‐runoff model show that the soil hydraulic properties calculated at the local scale give acceptable results when used to model runoff at the plot scale directly, without any numerical calibration. Also, at the plot scale, no clear improvement of the results can be seen when using a double‐layer representation of the soil in comparison with a single homogeneous layer. This is due to the hydrological characteristics of Hortonian runoff, which is triggered by a rainfall intensity exceeding the saturated hydraulic conductivity of the soil surface. Consequently, the rainfall‐runoff model is more sensitive to rainfall than to the subsoil's hydrodynamic properties. Therefore, the use of a double‐layer soil model to represent runoff on a crusted soil does not seem necessary, as the increase of precision in the soil discretization is not justified by a better performance of the model. Copyright © 2005 John Wiley & Sons, Ltd.     

广义反距离加权空间推估法

提出了k阶广义反距离加权空间推估法。首先,估计空间函数在已知点处的各阶偏导数;然后,利用泰勒级数逼近原理计算待估点处的函数值,根据空间函数偏导数估值的方差协方差评价推估精度;最后,对待估点处的这些推估值进行最小二乘平差。引入了确定k阶广义反距离加权空间推估法阶次k的BIC准则,并以GPS水准推估为例进行了实验。     

Inverse methods and kinetic models of hydrocarbon generation: III. Case histories for autoclave product analysis

From a determination of the transformation matrix for three pyrolysis product experimental data sets, an examination is given of both the applicability of the laboratory experimental data to the modeling of oil cracking in a sedimentary basin, and of the appropriateness of an inverse model. The results of the laboratory experimental data sets, which were done under different thermodynamic conditions and using different sources, show that the transformation matrix varies over each data set and also with time. Therefore, it is necessary to check the data sets before applying them to a basin for hydrocarbon modeling. The laboratory experimental data taken at lower temperature and over longer times appear more pertinent for the construction of an oil-cracking kinetic model suitable for geologic conditions.     

Analytical calculations related to facies simulation

The paper presents calculations related to the volumes of different facies types. Asymptotic distributions are derived approximating the volumes of facies proportions. The number of bodies required to meet net-gross specifications are shown to follow an inverse gaussian distribution. The calculations typically are performed prior to simulations of reservoir facies and may add substantially to the understanding of the model. In particular, analytical results are useful when it comes to assessing parameters of complex simulation models. The theoretical distributions agree well with results based on marked point process simulations of a turbiditic reservoir.     

Geostatistical analysis of field hydraulic conductivity in compacted clay

Hydraulic conductivity(K) of fractured or porous materials is associated intimately with water flow and chemical transport. Basic concepts imply uniform flux through a homogeneous cross-sectional area. If flow were to occur only through part of the area, actual rates could be considerably different. Because laboratory values ofK in compacted clays seldom agree with field estimates, questions arise as to what the true values ofK are and how they should be estimated. Hydraulic conductivity values were measured on a 10×25 m elevated bridge-like platform. A constant water level was maintained for 1 yr over a 0.3-m thick layer of compacted clay, and inflow and outflow rates were monitored using 10×25 grids of 0.3-m diameter infiltration rings and outflow drains subtending approximately 1×1 m blocks of compacted clay. Variography of inflow and outflow data established relationships between cores and blocks of clay, respectively. Because distributions of outflow rates were much less and bore little resemblance to the distributions of break-through rates based on tracer studies, presence of macropores and preferential flow through the macropores was suspected. Subsequently, probability kriging was applied to reevaluate distribution of flux rates and possible location of macropores. Sites exceeding a threshold outflow of 100×10^–9 m/s were classified as outliers and were assumed to probably contain a significant population of macropores. Different sampling schemes were examined. Variogram analysis of outflows with and without outliers suggested adequacy of sampling the site at 50 randomly chosen locations. Because of the potential contribution of macropores to pollutant transport and the practical necessity of extrapolating small plot values to larger areas, conditional simulations with and without outliers were carried out. Simulated scenarios based on all available data compared well with conditional simulations based on randomly chosen locations.This paper was presented at MGUS 87 Conference, Redwood City, California, 14 April 1987.