Morphology,sedimentary infill and depositional environments of the Early Quaternary North Sea Basin (56°–62°N)

The North Sea Basin has been subsiding during the Quaternary and contains hundreds of metres of fill. Seismic surveys (170 000 km²) provide new evidence on Early Quaternary sedimentation, from about 2.75 Ma to around the Brunhes-Matuyama boundary (0.78 Ma). We present an informal seismic stratigraphy for the Early Quaternary of the North Sea, and calculate sediment volumes for major units. Early Quaternary sediment thickness is > 1000 m in the northern basin and >700 m in the central basin (total about 40 000 km³). Northern North Sea basin-fill comprises several clinoform units, prograding westward over 60 000 km². Architecture of the central basin also comprises clinoforms, building from the southeast. To the west, an acoustically layered and mounded unit (Unit Z) was deposited. Remaining accommodation space was filled with fine-grained sediments of two Central Basin units. Above these units, an Upper Regional Unconformity-equivalent (URU) records a conformable surface with flat-lying units that indicate stronger direct glacial influence than on the sediments below. On the North Sea Plateau north of 59°N, the Upper Regional Unconformity (URU) is defined by a shift from westward to eastward dipping seismic reflectors, recording a major change in sedimentation, with the Shetland Platform becoming a significant source. A model of Early Quaternary sediment delivery to the North Sea shows sources from the Scandinavian ice sheet and major European rivers. Clinoforms prograding west in the northern North Sea Basin, representing glacigenic debris flows, indicate an ice sheet on the western Scandinavian margin. In the central basin, sediments are generally fine-grained, suggesting a distal fluvial or glacifluvial origin from European rivers. Ploughmarks also demonstrate that icebergs, derived from an ice sheet to the north, drifted into the central North Sea Basin. By contrast, sediments and glacial landforms above the URU provide evidence for the later presence of a grounded ice sheet.     

Parsimonious wave‐equation travel‐time inversion for refraction waves

We present a parsimonious wave‐equation travel‐time inversion technique for refraction waves. A dense virtual refraction dataset can be generated from just two reciprocal shot gathers for the sources at the endpoints of the survey line, with N geophones evenly deployed along the line. These two reciprocal shots contain approximately 2N refraction travel times, which can be spawned into refraction travel times by an interferometric transformation. Then, these virtual refraction travel times are used with a source wavelet to create N virtual refraction shot gathers, which are the input data for wave‐equation travel‐time inversion. Numerical results show that the parsimonious wave‐equation travel‐time tomogram has about the same accuracy as the tomogram computed by standard wave‐equation travel‐time inversion. The most significant benefit is that a reciprocal survey is far less time consuming than the standard refraction survey where a source is excited at each geophone location.     

Towards quantitative evaluation of gas injection using time‐lapse seismic data

Of particular concern in the monitoring of gas injection for the purposes of storage, disposal or improved oil recovery is the exact spatial distribution of the gas volumes in the subsurface. In principle this requirement is addressed by the use of 4D seismic data, although it is recognized that the seismic response still largely provides a qualitative estimate of moved subsurface fluids. Exact quantitative evaluation of fluid distributions and associated saturations remains a challenge to be solved. Here, an attempt has been made to produce mapped quantitative estimates of the gas volume injected into a clastic reservoir. Despite good results using three accurately repeated seismic surveys, time‐delay and amplitude attributes reveal fine‐scale differences though large‐scale agreement in the estimated fluid movement. These differences indicate disparities in the nature of the two attributes themselves, which can be explained by several possible causes. Of most impact are the effects of processing and migration, wave interference effects and noise from non‐repeatability of the seismic surveys. This subject highlights the need for a more careful consideration in 4D acquisition, amplitude processing and use of true amplitude preserving attributes in quantitative interpretation.     

Repeatability analysis of land time‐lapse seismic data: CO2CRC Otway pilot project case study

Time‐lapse seismics is the methodology of choice for remotely monitoring changes in oil/gas reservoir depletion, reservoir stimulation or CO₂ sequestration, due to good sensitivity and resolving power at depths up to several kilometres. This method is now routinely applied offshore, however, the use of time‐lapse methodology onshore is relatively rare. The main reason for this is the relatively high cost of commercial seismic acquisition on land. A widespread belief of a relatively poor repeatability of land seismic data prevents rapid growth in the number of land time‐lapse surveys. Considering that CO₂ sequestration on land is becoming a necessity, there is a great need to evaluate the feasibility of time‐lapse seismics for monitoring. Therefore, an understanding of the factors influencing repeatability of land seismics and evaluating limitations of the method is crucially important for its application in many CO₂ sequestration projects. We analyse several repeated 2D and 3D surveys acquired within the Otway CO₂ sequestration pilot project (operated by the Cooperative Research Centre for Greenhouse Technologies, CO2CRC) in Australia, in order to determine the principal limitations of land time‐lapse seismic repeatability and investigate the influence of the main factors affecting it. Our findings are that the intrinsic signal‐to‐noise ratio (S/N, signal to coherent and background noise levels) and the normalized‐root‐mean‐square (NRMS) difference are controlled by the source strength and source type. However, the post‐stack S/N ratio and corresponding NRMS residuals are controlled mainly by the data fold. For very high‐fold data, the source strength and source type are less critical.     

ICESat-2水深提取方法及其反演应用

相比于传统水深测量方法, 结合主动激光卫星ICESat-2数据与被动光学遥感影像的主被动融合水深反演方法具有低成本、覆盖面积广、不受地域限制、无需地面实测数据等优点.研究选取澳大利亚大堡礁中部的John Brewer Reef和南部的Fitzroy Reef为实验区, 开展主被动融合水深反演实验.针对ICESat-2光子点群信噪比低、光子点水深信息提取和改正难的问题, 本文提出基于光子点群高斯分布特征的海面光子点群提取方法和基于深度自适应DBSCAN的海底光子点群提取方法, 并顾及海面倾斜引起的入射角偏移, 改进了Parrish水体折射改正几何模型.实验表明, 基于光子点群高斯分布特征的海面光子点群提取方法和基于深度自适应DBSCAN的海底光子点群提取方法可以高效准确地提取出有效的海面、海底光子点群; 经过地球物理改正和基于改进模型的水体折射改正后, 水深变浅, 且改正项与改正前水深呈线性关系, 线性系数约为-0.2550.实验将提取并改正后的ICESat-2光子点水深与经大气校正、耀斑校正、空间域滤波等预处理后的Sentinel-2影像进行空间匹配构建水深反演样本集, 并利用单波段模型(SB)、比值模型(BR)、Lyzenga多项式模型(LP)、二次多项式模型(QP)、三次多项式模型(CP)、支持向量回归模型(SVR)、多层感知器模型(MLP)、随机森林回归模型(RF)等水深反演模型对实验区全域水深进行反演.通过对比分析多种水深反演模型在不同地理区域和ICESat-2水深样本分布条件下的表现, 研究为不同场景下的主被动融合水深反演策略提供了一定的科学依据.

     

从地震层析速度确定原位岩土力学参数

地震纵横波层析成像速度与岩土弹性力学参数有关,多道面波分析（MASW）采用类似二维反射地震处理,得到共中点互相关道集（CMPCC）提取面波频散曲线,由面波频散曲线反演二维速度结构,在岩土结构分层等方面有良好的效果。通过同一剖面的折射地震波层析成像得到纵波速度结构,对岩土层的动弹模量和泊松比等力学参数进行分析,为下一步地质工程设计及岩土力学试验提供参考。     

南海北部新生界地层的远偏移距折射波速度特征

南海北部为张裂大陆边缘,在新生代时期经历了复杂的地质演化,具有良好油气勘探前景.本文首次利用南海北部某区域长排列地震数据的远偏移距折射波速度计算方法获得南海北部地层层速度,并进行了分析,此外,还提出了非初至折射波的概念.本方法在共偏移距剖面上识别、拾取折射波层位,得到折射波走时的横向连续变化数据,利用不同共偏移距同一折射波层位的到达时差计算该折射层位的速度结构.该方法在南海北部陆坡某区域地震资料的应用,揭示了新生界地层层速度的整体结构,反映了新生代不同时代地层的速度变化特征,为进一步利用折射层速度资料探讨南海新生代地质演化提供了新的视角.研究表明,与常规初至折射波法相比,远偏移距折射波能够提供更多的地质信息;该方法的初步应用能够对远偏移距折射波的研究提供进一步的视角.     

基于相邻虚拟道叠加的超虚折射干涉法及其在广角OBS折射波增强中的应用

广角地震的远偏移距折射初至含有地下深层的信息,提高远偏移距折射波信噪比能够有效提高初至拾取精度,对于深部结构的层析速度建模十分有利.超虚折射干涉法基于干涉原理对远道折射波进行增强,它在炮点和检波点都非常密集的情况下效果较好.然而,在应用于台站间距较大且环境噪声较强的广角海底地震仪（OBS）观测时,该方法对折射波的增强能力不足,而且容易产生虚假波形,造成增强后的折射波信噪比仍然较低.针对这种情况,本文提出基于相邻虚拟道叠加的超虚折射干涉法,通过叠加相邻虚拟道来提高远道与近道互相关的准确度,以达到稳定增强远道折射波信噪比的目的.理论实验和实际资料测试均显示,基于相邻虚拟道叠加的超虚折射干涉法在台站间距较大和信噪比较低的情况下能够准确构建虚拟道,且增强后的波形同相轴连续程度和信噪比均较高,有利于拾取高精度的折射初至到时.本文方法也可用于增强陆上炮间距较大的广角地震数据折射波.

     

煤田三维地震叠前深度偏移技术及其应用

三维地震叠前深度偏移可实现反射点的正确空间归位和真正的共反射点叠加,减小菲涅尔带的影响范围,大大提高地震资料的分辨率。论述了三维地震叠前深度偏移的成像原理和关键技术,并应用于实际地震资料的处理中。结果显示,小断点显示较叠前时间偏移剖面更为清晰,提高了煤田三维地震探测细微构造的能力,取得了较好的地质效果。     

CTIO1m反射望远镜的色折射改正和星等差

利用ＣＴＩＯ（ＣｅｒｒｏＴｏｌｏｌｏＩｎｔｅｒ－ＡｍｅｒｉｃａｎＯｂｓｅｒｖａｔｏｒｙ）１ｍ反射望远镜观测球状星团ωＣｅｎ区域的Ｖ和Ｂ波段的照相底片，分析研究了色折射改正；描述了测定色折射改工项系数的方法，得到该望远镜在Ｖ和Ｂ波段照相观测的色折射改正项的系数分别为１．３５２１″和４．５４０３″．随后分析讨论了与星等有关的系统误差，结果表明，对大部分照相压片，星等差是很明显的．在利用底片常数归算时，星等的线性和二次项必须包括在内．