正逆累积法和分布函数法在研制储层电性标准中的应用

针对储层电性标准研制时采用的交汇图法--单一的 "L"型划分方法存在的不直观性和人为因素的影响,将正逆累积法和分布函数法首次应用到确定电性标准图版中,并以喇嘛甸油田气顶区为例,进行了气砂(含气砂岩)电性标准的研制.误差统计表明,两种划分方法均达到了国家储量规范的要求.这使得电性标准划分更加合理、精确与直观,为储层电性标准的研制提供了新的方法与思路.     

氦在矿物中的扩散及其在地震研究中的应用前景

研究氦(He)在地表和地球深部矿物中的扩散特征与封闭温度,对了解地震前兆信息和构造活动有重要意义,对了解成矿过程及规律也有重要价值。本文综述了He在石英等典型矿物中的扩散特征研究成果,主要包括：(1) He在锆石中的扩散;(2) He在磷灰石中的扩散;(3) He在石英和柯石英中的扩散;(4) He在碳酸盐中的扩散;(5) He在角闪石等含水矿物中的扩散,分析了He从矿物中扩散逃逸行为与地震活动过程中He浓度异常的关系,对未来的研究方向提出了展望。     

各向异性扩散在电磁测深曲线平滑中的应用

电磁类地球物理方法由于极易受到各类噪声的干扰,使得估计的视电阻率曲线或相位曲线发生畸变,严重地影响了反演解译精度,如何对这类观测曲线进行合理有效地平滑,是目前数据处理中的重点,也是难点.本文把各向异性扩散（anisotropic diffusion AD）引入到曲线平滑中,提出了各向异性扩散的曲线平滑方法,以平滑点梯度值的降函数作为扩散速度,在梯度较大的位置予以较小的扩散速度以保护特征,在梯度较小的位置予以较大的扩散速度平滑噪声引起的扰动.同时为了减小“飞点”(outliers)对平滑的干扰,采用局部单调扩散进行预处理.实验结果显示,本文提出的平滑方法有效地平滑了噪声与“飞点”的干扰,恢复了曲线的基本形态,有效地保持了曲线基本特征.     

Kirchhoff偏移法在探地雷达正演图像处理中的应用

本文首先通过对波动方程的分析,得出了声波波动方程和雷达波波动方程形式一致性,从而说明了把广泛应用在地震数据处理中的偏移技术引入到GPR资料处理中的可行性;其后说明了时域有限差分法（FDTD）的原理,并用它合成了几种常见的雷达正演剖面;最后利用Kirchhoff积分偏移法对正演所得的雷达剖面进行偏移处理,通过对比偏移处理前后的雷达正演剖面,可知Kirchhoff积分偏移法能使雷达正演剖面中的反射波的归位,绕射波收敛,从而大大地提高了雷达正演剖面的分辨率,更好地指导GPR剖面的地质解释和验证偏移方法的有效性.     

多模态影像技术在临床中的应用进展

自人类进入20世纪后,知识结构呈现信息爆炸征象,各项技术迅速发展。PET、PET/CT、PET/MRI、SPECT、SPECT/CT等仪器、技术得以广泛应用,并在疾病诊治等方面取得重大进展。因疾病发生机制较为复杂,各种模态信息量相对单一及片面,存在其自身局限性,使得即使拥有如此多模态的影像技术,仍无法完全实现对复杂疾病的精确评估。故此不同模态图像的信息融合使得信息互补及交叉验证成为可能。该文主要介绍多模态影像技术的发展及应用进展。     

边界元算法在电法勘探正演中的应用综述

边界元算法是一种计算量小,方便有效的处理开域电、磁场的数值模拟方法,被广泛应用于电法勘探正演领域.本文详述了用边界元算法对电法勘探中电、磁场进行的正演数值模拟,包括稳定电流场的直流电剖面法、直流电测深法和充电法,大地电磁场和人工源谐变电磁场.其中涉及了2D、3D地电体的正演计算和地形校正.对于目前应用广泛的人工源电磁法的3D边界元数值模拟,提出了六点具体的研究发展方向,希望对今后研究有一定帮助.     

正演模拟技术在地震采集设计中的应用

随着地震勘探开发的不断深入和发展,地震勘探的主战场逐渐向复杂地区转移.复杂地区既指地表条件复杂的地区,也指地下地质构造和地层条件复杂的地区,这些都对地震勘探提出了新的挑战和更高的要求.地震正演模拟正是开展此类问题研究的一个重要手段和方法.目前市场上具有正演模拟功能的软件大多是根据射线理论采用射线追踪的方法来完成正演模拟的,这种方法不能很好地反映地震波的动力学特征,特别是在复杂地区难以得到正确的结果.本文利用高阶有限差分有效克服了常规有限差分算法求解波动方程的频散问题,并以高效的OpenMP并行计算模式进行了并行优化,较大程度上提高了正演计算的速度和精度;同时实现了二阶Higdon边界条件,改善了边界吸收效果;也在一定程度上提高了计算速度.塔中地区主要目的层埋藏深,逆断层发育,地震反射特征复杂,增加了地震勘探的难度.本文依据该地区的地质模型,利用波动方程正演技术论证了该地区的地震采集观测系统,为该区地震采集观测系统的设计提供了科学依据.     

时域有限差分法在地质雷达二维正演模拟中的应用

本文从地质雷达正演原理着手,分析了差分格式中半空间步长与半时间步长的实现方法,同时通过分析数值频散的产生,进而推导出了理想频散关系和超吸收边界条件．最后文中还对比了有无边界条件的雷达正演模拟效果和精度．     

快速汉克尔变换及其在正演计算中的应用

在地球物理电(磁)法正演计算中常涉及贝塞尔函数积分,快速汉克尔数字滤波是最常用的计算方法之一.然而不同的汉克尔滤波系数,计算结果的精度也有所差异.本文对比了五组高精度快速汉克尔长滤波系数的计算精度,分析了不同核函数对计算精度的影响;同时将五种滤波系数应用于电阻率测深法、频域电磁测深法、瞬变电磁测深法的正演计算中,比较了不同滤波系数对响应和视电阻率模拟时精度的变化关系.经过大量的模型计算表明:滤波的精度并不单纯的只与滤波系数个数有关,而是多种因素综合作用的结果.本文的研究对于电(磁)法正演计算汉克尔滤波系数的选择具有指导作用,有利于节省计算机资源,提高计算效率.     

地震正演模拟在高分辨率隐伏断层地震勘探中的应用

应用波动方程有限差分方法人工合成地震记录进行城市隐伏断层的地震勘探.采用适当的差分算法、震源子波以及有效地边界条件和频散消除方法来提高正演模拟的精度和分辨率.结果表明地震数值模拟可以用来分析地层厚度、断层倾角、断层深度以及断裂带宽度等参数对地震记录的影响;通过实际试验数据和正演模拟的合成单炮地震记录的对比,可以快速判断断层的大致位置和断裂带的大致范围以及目标地层的大致深度.该方法能有效辅助实际野外工作中的数据采集参数的估计,提高勘探效率和精度.     

地震边坡动态临界加速度计算中抗剪强度参数概率分布分析

基于动态临界加速度计算地震边坡永久位移的方法应用了边坡分布的不均匀性,不均匀分布程度影响着永久位移计算结果。研究Newmark滑块底部或潜在滑动面上抗剪强度参数的分布特征,有利于了解地震过程中临界加速度的分布和变化过程,为地震边坡永久位移计算参数选取提供依据。通过计算发现：边坡土体标准差的线性关系比较微弱,在缺乏实验数据时,仅可参考使用;黏聚力随机数个数达到200时所计算的永久位移数值不再有波动;黏聚力无论设置为正态分布还是对数正态分布均不影响永久位移计算结果;黏聚力标准差的数值大小对永久位移的大小和离散性影响很大。在地震边坡计算时,尽量将其设置为符合正态分布;尽量根据实测数据设定黏聚力标准差,若缺乏实测资料,可参考两者的线性关系并结合经验进行标准差设定。     

基于离散灾情信息的地震烈度分布快速判定方法研究*

根据历次破坏性地震实际灾情反馈特点的分析总结,结合地震灾害的演变机理和分布特点,本文尝试利用3种基于离散灾情信息点推漫灾情分布面的方法在不同灾情获取阶段快速判定地震烈度分布。从计算原理、数据基础、方法特点和判定结果看,三种方法具有不同的应用场景。泰森多边形判定法简单、易行,可用于震后初期灾情量较少时使用,判定结果较粗糙,必须再结合震源特性、破裂过程以及地形地貌条件等对震害分布的影响进行综合判定后才能得出相对合理的震害分布预判结果;克里金插值判定法是一种基于地理统计方法的空间插值方法,对数据分布有明确要求,其判定结果对灾区范围整体把握较好,可用于已经获取一定量灾情信息但灾区情况尚未明朗时使用,其插值结果再结合地震破裂过程以及震后余震分布信息综合调整后,可作为灾害评估修正和应急决策参考;标准差距离椭圆拟合方法是基于地震烈度分布呈椭圆形的假设提出的,对数据的敏感度较高,可用于对可信度较高的或由专业途径上报的灾情信息进行拟合,其结果可作为现场灾害调查和确定最终烈度分布的参考。     

长波辐射(OLR)信息预测华北地震的方法研究

利用卫星对地遥测长波辐射(OLR)信息,分析了华北地块OLR信息场的空间分布及其变化特征,结合地震活动分布特点进行了短期预测中强地震判定指标的方法研究。结果表明该方法可对华北地区短期内有无可能发生M5以上的显著地震做出较为明确的诊断。经近一年的实际检验效果良好,对华北地区震情监测预报有实用价值。     

基于片法对地震低活动水平区域的构造分布特征研究

文中应用片法研究了首都圈地区在39.8°～40.4°N和114.0°～115.0°E范围内的地震活动分布。 片法主要是针对地震活动分布比较稀疏、断层无法直观识别的地区,用空间扫描的方法,并结合相应的统计检验方法,来评价地震震中分布的线性特征是否显著,从而进一步识别地震断层。 该文选定的研究区域,在过去30多年的时间里共发生了1.0～5.0级地震906次,地震的活动不活跃比较符合片法的应用研究。 将该方法用于所选定区域地震数据的研究表明,其识别能力还是较好的,可以整体上刻画比较真实的地质构造特征。 对由于种种原因无法进行实测,且地震活动较低的构造区域,片法是分析其构造分布的有效补充。     

总体样本风速风向联合概率分析方法

提出了针对总体样本的风速风向联合分布的实用统计分析方法。假定任意给定风向上的风速符合某一概型分布、风向间的关联以风速资料的相对发生频率反映、不同风向间的概型参数及频率分布的变化满足谐波函数,在此基础上,利用风速风向资料的总体样本可以得到一个能够反映任意方向的风速风向联合分布函数。利用澳门地区的百年风资料与某风站连续5年的风速观测记录,对该方法的可行性进行了验证。分析表明,该方法具有一定的适用性,便于结构的抗风可靠性分析,具有一定的工程应用价值。     

松辽盆地三肇凹陷延长排列法求取地壳反射纵波平均速度

为研究Moho几何特征，提出采用常规的延长排列法求取地壳内反射纵波平均速度，进而计算盆地基底平均波速。经过综合误差分析，得出松辽盆地内地壳平均纵波速率约为6.197km/s（相应的0为10．676s)，基底平均纵波速度约为6.75km/s。通过与全球地台区，陆－陆碰撞带，裂谷带盆地之下地壳平均波速数值的比较，认为在松辽盆地内存在上地幔物质上涌过程，上涌物质的一部分进入中部地壳，引起基底平均波速增高。讨论用延长排列法求取地壳平均波速的条件，包括Moho震相的条件，地表条件和接收排列长度的条件等。     

应用CT成像技术研究岩心孔隙度分布特征

清楚地识别岩心内部孔隙度分布特征,并对其进行定量描述,有助于深入研究岩心的渗流特性。应用CT成像技术,通过对干岩心和饱和地层水岩心扫描,重建了岩心不同扫描层面、沿某正交切面以及三维孔隙度分布图。应用孔隙度频率分布和累积频率分布曲线得到了岩心孔隙度分布的统计特征。将岩心孔隙度划分为“层面级”和“岩心级”两个级别,选用7种参数对岩心孔隙度分布的集中趋势、离散程度和分布形态进行表征。研究实例表明,同一岩心不同扫描层面孔隙度分布存在明显差异,其集中程度、离散程度和分布形态均有所不同。采用这7种表征参数可精细刻画岩心孔隙度分布特征。     

基于网络的地震灾情信息时-空分布规律及启示

以美国地质调查局基于网络的社区烈度调查系统（＂Did You Feel It？＂）1998—2008年给出的美国1510次有感地震灾情信息数据为数据源,经统计分析,初步得出地震灾情信息的时-空分布规律。我们认为利用此规律,可为地震应急灾情快速识别与评估提供一定的理论支持。     

基于GIS的震后修复建筑物空间分布格网化方法研究

针对传统格网化方法存在计算误差较大、计算效率较低等问题,提出一种新的方法——基于GIS的震后修复建筑物空间分布格网化方法。利用DEM数据、土地数据以及基础地理数据等,提取地形、地貌等各种影响因子。利用建筑面积总和以及土地利用分类数据进行多元线性回归建模,获取50 mm格网尺度下的建筑分布,利用POI数据与建筑空间的相关性,引入四叉树获取最高中心点集,在上述基础上构建格网震后修复建筑物空间分布模型,完成震后修复建筑物空间分布。实验结果表明,所提方法可有效降低计算误差,提高计算效率。     

基于Bootstrap方法确定岩土参数的标准值

提出一种基于Bootstrap重抽样技术的岩土参数标准值确定方法,将置信概率为95%对应的重抽样均值作为标准值,并通过一组实测土体的抗剪强度参数阐述标准值的确定方法和过程。对两个实例进行研究,分别探讨Bootstrap重抽样次数(即Bootstrap样本集大小)和原始数据的变异性对标准值精度和收敛性的影响。结果表明,基于Bootstrap方法确定的岩土参数标准值具有较好的收敛性和精度,当重抽样次数不小于100时,Bootstrap重抽样次数和原始样本的变异性对标准值的收敛性和精度不会产生显著影响。该方法能够有效地处理岩土参数的不确定性问题。