Adaptive Multiscale Permeability Estimation

With multiscale permeability estimation one does not select parameterization prior to the estimation. Instead, one performs a hierarchical search for the right parameterization while solving a sequence of estimation problems with an increasing parameterization dimension. In some previous works on the subject, the same refinement is applied all over the porous medium. This may lead to over-parameterization, and subsequently, to unrealistic permeability estimates and excessive computational work. With adaptive multiscale permeability estimation, the new parameterization at an arbitrary stage in the estimation sequence is such that new degrees of freedom are not necessarily introduced all over the porous medium. The aim is to introduce new degrees of freedom only where it is warranted by the data. In this paper, we introduce a novel adaptive multiscale estimation. The approach is used to estimate absolute permeability from two-phase pressure data in several numerical examples.     

鲁豫交界地区尾波Q_c值分析研究

该文对濮阳台记录的21次地震进行了地震波衰减研究,讨论了地球介质品质因子随频率变化的关系;利用不同地震采用不同流逝时间的尾波品质因子测定方法,有效去除了由于不同地震或同一地震的地震散射波在纵向介质中穿透深度的不同而导致对品质因子测量结果的影响,从而有利于扑捉到品质因子的异常变化特征与地震孕育过程的关系。利用改进后的模型,对鲁豫交界地区尾波值进行了计算,得到了该地区品质因子随频率的变化关系和区域介质衰减的时间变化特征。     

双重介质中核素迁移模型的解析解及其应用

给出了具有单裂隙的双重介质中，裂隙一端有分时注入源情形下的迁移模型的解析解及一些简化条件下的解析解，并利用解析解研究了核素在双重介质中的迁移问题。     

Bulk-rock Major and Trace Element Compositions of Abyssal Peridotites: Implications for Mantle Melting, Melt Extraction and Post-melting Processes Beneath Mid-Ocean Ridges

This paper presents the first comprehensive major and traceelement data for 130 abyssal peridotite samples from the Pacificand Indian ocean ridge–transform systems. The data revealimportant features about the petrogenesis of these rocks, mantlemelting and melt extraction processes beneath ocean ridges,and elemental behaviours. Although abyssal peridotites are serpentinized,and have also experienced seafloor weathering, magmatic signaturesremain well preserved in the bulk-rock compositions. The betterinverse correlation of MgO with progressively heavier rare earthelements (REE) reflects varying amounts of melt depletion. Thismelt depletion may result from recent sub-ridge mantle melting,but could also be inherited from previous melt extraction eventsfrom the fertile mantle source. Light REE (LREE) in bulk-rocksamples are more enriched, not more depleted, than in the constituentclinopyroxenes (cpx) of the same sample suites. If the cpx LREErecord sub-ridge mantle melting processes, then the bulk-rockLREE must reflect post-melting refertilization. The significantcorrelations of LREE (e.g. La, Ce, Pr, Nd) with immobile highfield strength elements (HFSE, e.g. Nb and Zr) suggest thatenrichments of both LREE and HFSE resulted from a common magmaticprocess. The refertilization takes place in the ‘cold’thermal boundary layer (TBL) beneath ridges through which theascending melts migrate and interact with the advanced residues.The refertilization apparently did not affect the cpx relicsanalyzed for trace elements. This observation suggests grain-boundaryporous melt migration in the TBL. The ascending melts may notbe thermally ‘reactive’, and thus may have affectedonly cpx rims, which, together with precipitated olivine, entrappedmelt, and the rest of the rock, were subsequently serpentinized.Very large variations in bulk-rock Zr/Hf and Nb/Ta ratios areobserved, which are unexpected. The correlation between thetwo ratios is consistent with observations on basalts that D_Zr/D_Hf< 1 and D_Nb/D_Ta < 1. Given the identical charges (5⁺ forNb and Ta; 4⁺ for Zr and Hf) and essentially the same ionicradii (R_Nb/R_Ta = 1·000 and R_Zr/R_Hf = 1·006–1·026),yet a factor of 2 mass differences (M_Zr/M_Hf = 0·511 andM_Nb/M_Ta = 0·513), it is hypothesized that mass-dependentD values, or diffusion or mass-transfer rates may be importantin causing elemental fractionations during porous melt migrationin the TBL. It is also possible that some ‘exotic’phases with highly fractionated Zr/Hf and Nb/Ta ratios may existin these rocks, thus having ‘nugget’ effects onthe bulk-rock analyses. All these hypotheses need testing byconstraining the storage and distribution of all the incompatibletrace elements in mantle peridotite. As serpentine containsup to 13 wt % H₂O, and is stable up to 7 GPa before it is transformedto dense hydrous magnesium silicate phases that are stable atpressures of 5–50 GPa, it is possible that the serpentinizedperidotites may survive, at least partly, subduction-zone dehydration,and transport large amounts of H₂O (also Ba, Rb, Cs, K, U, Sr,Pb, etc. with elevated U/Pb ratios) into the deep mantle. Thelatter may contribute to the HIMU component in the source regionsof some oceanic basalts. KEY WORDS: abyssal peridotites; serpentinization; seafloor weathering; bulk-rock major and trace element compositions; mantle melting; melt extraction; melt–residue interaction; porous flows; Nb/Ta and Zr/Hf fractionations; HIMU mantle sources     

用随机介质模型方法描述孔洞型油气储层

孔洞型储集层具有强烈的非均质性,在空间分布上变化剧烈,规律性差,其地球物理参数(如速度、密度、弹性参数等)在空间上的变化很难用传统的层状介质模型进行刻画.本文借助随机模型方法,研究孔洞型储层中不同孔洞尺度、分布密度与模型参数的关系,建立反映实际孔洞油气储层介质地球物理参数空间分布统计特征的随机介质模型,以更方便地对孔洞型油气储层进行描述,并为碳酸盐岩孔洞型储层地震反射特征的研究提供有效的方法途径.     

A control-volume finite-element method for three-dimensional multiphase basin modeling

In this paper we describe a 3D control-volume finite-element method to solve numerically the coupled partial differential equations (PDEs) governing geological processes involved in the evolution of sedimentary basins. These processes include sediment deposition and deformation, hydrocarbon generation, multiphase fluid flow, and heat transfer in deforming porous media.     

Clogging of saturated porous media by silt‐sized suspended solids under varying physical conditions during managed aquifer recharge

Managed aquifer recharge is an effective method for utilizing excess flood flows, but clogging of porous media is a limiting factor in the implementation of this water storage technique. In recent years, much research on the physical clogging of porous media during artificial recharge has been conducted. However, the understanding of clogging due to silt‐sized suspended solids (SS) is still inadequate, especially under varying physical conditions. Here, we subjected sand columns to controlled rates of flow and SS suspensions to investigate the influence of media size, SS size, SS concentration, and flow velocity on the clogging of porous media by silt‐sized SS. The results show that the diameter ratio of SS particles to sand grains is the dominant factor influencing the position of physical clogging. As pore velocity increased, the mobility of silt‐sized SS was enhanced and retention in the porous media decreased noticeably. The spatial retention profiles in the porous media were found to vary greatly at different flow velocities. The SS concentration of the infiltrating suspension also dramatically influenced the mobility and deposition of silt‐sized SS particles, such that high concentrations accelerated the clogging process. As the different physical factors changed, the breakthrough curves and retention profiles of silt‐sized SS particles changed obviously and the mechanisms of retention differed. On the whole, clogging position is mainly determined by particle size ratio, but clogging rate is dominated by a variety of factors including particle size ratio, SS concentration, and flow velocity.     

Point dilution tests to calculate groundwater velocity: an example in a porous aquifer in northeast Italy

ABSTRACT

The point dilution test is a single-well technique for estimating horizontal flow velocity in the aquifer surrounding a well. The test is conducted by introducing a tracer into a well section and monitoring its decreasing concentration over time. When using a salt tracer, the method is easy and inexpensive. Traditionally, the horizontal Darcy velocity is calculated as a function of the rate of dilution and is based on the simple assumption that the decreasing tracer concentration is proportional both to the apparent velocity into the test section and to the Darcy velocity in the aquifer. In this article, an alternative approach to analyse the results of point dilution tests is proposed and verified using data acquired at a test site in the middle Venetian plain, northeast Italy. In this approach, the one-dimensional equilibrium advection–dispersion equation is inverted using the CXTFIT model to estimate the apparent velocity inside the test section. Analysis of the field data obtained by the two approaches shows good agreement between the methods and suggests that it is possible to use the equilibrium advection–dispersion equation to estimate apparent velocity over a wide range of velocities.

Editor D. Koutsoyiannis; Associate editor K. Heal