河南义马地区页岩气储层孔隙非均质性

运用高压压汞、液氮吸附及渗透率测试等实验,利用极差、突进系数、变异系数等参数,表征煤系页岩孔渗特征及孔隙层间非均质性,分析孔隙层间非均质性主控因素。结果表明:河南义马地区上、下石盒子组泥页岩微孔与小孔较为发育,山西组大孔较为发育,分别提供了气体吸附附着面积和储存运移空间;上石盒子组孔容与孔径相对偏差最小,孔隙分布均匀,山西组比表面积相对偏差最小,表面积分布均匀;随围压增大,渗透率不断降低,且满足负指数相关关系;渗透率级差、突进系数和变异系数显示下石盒子组渗透率非均质性较弱,山西组渗透率非均质性很弱,更易于压裂开发。非均质性宏观上主要受沉积物质组成和构造改造作用影响,微观上受成岩演化影响。     

基于集合卡尔曼滤波的多相流模型参数估计——以室内二维砂箱中重质非水相污染物入渗为例

重质非水相有机污染物(DNAPL)泄漏到地下后,其运移与分布特征受渗透率非均质性影响显著。为刻画DNAPL污染源区结构特征,需进行参数估计以描述水文地质参数的非均质性。本研究构建了基于集合卡尔曼滤波方法(EnKF)与多相流运移模型的同化方案,通过融合DNAPL饱和度观测数据推估非均质介质渗透率空间分布。通过二维砂箱实际与理想算例,验证了同化方法的推估效果,并探讨了不同因素对同化的影响。研究结果表明:基于EnKF方法同化饱和度观测资料可有效地推估非均质渗透率场;参数推估精度随观测时空密度的增大而提高;观测点位置分布对同化效果有所影响,布置在污染集中区域的观测数据对于参数估计具有较高的数据价值。     

考虑空间效应的地下洞室爆破开挖松动区参数场位移反分析

将地下洞室爆破开挖松动区视为一个随开挖过程演变的非均匀、非稳定三维扰动场,松动区内岩体力学参数则是一个具有时空演化特性的参数场。考虑爆破开挖扰动空间效应和岩体真实变形响应,提出了真实工作状态下开挖松动区岩体参数场辨识的位移反分析方法。基于局部监测变形空间插补得到的空间位移场,通过分析洞室爆破开挖围岩变形扰动机制,建立了开挖松动区岩体变形模量参数场数值演化模型,并进行了模型适用性和参数敏感性分析。在此基础上,以变形模量参数为例,结合围岩实测位移信息,提出了开挖松动区参数场位移反分析的动态实现过程。将该方法应用于溪洛渡地下洞室群施工期参数场反演和围岩稳定动态反馈评价及预测,结果表明,该方法合理有效,在大型地下洞群施工开挖与快速监测反馈方面具有显著的工程适用性及实用性。     

沿程渗漏河道的洪水流量演算模型

从质量守恒原理出发，推导了明渠非恒定流在考虑入渗条件下的连续方程。将水力学原理与水文学方法结合起来，研制了适应海河南系平原河道现状条件的有沿程渗漏的明渠非恒定流河道洪水流量演算模型，并对海河地系７个河段９次洪水过程进行了模型参数的率定检验，较好地解决了北方平原河道沿渗漏的洪水流量演算问题。     

浅层高渗砂岩油藏分级相控地质建模

由于裂谷盆地浅层疏松砂岩储层非均质性强,渗透率级差变化范围巨大,常规相控建模无法表征储层内部非均质性,导致渗透率模拟结果与测井及试井解释难以匹配。为了解决这一难题,在沉积微相研究的基础上,利用流动分层指标(FZI)、孔隙喉道半径(R35)2个参数,将储层细分为4类流动单元,依据各流动单元的孔渗对应关系重新计算测井渗透率。采用分级相控建模的思路,首先建立沉积微相模型,在微相模型的控制下,建立流动单元模型。在此基础上,根据各流动单元的孔渗分布特征,建立各流动单元的孔隙度和渗透率场。该方法建立的三维地质模型,既精细刻画了储层内部非均质性,又能定量表征储层物性空间分布特征,保证渗透率模拟结果与测井及试井结果一致,为油藏数值模拟提供了准确的数据基础。     

非均质综合指数法在砂砾岩储层非均质性研究中的应用——以双河油田V下油组为例

影响储层非均质性形成与分布的因素较多,若只从某一个或某几个因素去考虑,可能会影响对储层非均质程度的认识(宏观和微观),目前多局限于从静态或动态角度研究储层的非均质性,而将动态和静态资料相结合的研究较少.为了充分考虑影响储层非均质性分布的因素,在研究双河油田V下油组储层的非均质性时,采用了主因子分析法选出了沉积微相、砂体厚度、孔隙度、渗透率、油层的顶底面微构造、目前的综合含水率和流动带指标7个参数作为反映储层非均质性的参数.利用波的叠加原理在平面上对每个参数进行叠加从而得出了非均质综合指数值.现场实践证明,其结果比较真实地反映了储层非均质性的地下分布状况,符合率在90%以上.     

犁底层土壤入渗参数的空间变异性

通过野外布点试验,应用双环入渗仪测得了犁底层土壤各点的稳渗率fc和入渗经验指数α.试验结果分析表明,fc和α具有较大的空间变异性,并服从对数正态分布.根据fc和α在该田块内的等值线分布图,分析了入渗参数的空间分布特征.     

非均质介质中重非水相污染物运移受泄漏速率影响数值分析

重非水相污染物(DNAPL)在地下介质中运移和分布受多种因素控制,包括DNAPL本身的物理化学性质,土的性质,泄漏条件等等。由于介质的非均质性,使得多相流运移行为更为复杂。基于地下水随机理论构建渗透率随机场,采用蒙特卡罗方法探讨泄漏速率对非均质饱和介质中DNAPL运移的影响。数值结果表明,在泄漏总量一定的情况下,泄漏速率越低,介质非均质性对DNAPL运移的影响程度越高。反之,DNAPL的渗漏速率越高,小尺度地层的非均质性影响越低。由于DNAPL运移过程中在垂直方向受重力的影响,污染羽在空间上的质心位置(一阶矩)以及展布范围(二阶矩)在垂直方向上的变异程度要高于水平方向。     

地下洞室群对地面运动影响问题的级数解答:SH波入射

采用波函数展开法给出平面SH波入射下半空间中洞室群对地面运动影响问题的一个级数解答，并分析了入射波频率和入射角以及两洞室之间距离等参数对地面运动的影响。数值结果表明，地下洞室群对沿线地面运动具有显著的放大作用：两个洞室情况地面运动可以达到单个洞室情况的1．6倍、无洞室情况的3．5倍以上。建议地铁等地下工程在设计和规划时应考虑工程建设后对沿线设计地震动的影响。     

非均质裂隙介质中重非水相流体运移

由于裂隙介质具有强烈非均质性,使得重非水相流体(DNAPLs)在裂隙介质中的运移行为较孔隙介质更为复杂.基于指示模拟算法模拟裂隙介质的二元变量(粗糙面接触和裂隙开口),综合模拟退火算法构建裂隙渗透率随机场.采用T2VOC模拟DNAPLs在裂隙介质的运移,探讨粗糙面接触的相关长度、各向异性比、倾角以及非均质性程度对DNAPLs运移分布的影响.数值分析结果表明,粗糙面接触的相关性越差,DNAPLs污染范围越大;粗糙面接触的各向异性比和倾角增大,将导致更多DNAPLs残留;而相关长度和各向异性比增大以及倾角减小,都会引起DNAPLs锋面运移速率增大;随着渗透率非均质性增强,会导致局部粗糙面接触上蓄积的DNAPLs饱和度增大,运移速率减小.     

Geochemistry of Thermal Waters of Continental Margin of Far East of Russia

Studied waters belong to warm(T=30-50℃),alkaline(pH=8.9-9.3),low mineralized(TDS235 mg/1)Na-HCO_3 or Na-SO_4-HCO_3 thermal waters with high content of SiO_2(up to 81 mg/l)and F(up to 3.9 mg/1),occur on modern volcano-tectonic rejuvenated areas of Eastern Sikhote-Alin orogenic belt.Low~3He concentration as well as N_2/O_2 and N_2/Ar ratios exclude influence of deep mantle fluid.New rare earth element data constrain our understandmg of water-rock interaction occurring in the water source region.Meteoric origin of waters is proved by stable isotope values varying from-71‰to-136.1‰and from-10.8‰to-18.8‰forδ~2U andδ~(18)O respectively.REE patterns reflect high pH,resultfing from water-rock interaction and oxidative conditions.Calculations of deep aquifer temperature using Na-K and quartz geothermometers show 116.8-131.1°C and 82.2-125.8℃respectively.Presence of deep faults both with abnormal thermal gradient(~45-50 K/km)define unique geochemical shape of thermal waters of Sikhote-Alin,area,where no present volcanic activity is registered.     

黄土路基沉降曲线特征分析

在分析郑西客运专线试验段路基沉降实测资料的基础上,提出沉降变形常规预测方法受观测时间及单次观测结果的影响比较大,不利于用来预测最终沉降。通过分析现场沉降变形与时间关系曲线,提出路基沉降的发展主要受固结的影响,其发展是一个从沉降增加到趋缓的过程,可以明显分为固结变形及流变变形,并用非饱和土等效固结理论加以解释,提出了一种利用s- 曲线特征预测最终沉降量的方法,经检验对比发现其预测结果受观测时间及单次观测结果的影响比较小,优于其它常规预测方法     

沉管隧道基槽爆破施工对既有堤岸稳定性影响的数值仿真分析

结合爆破数值模拟中爆破荷载的特点,分析了各种荷载施加方法的特点及爆破震动场模拟的要求,基于圣维南荷载等效原理,提出了爆破荷载的等效施加方法,即将作用在炮孔壁上的爆破荷载等效后,施加在同排炮孔中心线或者面上来模拟爆破荷载的作用。结果表明,等效施加方法在近区存在一定差异,但在中远区吻合得较好。利用该等效荷载施加方法,采用三维实体模型计算分析了某沉管隧道基槽爆破施工对既有堤岸稳定性影响,计算了周围不同位置、不同深度位置处的振动响应,并在此基础上,分析研究了减小爆破振动对既有结构物影响的减振措施,可为类似工程提供借鉴。     

Structure of sediments of kaolinite

In an assembly of clay particles placed in a fluid, each particle is typically subjected to: (1) double-layer repulsive forces; (2) van der Waals attractive forces; and (3) contact mechanical forces. The study presented here concerns an approximate, quantitative analysis of clay suspensions, with considerations to the first two - the physico-chemical forces. Using recent theories to calculate the physico-chemical forces between two clay particles in an approximate model of an assembly, the equilibrium void ratio of a clay suspension is computed. The mechanical forces are ignored in the analysis. The results serve to verify the validity of physico-chemical theories employed and help interpret experimental data more fundamentally in terms of the system variables.     

不同形式诱导缝的等效强度承载性能的数值分析

将考虑尺寸效应规律的混凝土等效断裂模型应用到不同形式诱导缝等效强度公式中,得到了修正后的不同形式诱导缝的等效强度公式。同时,得到了诱导缝等效强度的变化曲线,并比较了单双向间隔诱导缝的优劣,得到双向间隔诱导缝优于单向间隔诱导缝的结论,为高拱坝中诱导缝的设计提出了建议。     

土的可恢复剪胀的一种解释

从土的各向异怀角度对土的可恢复剪胀现象进行了解释。基于各向异性情况下的土体弹性本构关系理论分析，认为土的可恢复剪胀现象可部分归因于土的各向异性引起的弹性剪胀。借助有关土体弹性参数实验结果，研究了应力诱导各迥异性对土体弹性剪胀的影响，结果表明：随土体应力诱导各向异性的增大，土体的弹性剪胀也增大。从土体弹性剪胀角度研究了土的卸荷体缩条件，认为土体卸荷体缩取决于加载应力路径的应力增量比，给出了土体出现卸荷体缩的区域。     

岩石坚硬程度对核电厂地基-基础地震响应特征的影响分析

考虑核电厂地基-基础的动力相互作用,应用显式动力有限差分法分析了地震作用下极软岩、较软岩、坚硬岩上核电厂建筑结构基础的地震响应特征,比较了岩石坚硬程度对基础加速度反应谱的影响。研究表明：随着岩石坚硬程度的提高,核电厂建筑物结构基础的地震响应有增加的趋势;在周期轴上,基础处的加速度反应谱曲线会随着岩石坚硬程度的提高逐渐向短周期（高频段）方向移动。在高频段,建造于较坚硬岩石上基础结构的加速度反应谱值偏大;在中等频段,建造于较软岩石上基础结构的加速度反应谱值偏大;在低频段,岩石坚硬程度对加速度反应谱的影响不显著。     

Strains of scattering of near-field of a point source

Three dimensional scattering of near-field is studied for dilatation and rotation in the time domain. The perturbation method is applied to solve the equation of motion for the first order scattering from a weak inhomogeneity in an otherwise homogeneous medium. The inhomogeneity is assumed close enough to the point source so that the near-field intermediate wave is dominating over the far-field sphericalP andS pulses. The integral expressions are derived to relate dilatation and rotation of scattering to the radial fluctuations of velocities and density in the inhomogeneity. These integrals are solved to calculate the strains of scattering from (a part of) an inhomogeneous spherical shell of arbitrary curvature. Variable curvature may allow the shape of inhomogeneity volume element to change uniformly from spherical to rectangular. Rotation of scattering from a spherical shell is independent ofP wave velocity inhomogeneity. Dilatation of scattering does not involveS wave velocity inhomogeneity but its gradient. The back scattering results are obtained as a special case. Strains are computed numerically, for hypothetical models to study the effects of various parameters viz., velocity inhomogeneity, distance of source from inhomogeneity and from receiver, and thickness of inhomogeneity. The curvature of the spherical shell is varied to study the effects of the shape of inhomogeneous volume element on scattering.     

土石坝心墙水力劈裂分析

简要分析了水力劈裂的发生条件、力学机理。基于水压楔劈机理,利用有限元方法建立了一种水力劈裂的发生判定方法。该方法假定心墙预先存在局部渗透弱面（裂缝）,通过将裂缝位置的单元材料改为裂缝软材料,考虑库水进入裂缝后对裂缝周围土体的作用,建立水力劈裂分析的平面有限元模型,确定裂缝端部垂直于裂缝面的正应力,进而依据该正应力判断水力劈裂发生的可能性,该方法同时可模拟水力劈裂的发展过程。