Reverse time migration with Gaussian beams using optimized ray tracing systems in transversely isotropic media

Prestack depth migration is a key technology for imaging complex reservoirs in media with strong lateral velocity variations. Prestack migrations are broadly separated into ray-based and wave-equation-based methods. Because of its efficiency and flexibility, ray-based Kirchhoff migration is popular in the industry. However, it has difficulties in dealing with the multi-arrivals, caustics and shadow zones. On the other hand, wave-equation-based methods produce images superior to that of the ray-based methods, but they are expensive numerically, especially methods based on two-way propagators in imaging large regions. Therefore, reverse time migration algorithms with Gaussian beams have recently been proposed to reduce the cost, as they combine the high computational efficiency of Gaussian beam migration and the high accuracy of reverse time migration. However, this method was based on the assumption that the subsurface is isotropic. As the acquired azimuth and maximum offsets increase, taking into account the influence of anisotropy on seismic migration is becoming more and more crucial. Using anisotropic ray tracing systems in terms of phase velocity, we proposed an anisotropic reverse time migration using the Gaussian beams method. We consider the influence of anisotropy on the propagation direction and calculate the amplitude of Gaussian beams with optimized correlation coefficients in dynamic ray tracing, which simplifies the calculations and improves the applicability of the proposed method. Numerical tests on anisotropic models demonstrate the efficiency and accuracy of the proposed method, which can be used to image complex structures in the presence of anisotropy in the overburden.     

利用互连网发布海量地震应急信息方案

通过比较几种Internet GIS不同开发方案,详细探讨了ArcIMS地图服务软件平台的特点和功能。在Internet上开发高效的海量地理信息数据发布系统时候,涉及到超大容量的影像/图像数据,提出了解决这个问题的方案,并以构建地震应急指挥系统中的信息发布子系统为例,介绍了使用ArcIMS进行海量信息发布的实现方案。     

莱州湾海域郯庐断裂带活断层探测

利用浅地层剖面仪对郯庐断裂带莱州湾段进行了活断层探测,发现郯庐断裂带主干断裂在第四纪晚期以来具有明显的活动,继承了晚第三纪以来的主要构造活动特点,仍是这一区域的主导性构造. 西支KL3断裂由多条高角度正断裂组成,最新活动时代为晚更新世晚期至全新世早期,受到一系列错断晚更新世晚期沉积的北东或近东西向断裂的切割;东支龙口断裂由两段右阶斜列的次级断层组成,沿断裂带不但有明显的晚第四纪断错活动,而且还发育北北东向晚第四纪生长褶皱,表现出明显的晚更新世晚期至全新世活动特征. 在山东陆地区也发现了与龙口断裂相对应的安丘——莒县断裂,安丘段由一系列右阶斜列的次级断层组成. 从安丘向北至莱州湾凹陷,郯庐断裂带东支活断层构成了一条右旋单剪变形带,每一个次级活断层段相当于带内理论上次级压剪面,在第四纪晚期以来仍以右旋走滑活动为主要特征.      

深切怀念地球物理研究所老所长、著名地球物理学家顾功叙院士

2022年1月14日是国际著名地球物理学家、我国勘探地球物理事业的奠基人与开拓者之一、也是地震预报研究的开创者之一顾功叙先生逝世30周年纪念日。1992年的今天,让人尊敬爱戴的、德高望重的地球物理研究所名誉所长顾功叙先生因病与世长辞。顾老虽然已离开我们整整30年了,但是他崇高的品德风范、高瞻远瞩的学术思想以及事业上的杰出功绩,至今依然被我们所铭记,并以崇敬的心情深切缅怀他。     

1958-2017年太湖蒸发量年际变化趋势及主控因子

湖泊蒸发对气候变化非常敏感,是水文循环响应气候变化的指示因子,因此研究湖泊蒸发的控制因素,对于理解区域水文循环有重要意义.本文利用太湖中尺度涡度通量网避风港站观测数据校正JRA-55再分析资料,驱动CLM4.0-LISSS模型,并利用2012-2017年涡度相关通量数据和湖表面温度数据检验模型模拟蒸发结果,验证了该模型在太湖的适用性;估算了1958-2017年间太湖的湖面蒸发量,并利用Manner-Kendall趋势检验分析了湖面蒸发的变化趋势,寻找太湖实际蒸发的年际变化的主控因子.结果如下:校正后的JRA-55再分析资料模拟的太湖蒸发与观测值之间存在季节偏差,但是季节偏差在年尺度上相互抵消,再分析资料可用于年际尺度太湖蒸发变化的模拟;1958-2017年间太湖蒸发量以1977年为界,先下降(-3.6 mm/a),后增加(2.3 mm/a);多元逐步回归结果表明,向下的短波辐射是太湖1958-2017年间太湖蒸发变化的主控因子,向下的长波辐射、气温、比湿也对湖泊蒸发年际变化有一定影响,但是风速对蒸发量的年际变化影响不大.     

Revisiting Force Model Error Modeling in GRACE Gravity Field Recovery

Surveys in Geophysics - The gravity field recovery from GRACE (Gravity Recovery and Climate Experiment) mission data is contaminated by both observation noise and dynamic force errors, especially...     

Mineralogical and geochemical study of a newly discovered mafic granulite, northwest China: Implications for tectonic evolution of the Altay Orogenic Belt

Abstract A mafic granulite body was newly discovered in the Altay Orogenic Belt, northwest China. The rocks comprise a suite of coarse‐grained and fine‐grained granulites. Orthopyroxenes (hypersthenes) in the rocks have high X_Mg and low Al₂O₃ contents, whereas clinopyroxenes have low TiO₂ and Al₂O₃ contents. Amphiboles and biotites have a high Mg/(Mg + Fe²⁺) ratio and low contents of F and Cl. The peak metamorphic pressure–temperature (P–T) conditions are estimated as 750–780°C and 6–7 kbar, and retrograde P–T conditions are in the range of 590–620°C and 2.3–3.7 kbar, indicating significant decompression. Metamorphic reactions and P–T estimates define a clockwise P–T path. Geochemically, the rocks are high in Mg/(Mg + Fe) and Al₂O₃, depleted in U, Th, K and Rb, and characterized by light rare earth element enrichment and a weak positive Eu anomaly. The Altay mafic granulite shows depleted Nb, P and Ti contents in the mid‐oceanic ridge basalt normalized spider diagram. The geochemical characteristics suggest that the protolith of the Altay mafic granulite was calc‐alkaline basalt and andesite with an island‐arc affinity. The rock has a high ¹⁴³Nd/¹⁴⁴Nd ratio with ?Nd(0) > 0, indicating derivation from a mantle‐depleted source. In the present study, a two‐stage model for the evolution of the Altay mafic granulite is proposed: an early stage in which calc‐alkaline basalt and andesite with island‐arc affinity were subducted into a deeper level of the crust and subjected to granulite‐facies metamorphism generating the mafic granulite, followed by the later stage exhumation of the system into the upper crust by the late Paleozoic thrusting.     

Preliminary results of the close-off depth and the stable isotopic records along a 109.91 m ice core from Dome A,Antarctica

A 109.91 m ice core was recovered from Dome A (or Dome Argus), the highest ice feature in Antarctica, during the 2004/05 austral summer by the 21st Chinese National Antarctic Research Expedition (CHINARE-21). Both methane profile along the core and firn densification model calculation suggest that the close-off depth is at about 102.0 m with an ice age about 4200 a. Stable isotopes (δ18O and δD) of the chips samples produced during each run of ice core drilling at Dome A, together with those of the other co...     

基于二次场的可控源电磁法三维有限元-无限元数值模拟

本文提出了一种新的混合有限元-无限元三维可控源电磁法（CSEM）问题快速高精度正演模拟算法.首先从电场双旋度方程出发,推导了水平电偶极子源的二次场边值问题,采用无限元代替截断边界条件和有限元离散内部计算区域的新策略,达到减小计算区域的目的,基于并行直接求解技术,实现多源CSEM问题的快速精确求解.其次,通过层状解析模型测试,一方面验证了新算法的正确性,另一方面通过与其他三种已知CSEM问题求解策略进行对比,表明了本文提出的基于二次场有限元-无限元算法具有离散区域小、求解速度快和计算精度高等优点.最后,通过3D模型计算,清晰直观地模拟了场源阴影效应,为野外数据的处理与解释提供指导.

     

超声波处理对微囊藻群体光合活性和沉降过程的影响

关注微囊藻群体在超声处理下光合活性变化及超声解除后的浮力恢复情况,可以为超声波控藻技术提供理论依据和数据支撑,具有重要的实践意义.本研究采用处理后水柱中叶绿素a浓度的垂向分布比例和叶绿素荧光变化情况分析微囊藻群体在超声波处理下的光合活性变化和沉降过程.结果表明,适量的超声处理(35 kHz、0.0353 W/cm~3)能在避免破裂藻细胞的同时,显著抑制其光合活性;60 s的处理造成了45.5%的光合活性被抑制.然而,在适宜生长条件下,所有实验组藻细胞中受抑制的光合活性均在24 h内恢复至对照组的80%左右.此外,在上述超声条件下,5 s的超声处理能使水柱中叶绿素a浓度在短时间内(0.5 h)的去除率达79.5%.然而,当处理时间小于30 s时,大于90%的沉降藻细胞可在超声解除后的72 h内恢复浮力上浮;而当处理时间延长至60 s以后,藻细胞基本丧失了上浮能力.通过分析发现,超声处理后微囊藻群体的粒径分布对藻细胞沉降及上浮过程起决定性作用,并且还发现微囊藻群体在超声处理时首先表现为藻细胞失去浮力下沉和光合系统受损,进而发生大群体振散.