排土场土石混合体分形特征与抗剪强度特性研究

露天矿山剥离堆积的土石混合体颗粒尺度不一,且力学性质异常复杂。以西藏甲玛铜矿南坑排土场为研究对象,采用现场勘测的方法测定堆积体的粒度组成及其分布特征,开展室内重塑土样的固结排水三轴压缩试验,以分析土体屈服强度特性与剪胀性,以及抗剪强度参数的取值。研究结果表明,自排土台阶坡顶至坡底,土石混合体的土粒减少且块石增多,颗粒的平均粒径呈指数关系增大;土体的粒度分布具有显著的分形结构特征,随块石含量的增大,其粒度分形维数减小。块石含量不同的土样在偏差应力作用下均表征出压密、屈服、塑性流动破坏的剪切特性,且围压条件与块石含量是影响土体剪切变形的两个关键因素,低围压时增加块石含量,土体更易产生剪胀。摩尔库伦强度准则能客观反映土体屈服的现象及过程,随块石含量的增加,土体的黏聚力呈线性减小且内摩擦角呈指数增大。     

沉积体系非均质性对致密油气储集层控制及其定量评价方法

以鄂尔多斯盆地延长组长8-长6致密储集层为例,采用地质统计数学(Levy Sim)和沉积地层正演模拟对沉积非均质性进行研究。结果表明,地质统计学Levy函数能够有效地表征非均质性参数,不仅可以表征油藏尺度(岩心、露头和测井曲线尺度),而且也可以有效刻画岩石微观尺度的非均质性。而三维正演地层模拟方法 Sedsim则能够有效地模拟浅水三角洲砂体的空间展布。沉积地层正演模拟与非传统地质统计学相结合,能够更真实地定量刻画不同尺度的沉积非均质性,该方法不仅有效地解决了不同尺度沉积储集层、烃源岩及盖层非均质性定量表征的难题,而且为非常规油气成藏研究提供了新思路与新的定量评价方法。     

不排水抗剪强度非平稳随机场模拟及边坡可靠度分析

边坡可靠度分析中通常假定采用平稳或准平稳随机场表征土体参数的空间变异性,然而大量现场试验数据表明,土体参数如不排水抗剪强度沿土体埋深常呈现明显的非平稳分布特征,即其均值和标准差均随埋深发生变化,因此亟需发展土体参数非平稳随机场模型及其模拟方法。针对目前不能有效单独模拟土体参数趋势分量和随机波动分量的不确定性,提出了一种有效的不排水抗剪强度参数非平稳随机场模型,并给出了土体参数二维非平稳随机场模拟方法计算流程,同时将新提出的模型与现有非平稳随机场模型及平稳随机场模型进行了系统比较。最后通过不排水饱和黏土边坡算例验证了提出模型的有效性,并揭示了不排水抗剪强度非平稳分布特征对边坡可靠度的影响规律。结果表明：提出模型能够有效地单独模拟土体参数趋势分量和随机波动分量的不确定性,考虑土体参数均值和标准差随埋深增加而增大的特性,可为表征土体参数非平稳分布特征提供了一条有效的途径。此外,与采用非平稳随机场模拟土体参数空间变异性相比,采用常用的平稳随机场模型会低估边坡失效概率,从而造成偏危险的边坡工程设计方案。     

基于数字图像分析的冰水堆积体结构建模与力学参数研究

冰水堆积体是一种典型的混合介质,由软弱“土”与坚硬“石”构成,其参数选取和分析方法困难。利用数字图像重构方法模拟冰水堆积体结构特征,可反映块石分布特征;以基于图像识别的块石分割和颗粒离散元建模分析方法,研究了堆积体内介质分布、土体颗粒强度对混合介质力学参数的影响规律,分析了冰水堆积体变形、峰值抗剪强度及峰后段抗剪强度的变化规律。结果表明：冰水堆积体混合介质在峰值前的变形曲线呈非线性硬化性,峰后呈现塑性特性,其综合抗剪摩擦角随基质土的黏结强度提高而降低;当块石含量小于25%时,抗剪切强度近似等于基质土强度,块石含量高于25%时,抗剪切强度随块石含量增加而提高;由于受块石空间分布影响,块石含量与摩擦角间的相关性要比其与黏结强度间的高;该成果应用于云南古水水电站坝前冰水堆积体发现,含块石混合介质内摩擦角比纯土高0°～8°,与工程判断较吻合。该方法能弥补室内试验不能考虑大尺寸块石的困难,可作为室内试验获取力学参数方法的有益补充。     

随机分形在刻划储层非均质特性中的应用

储集层的非均质性质普遍存在于砂泥岩储层、碳酸盐岩和火成岩储层中。线性数学手段难以精确描述非均质储层参数的空间分布特征。在详细分析砂泥岩地层参数非均质特性的基础上,将分形方法与克里金方法相结合来刻划储层密度和声波时差的分布规律,通过对两口井的相应参数进行插值产生的虚测井值与真实测井值吻合良好。在井距合理的情况下,统计分形方法能够较好地刻划储集层的非均质性。     

土石混合体渗透侵蚀特性数值模拟研究

土石混合体是由块石、砾石与砂土、黏土等固体颗粒组成,以块石镶嵌在细粒土体中为结构特征的松散介质。由于其颗粒松散特性,细颗粒在渗流作用下易于流失而形成渗流侵蚀。采用室内试验与数值模拟相结合的方法,扩展渗流侵蚀模型在土石混合体渗透侵蚀研究中的应用。选取云南某滑坡土石混合体为研究对象,采用自制的渗透侵蚀装置,模拟土体基质渗透侵蚀过程,获取土体基质侵蚀参数试验值。根据块石在土体基质中分布特征,建立不同块石含量特征的土石混合体计算模型。基于参数试验值及渗流侵蚀方程、耦合渗流控制方程及孔隙率控制方程,模拟土石混合体渗流侵蚀过程,分析块石含量对土石混合体渗透侵蚀特性的影响。模拟结果表明:在同等水力梯度及时间条件下,块石含量与土石混合体宏观平均流速初值、侵蚀颗粒密度值呈反比,与液化颗粒密度峰值呈正比;块石含量与渗流侵蚀方程中试验参数存在线性相关性。将块石含量作为土石混合体结构特征引入渗流侵蚀方程中,以构建适用于土石混合体的渗流侵蚀模型。     

储层非均质性研究方法进展

储层非均质性作为储层表征的核心内容之一,一直是研究者关注的焦点。储层非均质性的研究内容主要包括储层非均质性分类、储层非均质性评价、储层非均质性成因分析和储层非均质性对油气田开发的影响等。国外对于储层非均质性研究主要包括利用各种实验研究储层孔渗等物性特征及其对流体渗流规律的影响、利用野外露头资料开展地质建模和模拟研究进行储层非均质性定量分析。国内主要包括通过沉积学、成岩作用等地质成因分析研究储层非均质性和基于测井精细解释的渗透率等储层物性参数进行非均质性的定量评价等。研究方法有地质成因分析法、层序地层学方法、各种实验研究方法、地质统计学方法、各种数学计算方法、测井解释方法、地质建模方法、地震预测方法、数值模拟方法和生产动态分析法等。结合自身科研实践,分析了不同研究方法的优缺点。在此基础上指出了目前储层非均质性研究中存在的主要问题和今后的发展方向。     

二维随机介质模型正演指导复杂地层地震解释——以薄层砂泥岩交互地层为例

大港油田砂泥岩交互地层存在着非均质复杂地质形态,为保证其精细解释效果,根据随机过程理论,提出了选用二维指数型椭圆自相关函数为复杂非均质地质体进行简化建模的方法。正演算法使用交错网格有限差分法模拟弹性波在二维随机介质模型中的传播,激发方式使用自激自收模式。正演结果表明,随机介质模型模拟大港油田薄层砂岩互层地层具有合理性,交错网格有限差分算法对于该模型的正演数值模拟,算法合适。     

基于MODFLOW参数不确定性的地下水水流数值模拟方法

考虑到模型不确定性引起的地下水数值模拟不确定性对模拟过程的影响,在简要介绍含水层水文地质参数变异性研究进展和地质统计学的基础上,基于常用的确定性地下水流数值模拟软件MODFLOW开发了MODFLOW-Gslib软件,相较于传统的数值模拟方法,将地质统计学与数值模拟结合的方法能够模拟非均质含水层中的参数变异性问题。将MODFLOW-Gslib软件运用于模拟实例中,选择常见的不确定性因素进行模拟,并对其模拟产生的数据进行统计分析,结果表明,软件转化后的参数符合水文地质参数不确定性的相关特征;与原模拟结果进行对比,该软件能够更加真实地刻画含水层参数变异性特征。     

下寺湾油田长8段储层相控建模研究

运用Petrel储层建模软件,对下寺湾油田长8段砂组进行储层相控建模研究。首先利用序贯指示模拟方法建立能够表征储层较大规模非均质性的沉积微相模型;在沉积微相模型建立之后,采取相控物性参数建模方法,分小层分微相统计各种物性参数分布特征;然后利用序贯高斯模拟方法分别进行模拟,建立了研究区的孔隙度模型、渗透率模型和含油饱和度模型。结果证明按此方法建立的模型效果较好,客观反映了地质实际,为该区下一步滚动勘探开发提供了地质依据。     

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

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

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.     

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

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

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.     

土体再压缩变形规律的试验研究

结合大量不同土性土体的回弹再压缩试验、模型试验,提出再加荷比、再压缩比率的概念,并在此基础上得出了土体再压缩变形的基本规律：当再加荷量为卸荷量的20%时,土样产生的再压缩变形量已接近回弹变形量的40%~50%;当再加荷量为卸荷量的80%时,再压缩变形量与回弹变形量大致相等,则此时回弹变形完全被压缩;当再加荷量与卸荷量相等时,再压缩变形量大于回弹变形量,且再压缩变形的增大程度与土性有关。可见在土样的再压缩过程中,在初始阶段再压缩变形增长速率较大,之后增长速率随着加荷量的增加反而逐渐降低。由此得出土体的再压缩变形发展规律为两阶段线性规律,这一规律具有工程实用意义：基底以下土体的回弹再压缩变形对于减小主群楼之间的差异沉降是一种有利因素;再压缩变形的发展规律为建筑物基底以下土体的回弹再压缩变形而产生的沉降计算提供了依据     

土石坝心墙水力劈裂分析

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

土钉支护失稳分析

某基坑采用复合土钉墙支护,开挖至坑底后,基坑北侧发生失稳事故。本文就其事故发生原因作一探讨。     

地震次生地质灾害风险评估集对态势分析

5．12四川汶川8级强震引发了大量次生地质灾害,加剧了灾情并严重影响抗震救灾。地震次生地质灾害受多重因素制约,其孕育与致灾过程具有不确定性。利用集对分析原理与方法,对地震次生地质灾害进行同异反态势分析,拟定了风险评估的集对分析同一度、差异度、对立度等指标体系的构建原则与赋值标准,对承灾体系统的不确定性及其作用作了刻画与分析,建立了不同风险分区代表性的集对分析联系度表达式,为地震灾害研究提供了可资借鉴的新思路和技术方法。     

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

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