Investigating a time‐shift extended imaging condition in a Kirchhoff pre‐stack depth migration algorithm

Extracting true amplitude versus angle common image gathers is one of the key objectives in seismic processing and imaging. This is achievable to different degrees using different migration techniques (e.g., Kirchhoff, wavefield extrapolation, and reverse time migration techniques) and is a common tool in exploration, but the costs can vary depending on the selected migration algorithm and the desired accuracy. Here, we investigate the possibility of combining the local‐shift imaging condition, specifically the time‐shift extended imaging condition, for angle gathers with a Kirchhoff migration. The aims are not to replace the more accurate full‐wavefield migration but to offer a cheaper alternative where ray‐based methods are applicable and to use Kirchhoff time‐lag common image gathers to help bridge the gap between the traditional offset common image gathers and reverse time migration angle gathers; finally, given the higher level of summation inside the extended imaging migration, we wish to understand the impact on the amplitude versus angle response. The implementation of the time‐shift imaging condition along with the computational cost is discussed, and results of four different datasets are presented. The four example datasets, two synthetic, one land acquisition, and a marine dataset, have been migrated using a Kirchhoff offset method, a Kirchhoff time‐shift method, and, for comparison, a reverse time migration algorithm. The results show that the time‐shift imaging condition at zero time lag is equivalent to the full offset stack as expected. The output gathers are cleaner and more consistent in the time‐lag‐derived angle gathers, but the conversion from time lag to angle can be considered a post‐processing step. The main difference arises in the amplitude versus offset/angle distribution where the responses are different and dramatically so for the land data. The results from the synthetics and real data show that a Kirchhoff migration with an extended imaging condition is capable of generating subsurface angle gathers. The same disadvantages with a ray‐based approach will apply using the extended imaging condition relative to a wave equation angle gather solution. Nevertheless, using this approach allows one to explore the relationship between the velocity model and focusing of the reflected energy, to use the Radon transformation to remove noise and multiples, and to generate consistent products from a ray‐based migration and a full‐wave equation migration, which can then be interchanged depending on the process under study.     

瞬变电磁法深度探测正演模拟研究

针对当前瞬变电磁法探测能力有限的问题,利用一次磁场理论公式,分析相同磁矩下、不同发射线圈一次磁场分布及衰减规律,从而确定其信号深度影响范围;对水平层状介质相同磁矩下、不同发射线圈的二次场进行数值模拟,采用Guptasarma和Singh滤波算法及余弦变换法,利用均匀半空间解析公式进行验证,并对不同正演模型结果进行分析,总结发射磁矩决定探测深度大小。对磁矩相同条件下产生的结果数据进行拟合,得出有效探测深度与线圈边长2倍相当的结论,可为野外矿井生产与工程勘探提供施工布置与定量解释的指导依据,从而提高野外工作效率。     

青藏高原东北缘主要断裂滑动速率及其动力学意义

本文利用青藏高原东北缘地区1991—2015年的GPS速度场资料, 基于弹性球面块体模型获得了区域活动断裂的滑动速率, 并讨论了断裂滑动速率分配的动力学意义。 反演结果表明, 青藏高原东北缘地区主要块体以北东向并兼顺时针旋转运动为主; 区域断裂平均闭锁深度为17 km; 另外, 各主要断裂滑动速率也不尽相同。 其中, 阿尔金断裂、 东昆仑断裂左旋走滑速率为10～12 mm/a, 祁连—海原断裂左旋走滑速率为3～5 mm/a, 鄂拉山断裂、 拉脊山断裂右旋走滑速率为1～3 mm/a。 阿尔金断裂、 祁连—海原断裂、 东昆仑断裂的走滑速率被其端部的山脉隆起和逆冲断裂所吸收和转换, 鄂拉山断裂和拉脊山断裂则起到了调节块体间运动平衡的作用。     

由初至P震相重新定位2017年九寨沟地震序列的主震与ML≥3.0余震并分析发震构造

综合利用区域台站和流动台站(近台)的记录,基于初至P震相重新测定了2017年九寨沟序列M_S7.0主震和M_L≥3.0余震的震源位置,并利用较高精度的定位结果分析余震分布与地震构造的关系,解释发震断裂带的结构.获得的新认识有:(1)九寨沟主震震源深度为16km,位于余震带中段的南缘;余震主要分布深度为4～17km.(2)沿余震带的走向,余震分布与主震同震位错大小的分布明显相关.余震带中段8～16km深度存在的余震稀疏区与同震位错的高值区相吻合,应是发震断裂带主凹凸体的部位,也是主震时应变释放较充分的部位;余震带南东段10～18km深度的余震密集区对应了同震位错的亏损区之一,三次M_L≥5.0余震都发生于此;余震带西北段在5～10km之下既缺少余震,又属同震位错的亏损区,可能与那里多条断裂的交汇或合并造成的构造复杂性有关;余震带中-北西段3～5km深度的也缺少余震,也对应了浅部的同震位错亏损区.(3)证实了九寨沟地震的发震构造为虎牙断裂带北段,同时新揭示出发震断裂带表现为由主断裂和分支断裂构成的、向上分叉的花状结构,尺度约为4.5km宽(最大)、35km长,主断裂朝SW陡倾.这些反映主震破裂可能不只受控于单一的断裂,而有可能是沿主断裂发生主破裂,而沿分支断裂发生次要破裂.另外,本文对发震断裂带结构的分段解释,是遵循构造地质学原理去综合震源排列、震源机制解、地表断层已知位置、相邻剖面断层解释结果等信息的分析结果,而不仅仅依据余震的密集分布进行推断.     

阜阳M_S4.3地震震源深度的确定及地震成因分析

2015年3月14日在安徽阜阳地区发生了M_S4.3地震,随后发生3月23日M_s3.6余震.主震造成2人死亡13人受伤.房屋倒塌155间,受损1万多间.主震震级不大,而造成的灾害巨大.本文使用CAP方法反演了两次地震的震源机制解和震源深度,结果显示两次地震的震源机制解和深度一致.主震的机制解节面Ⅰ走向110°,倾角75°,滑动角—10°;节面Ⅱ走向202°,倾角80°,滑动角—164°;矩震级M_w4.3,余震矩震级M_w3.7,反演最佳深度均为3 km.最佳深度时波形拟合相关系数较高,表明反演结果是可靠的.使用sPn和sPL深度震相进一步分析了两次地震的震源深度.结果显示,选取的7个台站的sPn震相与Pn震相的平均到时差为1 s,对应的震源深度为3 km.震中距为36 km的利辛台的sPL震相与Pg震相到时差约为1.1 s,对应震源深度约3~4 km之间.两种深度震相分析的震源深度与CAP方法的结果一致,表明本文给出的阜阳地震震源深度为3 km左右基本是可靠的.本次地震造成较大灾害的原因很可能与地震震源较浅有关.阜阳地区地壳结构相对稳定,地质构造演化形成3 km厚的沉积层,本次地震可能是区域应力作用下发生在沉积层里的一次地震.     

MATCH AND LOCATE FOR SMALL EVENT DETECTION OF NINGXIA SHIZUISHAN EARTHQUAKE SWARM AND INVESTIGATION OF ITS SEISMOGENIC FAULT

In this paper, a new method for small event detection named Match & Locate(M&L)is used to detect and locate the small earthquakes that are missing in the catalogue of the February 28, 2014 Shizuishan M_L4.4 earthquake swarm. A total of 34 earthquakes were detected, which is nearly 3 times as much as the number(13)of earthquakes from Ningxia seismic network. The number of earthquake swarm sequence is greatly increased. Then, it provides the possibility for the fine study of the earthquake swarm activity and seismogenic fault. The best double couple solution of the main shock obtained by the cut and paste method is strike 354°, dip 70° and slip 166° for nodal plane I, and strike 89°, dip 77° and slip 21°for nodal plane Ⅱ. The main shock is a dextral strike-slip earthquake with a small amount of thrust component. And, the depth of the main shock is 7~8km, which is a shallow earthquake, derived from the results of the double difference relocation and the best fitting depth of focal mechanism. Together with the results of deep 3-D seismic tomography of the Yinchuan Basin, our results show that the main shock and the largest aftershock more likely occurred in the upper crust, and the rest of earthquakes mainly occurred at the bottom of sedimentary layer or on the top of the upper crust crystallization basement. We find some interesting phenomena on the pattern of time-space evolution of the earthquake swarm. The distribution of earthquake swarm is in the near north-south direction. Aftershocks are mainly concentrated in the north region of the main shock, which show an obvious trend of extending gradually from the south to the north. Also, the result shows the general trend of shallower focal depth with the development of aftershocks to the north. The results of distribution and depth profile of the earthquake swarm and the focal mechanism of the main shock all show that the sequence probably occurred in the fault at the east foot of Helan Mountain with an eastward dip and a larger dip angle. Surface projection image of the earthquake sequence shows that the epicenter distribution extends northward from the northern end of the fault. This may suggest that the deep part of the fault is likely to extend northward.     

云南洱海表层沉积物中C25高支链类异戊烯烃(HBIs)的特征及其生物地球化学指示意义

C₂₅高支链类异戊烯烃(C₂₅HBIs)通常被认为是指示硅藻的生物标志化合物,近年来在长江下游湖泊沉积物中被广泛检出。尽管云南湖泊分布广泛,但尚未有C₂₅HBIs检出的相关报道。本研究通过分析洱海不同湖区表层沉积物,首次指出其存在C₂₅HBIs化合物,并通过分析不同C₂₅HBIs化合物的空间分布特征解析其生物地球化学指示意义。结果显示,洱海表层沉积物中存在C_25:1 HBI和C_25.2HBI化合物,但均未检测到同分异构体。洱海沉积物中C₂₅HBIs浓度为11.18μg/g(TOC),显著高于富营养化的太湖,与澄湖相当。C_25:1 HBI浓度表现出南高北低的空间分布特征,而C_25:2HBI浓度则相反,因此C_25.1HBI可能表征了洱海南部湖区的优势硅藻。在C_25.2HBI浓度较低(<0.4μg/g)的沉积物中,C_...     

大庆地区地震精定位与震源机制解分析

依托大庆城市活动断层探测项目,开展大庆区域地震活动性研究。筛选黑龙江省及周边662个地震事件的震相到时,修订研究区域走时模型。对研究区148个地震运用交切法和双差法重新定位,对比发现,双差定位后地震定位精度显著提高,地震空间分布更加集中,条带状更为明显,震源深度表明,研究区孕震层基本位于地壳中上部。利用震源机制解资料,对工作区现代构造应力场进行分析,认为本区域构造应力场主压应力方向为NE向,以水平走滑作用为主。     

内蒙古阿拉善地区几次中等强度地震震源深度测定

震源深度对发震构造研究具有重要意义,内蒙古阿拉善地区近年来地震活动频繁,发生数次中强地震和大量中小地震。选取阿拉善地区8个中等强度地震,使用Pn、Pg联合测定深度方法(PTD)、sPn和Pn走时差方法和CAP方法测定震源深度,3种方法计算结果的一致,重新测定的震源深度范围为(18±7)km。     

采用CAP方法反演2016年广西苍梧M_S 5.4地震震源机制解

北京时间2016年7月31日广西梧州市苍梧县发生M_S 5.4地震,基于海南地震台网数字波形资料,采用CAP方法反演震源机制解。结果显示,此次M_S 5.4地震震源深度较浅,最佳深度为5.1 km,其中节面Ⅰ参数为:走向340°,倾角37°,滑动角-18°;节面Ⅱ参数为:走向85°,倾角79°,滑动角-125°。初步推断苍梧M_S 5.4地震破裂面运动以走滑为主,兼有正断性质,反演参数与中国地震台网中心结果较为一致。