鄂尔多斯盆地黄陵、东胜地区地温场对比

鄂尔多斯盆地黄陵、东胜铀矿区分别处于盆地南部渭北隆起的北侧边缘和盆地北部伊盟隆起的东部,赋矿层位都是中侏罗统直罗组。盆地南、北铀矿区在现今地温场及古地温场都存在明显差异,南部现今大地热值及热演化程度明显高于北部。对于下侏罗统延安组和石炭—二叠系煤层,黄陵地区镜质体反射率都高于东胜地区。通过镜质体反射率资料得出同一埋深的一套地层经历的最大古地温和对应的古地温梯度也有南部高于北部的现象。由于早白垩世后期盆地普遍整体抬升使得现今地温相对古地温降低,南部黄陵地区抬升剥蚀量大于北部东胜地区,导致古、今地温差异也大于后者。盆地南部庆阳—富县一带局部构造热运动,导致南部异常地温场的形成,使得南部热演化程度高于北部。     

西藏扎布耶盐田盐渍土的改性配比试验

根据盐渍土的物理化学性质,采用掺入不同比例石灰或水泥的方法进行改性处理。经击实、抗剪和抗压强度试验结果表明,盐渍土掺入石灰或水泥改性后,混合料的最大干密度、最优含水率、强度指标等均有增加,特别是内聚力迅速增长。研究表明,盐渍土混合料的性状得到了明显改善,改性后的混合料可用于修建盐田埝堤类相关工程。     

苏里格气田盒8气藏气水层测井识别研究

鄂尔多斯盆地苏里格气田盒8气藏为弹性驱动的河流相低效气藏,储层非均质性强,在开发中普遍见水。根据试采情况将气藏储层分为气层、水层、气水层三类,开展气、水层的测井识别研究。以区内26口井50个射孔层段各层的试油结果作为依据,选出了700个样本点作为判识建模的标本,以声波时差和中子孔隙度等8个原始测井值作为判识变量组合,采用模糊判别分析,建立了各层的判识模型。回判结果表明,模型对气层精度达到97％,水层为89％,气水层为94％。利用未参加建模的出水井射孔段对模型进行验证,判识结果基本符合实际情况,模型具有较高精度。     

邢东矿井构造特征及应力应变场分析

邢东矿井位于太行山隆起带与华北沉陷带的过渡带,通过区域构造应力-应变场分析,本区先后经历了由挤压-拉张的变形过程,变形强度由弱到强。海西期构造应力影响较弱,印支期-中燕山期构造应力作用较为强烈,形成了大量的挤压结构面和压缩构造．燕山晚期-喜山期是构造形成的最终阶段,以拉张伸展变形为主,现今的多数正断层的断裂结构面仍保存有先期的挤压结构面特征。区内构造以断裂构造为主要形式,断层组合形式多样,这些构造形成期次不同,北西向断裂对北东向断裂具有限制作用。     

松辽盆地宋芳屯油田扶杨油层孔隙结构及其影响因素

通过镜下观察,宋芳屯油田扶杨油层孔隙类型较多,以粒间孔为主。细管压力曲线上看,曲线较陡,分选差,孔隙结构特征参数总体上较差,与孔隙度、渗透率有一定的关系,大的孔隙对渗透率的贡献最大。扶杨油层的孔隙结构较复杂主要由于成岩作用影响,压实作用使孔隙空间显著缩小,胶结作用和自生粘土矿物生长使孔隙进一步缩小和复杂化,后期的溶蚀作用使孔隙空间增加。     

Reflection of waves from the boundary of a random elastic semi-infinite medium

In this paper the smooth perturbation technique is employed to investigate the problem of reflection of waves incident on the plane boundary of a semi-infinite elastic medium with randomly varying inhomogeneities. Amplitude ratios have been obtained for various types of incident and reflected waves. It has been shown that an incidentSH orSV type of wave gives rise to reflectedSH, P andSV waves, the main components beingSH andP, SV in the respective cases. The reflected amplitudes have been calculated depending upon the randomness of the medium to the square of the small quantity , where measures the deviation of the medium from homogeneity. An incidentP-type wave produces mainly aP component and also a weakSH component to the order of ². The reflected amplitudes obtainable for elastic media are also altered by terms of the same order. The direction of the reflected wave is influenced by randomness in some cases.     

地震综合效应场函数及其在地震预报中的应用──（一）地震综合效应场函数及其计算方法

在系统地分析了目前各种测震学地震预报方法科学思路的基础上,认为测震学地震预报方法基本上可以分为两大类。一类是以已经发生的一些地震作为未来可能发生的地震的“因”,即由于已经发生的地震对区域应力场的影响,导致未来发生较强地震。这一类包括的预报方法较多,如空区、条带、ｂ值、地震迁移、相关地震等等及其由此衍生出来的各种方法。另一类是把已经发生的一些地震作为区域应力场增强的“果”,即已经发生的地震是区域应力场增强过程中的一种反映,而未来地震不一定是已经发生的地震所导致的结果。这一类包括“地震窗口”、小震群活动等方法。针对第一类方法,各种预报方法都是力图从地震三要素中提取未来地震的信息,而具体作法又都是利用地震三要素这个多维空间的某个剖面。为了从地震活动诸要素的多维空间提取综合信息,我们对每个地震加入了破裂面方位,构成了地震第四要素,并依据地震４要素建立了地震综合效应场函数。地震综合效应场函数概括了多种测震学地震预报方法的科学思路和预报经验,从而可以形成测震学的综合预报方法。     

极性转换期间地球磁场形态学研究

本文对采自中国黄土高原西峰（35.7°N,107.6°E）和段家坡（34.2°N,109.2°E）两个剖面中黄土层L8和古土壤层S8的1281块定向古地磁样品做了详细的岩石磁学和古地磁学研究.证实了Matuyama-Brunhes（M-B）极性转换带位于L8的中下部.提出了下列观点:1.M-B极性转换过程与地球磁场方向变化相联系的持续时间为3600-4500a,而与地球磁场强度变化相联系的持续时间则为8000-9000a,即强度变化存在“超前和滞后” 效应;2.M-B转换场的形态是由三次快速倒转和一次不成功的倒转构成,或者说,转换场具有快速变换极性的振荡特征;3.M-B转换过程中地球磁场并不是以轴对称的非偶极子场为主,而是偶极子场至少与非偶极子场相当;4.中国黄土-古土壤沉积物所含磁性矿物的主要成分是磁铁矿,它是研究极性转换期间地球磁场详细结构的良好物质.     

Efficient generation of conditional simulations by chebyshev matrix polynomial approximations to the symmetric square root of the covariance matrix

Consider the problem of generating a realization y₁ of a Gaussian random field on a dense grid of points ₁ conditioned on field observations y₂ collected on a sparse grid of points ₂. An approach to this is to generate first an unconditional realization y over the grid =_{1 2}, and then to produce y₁ by conditioning y on the data y₂. As standard methods for generating y, such as the turning bands, spectral or Cholesky approaches can have various limitations, it has been proposed by M. W. Davis to generate realizations from a matrix polynomial approximations to the square root of the covariance matrix. In this paper we describe how to generate a direct approximation to the conditional realization y₁, on ₁ using a variant of Davis' approach based on approximation by Chebyshev polynomials. The resulting algorithm is simple to implement, numerically stable, and bounds on the approximation error are readily available. Furthermore we show that the conditional realization y₁ can be generated directly with a lower order polynomial than the unconditional realization y, and that further reductions can be achieved by exploiting a nugget effect if one is present. A pseudocode version of the algorithm is provided that can be implemented using the fast Fourier transform if the field is stationary and the grid ₁ is rectangular. Finally, numerical illustrations are given of the algorithm's performance in generating various 2-D realizations of conditional processes on large sampling grids.     

Estimation for partially observed Markov processes

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.