一阶弹性波方程交错网格高阶差分解法   总被引：150，自引：34，他引：116       下载免费PDF全文

董良国  马在田  曹景忠  王华忠  耿建华  雷兵  许世勇 《地球物理学报》2000,43(3):411-419

提高计算精度和运算效率是所有波场正演方法所追求的目标,本文通过将速度 （应力）对时间的奇数阶高阶寻数转化为应力（速度）对空间的导数,运用时间和空间差分精度 均可达任意阶的高阶差分法,通过交错网格技术,对一阶速度－应力弹性波方程进行了数值求 解．波场快照以及实际模型的正演结果表明,这种求解一阶弹性波方程的高阶差分解法,和 常规的差分法相比网格频散显著减小,精度明显提高,而且可以取较大的空间步长,提高计算 效率。  相似文献   

风化花岗岩声谱特征分析   总被引：7，自引：0，他引：7  

刘彤  徐鸣洁  胡德昭  耿建华 《高校地质学报》2000,6(4):588-594

岩石材料可视为声波的天然低通滤波器 ,声波在不同力学性质岩石材料中传播时 ,由于岩石对声波吸收程度不等因而高频声波滤除状况也不同。接收波中频率丰富程度与岩 石 的力学性质密切相关,岩石力学强度愈高,声谱中高频成分就愈丰富。反之,声谱中主要为 低频,高频分量因衰减而缺失。实验采用全数字化DB4型多波参数分析仪和计算机,对波形 数据进行FFT分析,通过接收谱和频率响应的特征来研究岩石的力学性质。声衰减研究可为 岩石、岩体强度和综合评价提供更多、更可靠的信息,是解决工程地质风化岩问题研究的一 种可行方法。  相似文献   

新元古代帽碳酸盐岩中帐篷状构造的成因     

刘芊  陈多福  冯东 《地学前缘》2007,14(2):242-248

帐篷构造是碳酸盐岩中的一种特殊沉积组构,因其倒“V”字形形态类似于帐篷(tepee)而得名,被认为是一种无成因意义的沉积组构。传统的帐篷构造在海相至陆相碳酸盐岩中都有发育,形态和胶结物的不同反映了沉积环境的变化,其成因为裂隙填充的胶结物结晶膨胀导致层面突起变形。新元古代帽碳酸盐岩中广泛发育有倒“V”字形的类似构造,但由于形态和成因上都和传统的帐篷构造有所区别,被称为“帐篷状构造”(te-pee-likestructure)。目前其成因解释主要有:“巨风暴潮波痕”、“甲烷气体渗漏”、“地下水侵位”和“晶体结晶”的假说。由于帐篷状构造的形成过程与机制和帽碳酸盐岩的成因密切相关,对帐篷状构造的进一步研究必将帮助我们对新元古代冰期结束机制的理解。  相似文献   

海底节点长期地震观测:油气田开发与CO2地质封存过程监测     

耿建华  董良国  马在田 《地球科学进展》2011,26(6):669-677

地震观测方法目前被广泛用于海底油气田勘探开发与CO2地质封存过程监测.比较了海面拖缆地震观测与海底节点地震观测方法的差异.海底节点地震观测方法能够采集多分量、全方位、高信噪比和高分辨率数据,同时具有数据定位精度高、采集可重复性以及数据采集性价比高等优势.在回顾海底节点地震观测方法发展过程与现状的基础上,展望了海底长期节...  相似文献   

低信噪比海底多分量地震数据波场分离方法          下载免费PDF全文

刘学义  程玖兵  王腾飞  耿建华  刘玉柱 《地球物理学报》2021,64(2):684-699

以多分量地震观测为基础,联合纵波和转换横波数据能更有效地估计地下介质的弹性和物性参数,提升地质构造成像与油气储层描述的精度.在海底多分量地震数据处理过程中,观测记录的上-下行波分解和P/S波分离可压制水层鸣震以及P与S波之间的串扰,对偏移成像和纵横波速度建模至关重要.但受海底环境、仪器与观测因素共同影响,许多海底多分量地震资料都无法基于现有的海底波场分离方法与流程取得合理的结果.本文以海底声波场与弹性波场分离基本原理为基础,通过对方法流程的修正,摆脱常规流程对中小偏移距直达波信号的依赖性.借助模拟数据实验讨论了波场分离对海底介质参数、噪声的敏感性.结合东海YQ探区海底多分量地震资料上-下行P/S波分离及其叠前深度偏移处理,验证了本文方法流程的可行性.  相似文献   

泥灰岩中自生方解石的稀土元素酸溶方法研究   总被引：1，自引：0，他引：1  

陈琳莹  李崇瑛  陈多福 《地球化学》2014,(6):647-654

碳酸盐岩稀土元素是沉积环境示踪的理想参数,但由于自生沉积碳酸盐矿物的稀土元素含量远低于陆源黏土等其他矿物,在样品处理过程中少量陆源黏土矿物的溶解,将影响沉积环境的示踪。因此,如何在样品处理过程中避免陆源黏土矿物中稀土元素的释放,从而获得海洋环境自生沉淀形成的碳酸盐矿物中的元素方法非常重要。本文应用地层中常见的同一泥灰岩样品,分别用5%醋酸、5%盐酸、2%和5%硝酸浸取,采用0.5 h、1 h、2 h、3 h、6 h和24 h的反应时间,用X射线衍射法（XRD）分析酸不溶残渣矿物组成、电感耦合等离子体原子发射光谱法（ICP-AES）测定酸可溶相溶液的主元素和微量元素含量、电感耦合等离子体质谱法（ICP-MS）测定酸可溶相溶液的稀土元素含量,结果表明所有酸溶均可完全溶解泥灰岩中的方解石矿物,并有少量黏土发生溶解,导致稀土总量和Ce/Ce＊增高,但对Eu/Eu＊、LaN/YbN和Gd/Gd＊等参数没有影响。泥灰岩中自生方解石的稀土元素理想溶样条件是5%醋酸反应1 h以内,溶液的ICP-MS分析结果适用于示踪沉积环境的氧化还原条件。  相似文献   

基于离散粒子理论地震波传播数值模拟网格剖分计算方法     

高伟  耿建华 《地球物理学进展》2013,28(3)

介绍一种基于离散粒子理论地震波传播数值模拟的网格剖分计算方法.根据离散粒子理论,将研究区域划分为由一系列相互作用的粒子组成的正六边形网格,这些粒子在它们的接触点处发生相互作用,并用Hooke定律和Newton定律描述.为解决六边形网格带来的网格交错而难以计算以及波场输出问题,将横向网格进行加密,加密处赋予假想的粒子,输出波场时选取偶数行偶数列点或奇数行奇数列点的波场值.均匀介质和层状介质模型的数值模拟结果表明,该网格剖分计算方法能够将离散粒子理论用于模拟弹性波在非均匀各向同性介质中地震波的传播.  相似文献   

青藏高原冻土带天然气水合物的形成条件与分布预测   总被引：44，自引：7，他引：37       下载免费PDF全文

陈多福  王茂春  夏斌 《地球物理学报》2005,48(1):165-172

冻土带是天然气水合物发育的两个重要地质环境之一.青藏高原平均海拔在4000m以上,多年冻土面积约1.4×106km2.本文根据青藏高原冻土层厚度和地温梯度特征,运用天然气水合物的热力学稳定域预测方法,确定中低纬度高海拔区冻土带天然气水合物的产出特征.青藏高原多年冻土带热成因天然气水合物形成的热力学相平衡反映,水合物顶界埋深约27~560m,底界埋深约77~2070m.初步计算表明,青藏高原冻土带水合物天然气资源约1.2×1011~2.4×1014m3.在冻土层越厚、冻土层及冻土层之下沉积层的地温梯度越小的地区,最有利于天然气水合物的发育.气温的季节性变化对天然气水合物影响不大.在全球气温快速上升的背景下,青藏高原天然气水合物将处于失稳状态,天然气水合物顶界下降、底界上升,与冻土带的退化相似,分布区逐渐缩小,最终将完全消失.  相似文献   

地质历史时期天然气水合物分解释放的地质地球化学证据     

王钦贤  陈多福 《现代地质》2010,24(3):552-559

巨量天然气水合物广泛赋存于大陆坡海底和高纬度永久冻土带沉积物中,是重要的环境影响因素。有些地质历史时期发生的重大地质灾变可能是由天然气水合物分解释放甲烷引起的。持续时间短暂且显著的全球性海相和陆相碳同位素的负漂移、甲烷起源的极负碳酸盐岩碳同位素值(特别是δ¹³C<-40‰)、碳酸盐岩溶解、指示甲烷浓度异常的生物标志化合物、生物成因重晶石富集、类似冷泉碳酸盐岩发育的特殊沉积组构、海底垮塌及碳酸盐结晶扇(如文石)等是支持天然气水合物分解“假说”的主要证据。但这些证据还存在一定的局限性,该“假说”仍然存在较大的争论,因此需要进一步深入研究。  相似文献   

天然气水合物似海底反射层(BSR)AVA特征:双相介质模型     

麻纪强  耿建华 《应用地球物理》2008,5(1):57-66

The bottom simulating reflector （BSR） in gas hydrate-bearing sediments is a physical interface which is composed of solid, gas, and liquid and is influenced by temperature and pressure. Deep sea floor sediment is a porous, unconsolidated, fluid saturated media. Therefore, the reflection and transmission coefficients computed by the Zoeppritz equation based on elastic media do not match reality. In this paper, a two-phase media model is applied to study the reflection and transmission at the bottom simulating reflector in order to find an accurate wave propagation energy distribution and the relationship between reflection and transmission and fluid saturation on the BSR. The numerical experiments show that the type I compressional （fast） and shear waves are not sensitive to frequency variation and the velocities change slowly over the whole frequency range. However, type II compressional （slow） waves are more sensitive to frequency variation and the velocities change over a large range. We find that reflection and transmission coefficients change with the amount of hydrate and free gas. Frequency, pore fluid saturation, and incident angle have different impacts on the reflection and transmission coefficients. We can use these characteristics to estimate gas hydrate saturation or detect lithological variations in the gas hydrate-bearing sediments.  相似文献