伪谱法地震波正演模拟的多线程并行计算

地震波的正演模拟,尤其是3D正演模拟,往往涉及大规模的数据存储和计算,问题的规模往往超出计算机的物理内存,或者计算时间让问题的求解者难以忍受,即使采用目前存储和计算能力很强的计算机,其计算费用仍然是十分昂贵的.本文提出一种基于多线程协同使用多CPU和计算域分割的正演模拟并行计算技术,使得问题的求解过程得以加快,大大地缩短了用户等待的时间.为了检验我们的并行算法的可行性,文中以傅利叶正演模拟技术为例,给出了声波和3D各向异性弹性波模拟的例子,并对不同版本（串行、并行）运算效率进行了比较,证实了方法的有效性.     

四川盆地构造沉降特征及成因机制分析

盆地的沉降过程能够反映盆地的演化历史及成盆机制。为深入分析四川盆地构造沉降特征,本文基于最新钻井资料和地震数据,通过回剥反演方法,进行去压实、沉积负载、古水深和海平面变化校正,重建了四川盆地不同构造单元的构造沉降史。同时根据瞬时均匀伸展模型和裂后热坳陷模型进行正演模拟,对盆地成因进行分析。构造沉降史的恢复揭示了四川盆地具有典型的克拉通盆地沉降特征。四川盆地的形成演化可划分为震旦纪—早古生代、石炭纪—三叠纪、侏罗纪—白垩纪3个构造沉降旋回,盆地经历了晚震旦世—早寒武世、早志留世、晚二叠世—早三叠世以及中晚侏罗世4幕快速沉降,第一幕和第三幕快速沉降期为岩石圈伸展减薄引起,另外两幕为前陆盆地发育过程中所引起的快速沉降。构造沉降正演结果表明四川盆地在寒武纪—奥陶纪和晚二叠世—三叠纪经历了两期“快速沉降—缓慢沉降”过程,快速沉降受控于岩石圈的伸展作用,缓慢沉降为岩石圈热冷却作用所主导。盆地在热冷却沉降阶段后进入前陆挠曲沉降,出现不同规模的剩余沉降。

     

基于叠前同时反演的致密砂岩储层预测及含气性识别——以苏里格S区块为例

苏里格S区块储层与围岩阻抗差异小、气水关系复杂,叠后纵波阻抗反演方法难以区分优质储层,通过研究区岩石物理分析可知,叠前参数纵横波比可以有效区分岩性和含气性。本文首先进行模型正演确定储层地震响应特征;然后通过适用于砂泥岩的Xu-White模型进行岩石物理建模及横波预测,建立岩石物理模板;其次采取线性去噪、剩余振幅补偿等处理方法对CRP道集进行优化;最后通过叠前同时反演定量预测储层厚度及含气性。结果表明：(1)储层顶界的地震反射特征为强波谷反射,底界反射不明显;(2)纵横波比小于1.68可有效划分砂岩,联合纵波阻抗小于12 200 g·cm^-3·m·s^-1预测储层,再由更低的纵横波比小于1.57识别含气性;(3)储层展布特征与含气性预测范围趋势相似,但在局部上存在差异,储层发育程度与其含气性不一定呈正相关。本文所提方法期望为下一步圈定有利含气储层面积及井位部署提供有力的技术支撑。     

老挝波乔省某磁铁矿区磁法勘探

勘查区位于老挝西北部之波乔省,地表风化严重、浮土掩盖厚、深林茂密、基岩露头少,虽然局部地段出露有大量的铁帽,但不敢判断该区是否有大规模的磁铁矿.为此,我们采用1∶10000地面高精度磁法测量以了解磁铁矿的有无、范围、延伸、规模.后期勘探阶段开展2.5D人机交互正反演工作以了解深部磁铁矿的有无.结果表明,在老挝西北部地质勘查程度低的地区开展高精度磁测扫面工作是寻找磁铁矿产行之有效的方法,能够解决基础地质问题,矿山勘查问题,后期运用2.5D人机交互正反演能够解决一定的矿层预测工作.     

声波方程频率域有限元参数反演

推导出频率域有限元声波正演方程,为了消除边界反射,将Clayton-Engquist旁轴波动方程吸收边界条件引入频率域,并对有限元刚度矩阵和质量矩阵进行压缩存储,利用广义共轭梯度法求解有限元方程获得正演解.在此基础上,推导出在某一频率下波场数据残差δU与单元物性参数修改量δλ之间关系的Jacobi矩阵,反演方法允许利用地面二维炮集全波场资料与给出初始模型参数的正演值的差值δU,迭代求得δλ.由于计算机内存的限制,方法计算不允许有过多数目的未知数个数,因此还提出了对同一介质物性单元的Jacobi矩阵元素进行压缩组装的措施,从而使反演的未知量个数减少,结合采用共轭梯度迭代法,使得只需利用有效波频段的少数一些频率即可进行迭代反演.正演和反演理论模型的数值模拟结果表明方法是有效的.     

弱各向异性介质弹性波的准各向同性近似正演模拟

准各向同性(QI)近似可用于弱各向异性介质的正演模拟.本文通过运用QI方法的零阶和一阶近似,计算了VTI介质模型的地震记录.得出的地震记录与标准各向同性射线理论(IRT)和基于伪谱法的三维地震正演模拟得出的地震记录作了比较,可以认为是精确的合成地震记录.     

Implications of the Precambrian Non-stromatolitic Carbonate Succession Making up the Third Member of Mesoproterozoic Gaoyuzhuang Formation in Yanshan Area of North China

A particular non-stromatolitic carbonate succession making up the third member of the Mesoproterozoic Gaoyuzhuang (高于庄) Formation might demonstrate that a stromatolite decline of the Mesoproterozoic occurring at ca. 1 450 Ma besides other three events of the Proterozoic,respectively,occurred at ca. 2 000 Ma,ca. 1 000 Ma,and ca. 675 Ma. The forming duration of this non-stromatolitic carbonate succession can be generally correlative to that of a similar depositional succession in North America,i.e. a non-stromatolitic carbonate succession made up by the Helena Formation of the Belt Supergroup,which suggests that the stromatolite decline occurring at ca. 1 450 Ma may be a global event. This information endows the non-stromatolitic carbonate succession making up the third member of the Gaoyuzhuang Formation in the Yanshan (燕山) area with important significance for the further understanding of Precambrian sedimentology. The Mesoproterozoic Gaoyuzhuang Formation in Yanshan area is a set of more than 1 000 m thick carbonate strata that can be divided into four members (or subformations). The first member (or the Guandi (官地) subformation) is marked by a set of stromatolitic dolomites overlying a set of transgressive sandstones; the second member (or the Sangshu'an (桑树鞍) subformation) is a set of manganese dolomites with a few stromatolites; the third member (or the Zhangjiayu (张家峪) subformation) is chiefly made up of leiolite and laminite limestones and is characterized by the development of molar-tooth structures in leiolite limestone; the fourth member (or the Huanxiusi (环秀寺) subformation) is composed of a set of dolomites of stromatolitic reefs or lithoherms. Sequence-stratigraphic divisions at two sections,i.e. the Jixian (蓟县) Section in Tianjin (天津) and the Qiangou (千沟) Section of Yanqing (延庆) County in Beijing (北京),demonstrate that a particularly non-stromatolitic succession making up the third member of the Mesoproterozoic Gaoyuzhuang Formation is developed in the Yanshan area of North China,in which lots of grotesque matground structures (wrinkle structures and palimpsest ripples) are developed in beds of leiolite limestone at the Qiangou Section and lots of molar-tooth structures are developed in beds of leiolite limestone at the Jixian Section. The time scale of the Gaoyuzhuang Formation is deduced as 200 Ma (from 1 600 Ma to 1 400 Ma). The duration of an obvious hiatus between the Gaoyuzhuang Formation and the underlying Dahongyu (大红峪) Formation is deduced as 50 Ma to 100 Ma,thus the forming duration of the GaoyuzhuangFormation is thought as 100 Ma (1 500 Ma to 1 400 Ma). Furthermore,the age of the subface of the third member of the Gaoyuzhuang Formation that is just in the mid position of the Gaoyuzhuang Formation can be deduced as about 1 450 Ma,which is the basis to infer a stromatolite decline of the Mesoproterozoic occurring at ca. 1 450 Ma. Importantly,several features of both the molar-tooth structure and the stromatolite,such as the particular forming environment,the important facies-indicative meaning,and the episodic distribution in the earth history,might express the evolutionary periodicity of the surface environment of the earth and can provide meaningful clues for the understanding of the Precambrian world,although their origin and forming mechanism is highly contentious. Therefore,like other three stromatolitic declines,respectively,occurring at ca. 675 Ma,ca. 1 000 Ma,and ca. 2 000 Ma,the identification of the stromatolite decline occurring at ca. 1 450 Ma during the Golden Age of stromatolites (2 800 Ma to 1 000 Ma) has important meaning for the further understanding of the evolving carbonate world of the Precambrian.     

论“偶层面上磁化强度的空间梯度”问题

本文指出B.K.Bhattacharyya和K.C.Chan在文献[1]中基本公式（9）的错误,从而“μ是磁化强度分布沿总磁场的方向u的空间梯度”也是错误的。     

下扬子黄桥地区构造特征及其与油气的关系

下扬子黄桥地区位于下扬子对冲构造带上,记录了多期构造变形,并且保存了较完整的中古生界海相地层,发育多 套烃源岩,具有良好的油气前景。文中采用三维地震工区构造精细解释及运动学正演模拟方法,分析了研究区的构造变形 特征及其变形过程,结果表明黄桥地区中上古生界和下中生界地层在自SE向NW的挤压应力作用下,以前展式发育了多个 断层相关褶皱构造。在构造前锋,断层相关褶皱反向发育,于工区中部偏北西侧形成明显的对冲构造,地层遭受强烈的褶 皱、隆升,构造总缩短率高于47.55%;之后,区域应力场转变为NW-SE向拉张,工区南东侧发育负反转断层,且自南西 向北东负反转程度逐渐变大,由上正下逆复合型转变为似正断层型;此次拉伸正断作用导致断层回滑,形成上白垩统填充 的半地堑构造,对冲构造也受到了一定的破坏。在此基础上,结合前人油气成藏的研究成果,探讨了黄桥构造变形与其油 气成藏的关系,并提出了对冲构造背斜、对冲构造下盘及南东侧的广陵构造背斜为黄桥地区的有利油气区块。研究结果对 进一步认识黄桥及下扬子区的构造变形特征和演化、以及该区的油气勘探提供一定的参考。     

基于共反射点叠加的观测参数选择

在地震勘探数据处理中，通常采用的共中心点(CMP)叠加技术对于复杂的地质构造难以取得好的效果。若通过从道间距、覆盖次数等观测参数与目标层地层倾角、共反射点离散度的关系出发，并结合模型正演，对复杂构造地层的观测参数选取进行优化，使之能对目标层尤其是逆断层下降盘的阴影层段进行较充分的“照明”，即可得到共反射点叠加剖面，这对共反射点叠加是一项有益的研究。