与温度全相关的有效骨架压力模型及其在库车坳陷的应用

鉴于世界上极大多数盆地模拟系统的超压史尚未达到实用水平,本文基于力平衡方程、质量守恒法则及达西定律,完整地推导出了与温度全相关的有效骨架压力史的数学模型,使之适用于盆地模拟。通过塔里木盆地库车坳陷的实际模拟,尤其是克拉2井和依南2井超压史的现今值接近于实测数据,即平均相对误差绝对值分别约为3%和5%,证明该有效骨架压力模型具有较高的精度,既能使超压数值模拟避免通过岩石力学实验来测试有效骨架压力,又能达到实用水平。库车坳陷的实例指出:1)在有效骨架压力计算中所用的孔隙流体的密度和黏度与温度相关是必要的,否则超压计算结果会远离实测压力数据。这也恰好与盆地模拟中这两个参数原本与温度相关一致了;2)影响超压的最敏感参数是渗透率,将实测压力数据作为约束条件,在容许的范围内对Kozeny-Carman渗透率公式中的Sa(骨架的颗粒比面)进行校准,这是获得正确的超压史之关键;3)剥蚀事件对地史、热史乃至超压史的模拟结果的影响。     

渤中坳陷中部地区构造应力场光弹模拟实验研究

在对渤中地区沉积构造详细分析的基础上,明确了其中部地区断裂的分布特征,采用以E44-6101型环氧树脂为基体,顺丁烯二酸酐为固化剂,邻苯二甲酸二丁脂为增塑所配制的材料制成渤中中部地质结构模型,对渤中坳陷中部地区构造应力场进行了光弹实验模拟,再现了近北北东向和近东西向断裂在晚近构造格局下的应力状态及分布。试验证明,在构造应力场的作用下,油气往往由高应力值区向低应力值区迁移,而且渤中地区自新近系以来,最大主压应力方向比较稳定。基于试验结果的分析,对北北东向和近东西向两组断裂在晚近应力场的作用下的构造特征及其对油气的控制作用进行讨论。由于渤中坳陷是中国的油气主要产地,通过研究盆地构造应力场的分布,可为盆地的油气勘探和开发提供依据。     

Regional fracture network permeability using outcrop scale measurements

Estimating the hydraulic properties of fractured aquifers is challenging due to the complexity of structural discontinuities that can generally be measured at a small scale, either in core or in outcrop, but influence groundwater flow over a range of scales. This modeling study uses fracture scanline data obtained from surface bedrock exposures to derive estimates of permeability that can be used to represent the fractured rock matrix within regional scale flow models. The model is developed using PETREL, which traditionally benefits from high resolution data sets obtained during oil and gas exploration, including for example seismic data, and borehole logging data (both lithological and geophysical). The technique consists of interpreting scanline fracture data, and using these data to generate representative Discrete Fracture Network (DFN) models for each field set. The DFN models are then upscaled to provide an effective hydraulic conductivity tensor that represents the fractured rock matrix. For each field site, the upscaled hydraulic conductivities are compared with estimates derived from pumping tests to validate the model. A hydraulic conductivity field is generated for the study region that captures the spatial variability of fracture networks in pseudo-three dimensions from scanline data. Hydraulic conductivities estimated using this approach compare well with those estimated from pumping test data. The study results suggest that such an approach may be feasible for taking small scale fracture data and upscaling these to represent the aquifer matrix hydraulic properties needed for regional groundwater modeling.     

Integrative assessment of hydrological, ecological, and economic systems for water resources management at river basin scale

This study presents a basin-scale integrative hydrological, ecological, and economic (HEE) modeling system, aimed at evaluating the impact of resources management, especially agricultural water resources management, on the sustainability of regional water resources. The hydrological model in the modeling system was adapted from SWAT, the Soil and Water Assessment Tool, to simulate the water balance in terms of soil moisture, evapotranspiration, and streamflow. An ecological model was integrated into the hydrological model to compute the ecosystem production of biomass production and yield for different land use types. The economic model estimated the monetary values of crop production and water productivity over irrigated areas. The modeling system was primarily integrated and run on a Windows platform and was able to produce simulation results at daily time steps with a spatial resolution of hydrological response unit (HRU). The modeling system was then calibrated over the period from 1983 to 1991 for the upper and middle parts of the Yellow River basin, China. Calibration results showed that the efficiencies of the modeling system in simulating monthly streamflow over 5 hydrological stations were from 0.54 to 0.68 with an average of 0.64, indicating an acceptable calibration. Preliminary simulation results from 1986 to 1995 revealed that water use in the study region has largely reduced the streamflow in many parts of the area except for that in the riverhead. Spatial distribution of biomass production, and crop yield showed a strong impact of irrigation on agricultural production. Water productivity over irrigated cropland ranged from 1 to 1640 USD/(ha·mm⁻¹), indicating a wide variation of the production conditions within the study region and a great potential in promoting water use efficiency in low water productivity areas. Generally, simulation results from this study indicated that the modeling system was capable of tracking the temporal and spatial variability of pertinent water balance variables, ecosystem dynamics, and regional economy, and provided a useful simulation tool in evaluating long-term water resources management strategies in a basin scale.     

黑河流域生态—水文过程综合遥感观测联合试验总体设计

介绍了"黑河流域生态—水文过程综合遥感观测联合试验"的背景、科学目标、试验组成和试验方案。试验的总体目标是显著提升对流域生态和水文过程的观测能力,建立国际领先的流域观测系统,提高遥感在流域生态—水文集成研究和水资源管理中的应用能力。由基础试验、专题试验、应用试验、产品与方法研究和信息系统组成。其中,①基础试验:搭载微波辐射计、成像光谱仪、热像仪、激光雷达等航空遥感设备,开展一系列关键生态和水文参量的观测;发展遥感正向模型及反演和估算方法。形成覆盖全流域的水文气象综合观测网,为流域生态—水文模型研究提供更有代表性的模型参数、驱动数据及更高精度的验证数据。构建无线传感器网络,度量生态水文模型所需的若干关键的驱动、参数和模型状态的空间异质性。开展航空遥感定标和地基遥感试验。依托传感器网络,并辅之以地面同步和加密观测,开展遥感产品真实性检验。②专题试验:开展"非均匀下垫面多尺度地表蒸散发观测试验",采用密集的涡动相关仪、大孔径闪烁仪与自动气象站的观测矩阵,为揭示地表蒸散发的空间异质性,实现非均匀下垫面地表蒸散发的尺度扩展,发展和验证蒸散发模型提供基础数据。③应用试验:在流域上、中、下游分别开展针对积雪和冻土水文、灌溉水平衡和作物生长、生态耗水的综合观测试验,将观测数据和遥感产品用于上游分布式水文模型、中游地表水—地下水—农作物生长耦合模型、下游生态耗水模型,通过实证研究提升遥感在流域生态—水文集成研究和水资源管理中的应用能力。加强试验将在2012年5月起按中游、上游、下游的顺序展开,全流域持续观测期为2013—2015年。在各类试验的支持下,开展全流域生态—水文关键参量遥感产品生产,发展尺度转换方法,建立多源遥感数据同化系统。     

鄂尔多斯北部盆地古流体动力场的演化特征

以鄂尔多斯北部盆地为例,针对古流体动力场的演化问题,应用三维数值模拟技术,恢复了上古生界不同时期地层厚度和地质、水文地质条件;推导了考虑多个要素的地下流体三维数学模型,并依据不同时期古流体动力场特征,初步分析了上古生界油气的运移聚集规律。结果表明,在鄂尔多斯北部盆地,运用三维数值模拟方法能够比较完善地描述古流体动力场的特征和比较准确地指明油气运移和聚集区带,为天然气勘探靶区的确定提供比较可靠的依据。      

沁水盆地煤样静水压力下渗透率实验及模型分析

为研究沁水盆地煤样渗透率演化规律,构建了煤样渗透率测定的瞬态压力脉冲法实验装置,使用N2和CO2在实验室开展了3种试验条件的渗透率测定,应用Connell模型对实验结果进行分析,并讨论了模型预测值和实验值之间差别的原因。结果表明,（1）在恒定孔压变围压条件下渗透率随有效应力增大而减小;在等有效应力条件下,渗透率随孔压增大而减少;在恒定围压变孔压条件下,随孔压增大,渗透率呈先减小后变大的趋势。（2）运用Connell模型预测的恒定围压变孔压条件的渗透率值大于实验值,原因可能是由于裂隙压缩性系数和吸附应变系数存在估计误差。通过开展实验室渗透率实验和模型分析,对指导实验室内二氧化碳封存和气体驱替实验及其模拟研究具有借鉴意义。     

海水入侵的实验研究

在海岸含水层中,含较多盐分的海水向陆地方向入侵形成淡-咸水交界面,通常称为海水入侵界面。该界面的动态变化受气象、水文地质和人类活动等因素的影响,是海水入侵问题的主要研究内容。本文介绍了海水向含水层中入侵的室内实验研究,以着色海水和图像处理方法获得入侵界面的形态及其动态变化。实验方法先进,所得数据整体性强,精度高,而且观测方法对流场没有干扰,比传统的提取水样进行化验的方法更有效。实验中模拟了海洋潮汐条件,揭示了入侵界面在不同实验条件下的动态特征。应用数值模拟方法揭示了流场的特征。计算表明,流场分两个区———淡水区和咸水区,由入侵界面分开。淡水区的水流相对活跃,淡水由入侵界面以上区域渗出边界,在界面以上形成淡水通道。     

基于MLS的FLAC3D复杂三维地质建模

针对岩土力学分析软件FLAC3D建模效率较低的不足之处,在分析FLAC3D建模命令以及网格模型＊.FLAC3D文件格式特点的基础上,介绍了空间坐标变换、MLS插值的方法与意义,并提出了新的建模方法,开发了相应的计算机辅助建模程序.该程序直接读取原始数据点,无需借助第三方软件,通过坐标变换、MLS插值,可快速生成具有复杂边界的地质模型,且便于材料参数和应力边界条件参数的赋值,可有效地提高工作效率.     

Optimization and rationalization of karst groundwater resource in Feicheng Basin,China

Feicheng Basin is an important energy base in Shandong Province, China. The Ordovician limestone is the main aquifer for water supply in the basin. As a result of the current practice in which individual water users control the exploitation of the karst water, the regional groundwater level has dropped more than 30 m over the last 20 years. The groundwater is overdrawn in approximately 300 km² because of the poor management from 1981. A numerical model with distributed parameters is developed for the karst aquifer to help properly manage the water resource and sustain its utilization. The responding matrix method is used to couple the prediction model and the management model. The objective of the modeling exercises is to obtain the maximum amount of water while keeping the water level drawdown within an acceptable level. The optimum management scheme is selected among two proposed designs. The recommended measures from this study include (1) improvement of water use efficiency in industries, (2) comprehensive utilization of the water pumped from the mines, and (3) improvement of the irrigation practice. When the domestic water supply and corporates water demand take the top priorities, the numerical modeling provides a useful tool to help with the rational distribution of the karst water resources among different water users.     

Immobilization of uranium in the presence of Fe0(s): Model development and simulation of contrasting experimental conditions

A geochemical model is developed for the immobilization of U in the presence of metallic Fe. Zero-valent iron (ZVI) serves as a reducing agent inducing the reductive-precipitation of U, and ZVI corrosion products can serve as absorbing agents for U. The numerical model developed allows the complex interactions of U in solution in differing concentrations to be examined, under variable pH and redox conditions, with or without carbonate, in the presence of ZVI of different size and surface area. It incorporates Fe corrosion, Fe(II) and Fe(III) corrosion product formation, reductive-precipitation of U from the soluble U(VI) valence to the poorly soluble U(IV) valence, adsorption/de-sorption of U onto the Fe oxide corrosion products, and aqueous speciation. The processes of Fe corrosion and reductive precipitation of U are simulated as non-equilibrium, an improvement over other geochemical models. The reductive-precipitation process may use either ZVI or Fe(II) as the reducing agent. The model is calibrated using 3 separate sets of experimental data from published literature that cover a wide range of redox conditions. Sensitivity of the model predictions to variations in input parameters is examined. The simulation results show that the different published experimental results can be explained by different solution chemistries in the studies, specifically O₂ and CO₂ availability and pH, and the amount and surface area of the metallic Fe. With this numerical model the behavior of U in ZVI containing systems over a range of conditions realistic for groundwater can be investigated. By synthesizing observations across several experimental studies, it will lead to a broader understanding of the processes controlling U immobilization under varied geochemical conditions.     

Hybrid analytical and finite element numerical modeling of mass and heat transport in fractured rocks with matrix diffusion

Quantification of mass and heat transport in fractured porous rocks is important to areas such as contaminant transport, storage and release in fractured rock aquifers, the migration and sorption of radioactive nuclides from waste depositories, and the characterization of engineered heat exchangers in the context of enhanced geothermal systems. The large difference between flow and transport characteristics in fractures and in the surrounding matrix rock means models of such systems are forced to make a number of simplifications. Analytical approaches assume a homogeneous system, numerical approaches address the scale at which a process is operating, but may lose individual important processes due to averaging considerations. Numerical stability criteria limit the contrasts possible in defining material properties. Here, a hybrid analytical–numerical method for transport modeling in fractured media is presented. This method combines a numerical model for flow and transport in a heterogeneous fracture and an analytical solution for matrix diffusion. By linking the two types of model, the advantages of both methods can be combined. The methodology as well as the mathematical background are developed, verified for simple geometries, and applied to fractures representing experimental field conditions in the Grimsel rock laboratory.     

层序地层作为沉积盆地识别标志的研究

论文以盆地充填动力学为指导，从层序界面与盆山转换、层序级别与盆地规模、层序格架与盆地类型、层序充填与盆地演化四方面，对层序与盆地的关系进行了初步研究，结果表明层序作为盆地的识别标志主要包括：1）不同成因类型和规模级别的界面作为划分相应成因级别层序的重要标志，他们反映地壳演化过程中的动力转换事件，其中超级、I级、Ⅱ级界面分别反映盆地域之间、盆地之间、盆山之间的应力转换机制，Ⅲ-Ⅳ级界面揭示了构造、海平面、气候、物源联合引起的盆地容纳空间和沉积机制的转变过程；2）作为板间、板缘、板内过程的沉积记录，层序与盆地之间存在特定的耦合关系：超级层序相当于沉积盆地域，I级层序相当于沉积盆地，Ⅱ级层序对应于盆地构造沉积幕，Ⅲ级层序对应于盆地充填幕，Ⅳ-Ⅴ级层序相当于盆地充填韵律；3）不同成因沉积盆地，层序成因格架及内部构型各异，因而根据层序特征研究，可判识盆地成因属性；4）在盆地演化的不同阶段，盆地性质、构造特征、沉积格局存在差异，所形成的层序构型、界面属性及组合特征等亦不同，因而通过层序充填序列的研究，可揭示盆地演化过程。     

A laboratory study of tracer tomography

A tracer tomographic laboratory study was performed with consolidated fractured rock in three-dimensional space. The investigated fractured sandstone sample was characterized by significant matrix permeability. The laboratory transport experiments were conducted using gas-flow and gas-tracer transport techniques that enable the generation of various flow-field patterns via adjustable boundary conditions within a short experimental time period. In total, 72 gas-tracer (helium) tests were performed by systematically changing the injection and monitoring configuration after each test. For the inversion of the tracer breakthrough curves an inversion scheme was applied, based on the transformation of the governing transport equation into a form of the eikonal equation. The reliability of the inversion results was assessed with singular value decomposition of the trajectory density matrix. The applied inversion technique allowed for the three-dimensional reconstruction of the interstitial velocity with a high resolution. The three-dimensional interstitial velocity distribution shows clearly that the transport is dominated by the matrix while the fractures show no apparent influence on the transport responses.     

Applying a modified conduit flow process to understand conduit-matrix exchange of a karst aquifer

Due to the high heterogeneity and complexity of water flow movement for multiple karst water-bearing mediums,the evaluation,effective development,and utilization of karst water resources are significantly limited.Matrix flow is usually laminar,whereas conduit flow is usually turbulent.The driving mechanisms of water exchange that occur between the karst conduit and its adjacent matrix are not well understood.This paper investigates the hydrodynamic characteristics and the mechanism of flow exchange in dual water-bearing mediums(conduit and matrix)of karst aquifers through laboratory experimentation and numerical simulation.A karst aquifer consisting of a matrix network and a conduit was proposed,and the relationship between the water exchange flux and hydraulic head differences generated from the laboratory experiments was analyzed.Two modes of experimental tests were performed with different fixed water level boundaries in the laboratory karst aquifer.The results indicate that the water exchange capacity was proportional to the square root of hydraulic head differences.The linear exchange term in the conduit flow process(CFP)source program was modified according to experimental results.The modified CFP and the original CFP model experimental data results were compared,and it was found that the modified CFP model had better fitting effects.These results showed that the water exchange mechanism between conduit and matrix is very important for solid-liquid interface reaction,water resource evaluation,and understanding of karst hydrodynamic behavior.     

基于局部重力场建模的Tikhonov正则化点质量核径向基函数方法

针对点质量核径向基函数应用于局部重力场建模中的设计矩阵严重病态问题,本文引入Tikhonov正则化方法对传统点质量核径向基函数方程进行改造,建立了相应的正则化模型。通过模拟数据进行仿真实验,以传统格网化方法作为对比试验,利用"标靶法"确定两种模型的最优结构。实验结果表明:正则化点质量核径向基函数可以直接利用离散数据进行局部重力场建模。在两种模型的最优结构下,当实测数据无污染时,正则化方法达到与传统格网化方法相当的精度;当实测值中加入3 mGal的高斯白噪声时,正则化方法的精度获得了27.9%的提升。这说明本文方法可以应用于局部重力场建模中,且模型结构更优,抗干扰能力更强。     

淮河中游洪水演进模型

针对淮河流域的特点,建立了以水动力学方法为主的河道一维、行洪区二维及对无实测资料地区采用水文学方法为辅的流域洪水演进数学模型.并提出了河道一维、行洪区二维耦合的求解方法.对占相当大比重的无资料地区径流流入干流作了技术处理,这对模型的成功建立起着重要作用.通过9次洪水的验算表明,模型具有较高的精度、较强的通用性、灵活性和实用性.为淮河流域的防洪规划提供了工具.模型的主要思路同样适用于其它流域.     

An efficient matrix-free algorithm for the ensemble Kalman filter

In this work, we present an efficient matrix-free ensemble Kalman filter (EnKF) algorithm for the assimilation of large data sets. The EnKF has increasingly become an essential tool for data assimilation of numerical models. It is an attractive assimilation method because it can evolve the model covariance matrix for a non-linear model, through the use of an ensemble of model states, and it is easy to implement for any numerical model. Nevertheless, the computational cost of the EnKF can increase significantly for cases involving the assimilation of large data sets. As more data become available for assimilation, a potential bottleneck in most EnKF algorithms involves the operation of the Kalman gain matrix. To reduce the complexity and cost of assimilating large data sets, a matrix-free EnKF algorithm is proposed. The algorithm uses an efficient matrix-free linear solver, based on the Sherman–Morrison formulas, to solve the implicit linear system within the Kalman gain matrix and compute the analysis. Numerical experiments with a two-dimensional shallow water model on the sphere are presented, where results show the matrix-free implementation outperforming an singular value decomposition-based implementation in computational time.     

Flow balances in St. Andrew Bay revealed through hydrodynamic simulations

The circulation patterns in St. Andrew Bay, Florida are revealed through the application of a well-tested, extensively used three-dimensional hydrodynamic model. A high resolution grid resolving both the horizontal and vertical directions is used with a systematically developed set of forcing functions to simulate conditions over a full year. Water levels at the three open boundaries are deduced from a year-long deployment of pressure gauges, and freshwater loadings are based upon drainage basin characteristics and precipitation measurements. Model validation involves comparisons with hydrographic casts taken at twelve stations distributed throughout the bay at monthly intervals. The relative average error between the observed and model-predicted salinity is 15% for the surface of the water column and 4% for the bottom. The annual net flow balance consists of an influx of water at the two Intracoastal Waterway open boundaries, with that water exiting to the Gulf of Mexico. An average of about 100 m³ s⁻¹ enters from East Bay and about 40 m³ s⁻¹ enters through West Bay. On shorter time scales, the flow balance is quite variable both in terms of magnitude and direction. This study also presents methods to overcome the paucity of data that is usually available for the development of such a model. These include techniques to take bottom pressure data sets with short gaps and create reliable sea surface elevation boundary conditions and to take precipitation data and drainage basin characteristics and produce estimates of freshwater inflows.