三种混合泥质砂岩有效介质通用电阻率模型比较

复杂泥质砂岩油气层的识别与评价一直是测井解释领域亟待解决的难点和热点问题,因此,针对复杂泥质砂岩油气层的特点,有必要研究能够描述这类储层导电规律的新导电理论。针对用于复杂泥质砂岩油气层解释的3种混合泥质砂岩有效介质电阻率模型,从理论和实验角度,进行了全面系统的比较。通过对3种有效介质通用电阻率模型建模理论基础的分析,可知混合泥质砂岩有效介质非对称模型基于HB方程和并联导电理论,对称模型基于对称导电理论和并联导电理论,而孔隙结合模型基于混合介质电导率积分方程;对称电阻率模型可用于多连续项存在的地层,而非对称电阻率模型和孔隙结合电阻率模型只允许水为连续项。对骨架不导电的纯砂岩进行模型化简,可知非对称电阻率模型和孔隙结合电阻率模型可以简化成阿尔奇公式,而对称电阻率模型能化简成类似阿尔奇公式。通过1组混合泥质砂岩人造岩样和2组分散泥质砂岩岩样以及1组骨架导电人造岩样的实验数据,验证3种有效介质通用电阻率模型均适用于混合泥质砂岩地层解释,且精度相近。     

裘布依公式在实际应用中的缺陷分析

主要通过对裘布依公式中：1）实际流量（Q）-降深（s）曲线与理论的差异,2）实际井径对流量的影响与理论的差异,3）井中最大流量,4）影响半径等问题的分析,指出裘布依公式的缺陷,并分别提出了一定的解决方案。同时还结合土力学的相关知识,从“达西定律的适用范围”,“渗透系数的不唯一性”,“地基沉降导致含水层变形”三个方面分析了公式的不足之处,更全面地对问题进行了剖析。     

深层岩体松动爆破中不耦合装药效应的探讨

基于冲击波在交界面两侧压力和速度必须各自相等的连续性条件,求解爆轰产物中适用的反射波方程和介质中适用的冲击波方程,得到药包周围介质中冲击波的初始参数.通过波的传播机理,把集中药包应力波随距离的衰减公式扩展到延长药包,并计算耦合与不耦合装药爆破时距爆心相同距离处岩石中冲击波的参数.由计算结果可知:(1)耦合装药爆破时形成的冲击波压力超过岩石抗压强度极限几十倍以上,药包周围岩石形成粉碎区;(2)与耦合装药爆破相比,不耦合装药爆破可以降低孔壁处岩石中冲击波初始压力,但可以增加孔壁后岩石中的冲击波压力.合理的不耦合系数,可使岩石不形成粉碎区,大幅度减少能量耗散;(3)水一般认为是非线性弹性介质,因此水介质成为炸药爆轰产物与岩石间的弹性缓冲层,增加了能量传递,延长了冲击波作用时间,加大了爆炸的作用范围.     

基于REOF的两种卫星降水产品(IMERG和MSWEP)在金沙江流域降水分区尺度精度评估

高时空分辨率降水产品的精度评估是卫星降水用于水文气象干旱等研究的前提。本研究提出在降水分区尺度下评估IMERG和MSWEP两种卫星降水产品的精度,并与不分区尺度(即流域尺度)进行比较。首先采用旋转经验正交函数(REOF)对金沙江流域(JSB)进行降水分区,通过贡献率得出8个分区较为适合。然后识别降水的空间分布特征,发现2种降水产品都可以很好地捕捉降水呈现出的从上游到下游逐渐增加的趋势。最后在日尺度、降水发生概率和极端降水探测能力3个方面对降水产品在分区尺度和不分区尺度的性能进行评估。结果表明,在日尺度上,MSWEP的精度在多数降水分区优于IMERG,被推荐5次(1、3、6、7和8区),集中在流域的中游。同时流域尺度也推荐MSWEP。在降水事件发生概率方面,MSWEP能再现不同等级降水强度的概率密度分布,但过高估计0.1~1 mm/d降水事件的发生概率;而IMERG过高估计小于0.1 mm/d降水事件的概率。在极端降水探测能力方面,流域尺度的KGE值都是正值,且IMERG优于MSWEP,但分区尺度上,KGE值在部分降水分区中存在负值,表明IMERG和MSWEP均不能很好地探测出该区的极端降水事件。本研究成果表明降水分区尺度是必需的,能够更加精细地评估降水产品。研究结果可为具有类似气候条件的卫星降水评估提供参考。     

应用遥感方法估算区域蒸散量的制约因子分析

遥感方法估算区域蒸散量应用比较广泛的有两种,一种是完全以地表热量平衡方程为基础．用遥感方法估算出净辐射、土壤热流量和显热通量,然后用余项法求出蒸散量;另一种是以Penman—Monteith方程为基础．结合地表热量平衡方程,直接估算出蒸散量。本文根据国内外的研究现状,对应用遥感方法估算区域蒸散量的制约因子进行了深入的分析,这些制约因子包括(1)图像信息源;(2)反照率、比辐射率和表面温度等地表参数的遥感反演精度;(3)空气动力学阻抗和表面阻抗模型;(4)估算结果的验证方法(5)时间尺度的扩展问题。分析结果表明：长期以来由于这些因子中所涉及到的许多关键技术和难点问题都没得到很好的解决,因而极大地制约着应用遥感方法估算区域蒸散量的发展。通过对制约因子的分析．可以有助于高精度遥感蒸散模型的建立。     

Short Communication: A numerical method for diffusive transport with moving boundaries and discontinuous material properties

A method for the numerical simulation of diffusive transport with moving boundaries is developed and tested. The variable domain is mapped onto a fixed region, which introduces a term of convective form to the transformed governing equation. The resulting convection/diffusion equation is solved by a finite-difference method. An ‘Immersed Interface’ Method (IIM) is introduced in order to retain second-order accuracy near discontinuities in material properties, where the solution is not smooth. The method performs well in benchmark calculations against an analytical solution. The IIM scheme is capable of treating a strong discontinuity in the gradient, and it is readily extended to two or three dimensions. The methods are illustrated through a calculation for the temperature profile in a growing continental ice sheet, in which the thermal properties are discontinuous at the rock/ice interface. © 1997 by John Wiley & Sons, Ltd.     

Roughness coefficients and velocity estimation in well-inundated sheet and rilled overland flow without strongly eroding bed forms

Even with the flow of water over a soil surface in which roughness elements are well inundated, and in less erosive situations where erosional bed forms are not pronounced, the magnitude of resistance coefficients in equations such as those of Darcy–Weisbach, Chezy or Manning vary with flow velocity (at least). Using both original laboratory and field data, and data from the literature, the paper examines this question of the apparent variation of resistance coefficients in relation to flow velocity, even in the absence of interaction between hydraulics and resulting erosional bed forms. Resistance equations are first assessed as to their ability to describe overland flow velocity when tested against these data sources. The result is that Manning's equation received stronger support than the Darcy–Weisbach or Chezy equations, though all equations were useful. The second question addressed is how best to estimate velocity of overland flow from measurements of slope and unit discharge, recognizing that the apparent flow velocity variation in resistance coefficients is probably a result of shortcomings in all of the listed resistance equations. A new methodology is illustrated which gives good agreement between estimated and measured flow velocity for both well-inundated sheet and rill flow. Comments are given on the predictive use of this methodology. Copyright © 1999 John Wiley & Sons, Ltd.     

Control law for variable damping device defined by a non-linear differential equation

This paper proposes a non-linear control law for a variable damping device (VDD) aimed at reducing structural seismic responses. The VDD is attached to the structure by an auxiliary spring element composing a non-linear Maxwell element. The VDD's damping coefficient is adjusted to control the reactive internal force in the non-linear Maxwell element. A large controlled force is thus produced with little external power required to adjust the VDD's damping coefficient. The proposed control law defines the rate or increment of the VDD's damping coefficient at a certain moment by a differential equation or its discretized form. The controlled force vs. deformation relation plots parallelogram-like hysteretic curves, which indicates quick action and energy dissipation. Fundamental characteristics of an SDOF model with the VDD controlled by the proposed law are examined for impulse, sin and seismic excitations. The law for the SDOF model is extended to one for an MDOF model. The control effect for a 3DOF model is examined by numerical experiments. Copyright © 1999 John Wiley & Sons, Ltd.