基于含气泡液体声波方程的海底冷泉数值模拟

海底冷泉羽状流与海底天然气水合物的分布密切相关,对水合物稳定带的边界具有指示作用,是未来能源勘探的重要领域.研究海底冷泉羽状流的地震响应特征,对确定天然气水合物的储集区域及成藏环境等均有重要意义.当前获得海底冷泉羽状流的地震响应主要通过数值模拟进行,然而该过程所依据的含气泡介质声速模型及随机介质理论不能完整描述海底冷泉的物理性质,采用的声波方程也不适用于高频地震波数值模拟.为了准确地实现海底冷泉羽状流地震波数值模拟,精确分析其地震响应特征,本文提出利用Keller-Miksis气泡振动模型来描述气泡在声波作用下的运动状态,同时考虑气泡间的相互作用,建立海底冷泉气泡模型.在此基础上,本文创新性地采用含气泡液体声波方程进行海底冷泉高频地震波数值模拟.数值模拟结果表明,本文提出的方法能够实现海底冷泉羽状流地震响应的高精度数值模拟.     

基于含气泡液体声波方程的海底冷泉数值模拟

海底冷泉羽状流与海底天然气水合物的分布密切相关,对水合物稳定带的边界具有指示作用,是未来能源勘探的重要领域.研究海底冷泉羽状流的地震响应特征,对确定天然气水合物的储集区域及成藏环境等均有重要意义.当前获得海底冷泉羽状流的地震响应主要通过数值模拟进行,然而该过程所依据的含气泡介质声速模型及随机介质理论不能完整描述海底冷泉的物理性质,采用的声波方程也不适用于高频地震波数值模拟.为了准确地实现海底冷泉羽状流地震波数值模拟,精确分析其地震响应特征,本文提出利用Keller-Miksis气泡振动模型来描述气泡在声波作用下的运动状态,同时考虑气泡间的相互作用,建立海底冷泉气泡模型.在此基础上,本文创新性地采用含气泡液体声波方程进行海底冷泉高频地震波数值模拟.数值模拟结果表明,本文提出的方法能够实现海底冷泉羽状流地震响应的高精度数值模拟.     

利用小波变换分析羽状流模型偏移剖面尺度特征

甲烷羽状流备受关注,因为它提供了海底气体泄漏的直接证据,并对天然气水合物的勘探识别起到了间接指示作用.为了深入研究冷泉活动区气泡羽状流的地震响应特征,根据已有研究成果——羽状流模型及其偏移剖面,通过小波变换分析单道地震信号及羽状流剖面的尺度特征.单道地震信号最佳尺度曲线表明不同气含量对应的最佳尺度也不同,而且不同深度的最佳尺度变化较大;羽状流最佳尺度剖面能够显示出羽状流区域的轮廓.羽状流气含量与最佳尺度呈较好的正相关,即随着气含量的增加,最佳尺度也随之增大.     

南海北部琼东南海域活动冷泉特征及形成模式

近年来,活动冷泉的研究越来越受到关注.本文利用多波束数据、多道地震数据以及底质取样结果研究琼东南海域活动冷泉系统,分析活动冷泉的羽状流特征、海底地貌与底质特征以及流体活动构造特征.多波束水体数据上,观测到多个延伸高度超过750 m的气泡羽状流,海底流体活动非常强烈;多道地震上识别出麻坑、流体运移通道、气烟囱等流体渗漏相关的构造,与其他海域观测到的反射特征不同,羽状流的下方流体运移通道呈强振幅"串珠"反射;重力活塞取样在两个站位上获得浅表层块状天然气水合物.其中一个站位位于活动冷泉附近,天然气水合物赋存于海底以下8 m左右.基于以上三方面的数据,笔者提出了一个用于描述活动冷泉系统的形成模式,游离气通过气烟囱向上运移到达浅层,一部分在天然气水合物稳定带内形成天然气水合物,另一部分穿透天然气水合物稳定带到达海底,形成活动冷泉的羽状流.     

海底冷泉的地震海洋学初探

海底冷泉活动在全球大陆边缘海域广泛分布.与传统利用高频声学方法探测海底冷泉羽状流不同,本文利用常规多道反射地震(地震海洋学)方法对海底边界附近水体进行成像,结合盆地流体逸散结构特征,圈定活动冷泉流体活动发育位置,分析其地震反射特征和流体活动特征.研究表明,活动冷泉流体渗漏、逸散活动在多道反射地震剖面上一般呈羽状、扫帚状和不规则状几何形态,内部反射杂乱,反射振幅偏弱,但也有振幅增强的情况出现,这可能是由于含有较多泥质和细粒颗粒物悬浮.所造成的.冷泉活动一般与盆地内部泥底辟、流体管道、断层和裂隙、气烟囱、海底麻坑和泥火山等流体逸散结构相关,反映了地层内部流体自深部向浅部运移,在海底渗漏、逸散形成了活动冷泉.但分析结果仍需要进一步实地观测和理论模拟等相关研究确认.     

三维各向异性裂缝介质正演模拟的三种交错网格适应性比较及Lebedev方法的改进

研究了三维各向异性裂缝介质的正演模拟的三种不同有限差分法,即标准交错网格(SSG)、旋转交错网格(RSG)和Lebedev方法(LS),详细对比分析这三种交错网格方法在模拟复杂各向异性介质的优势与局限.提出一种新的改进方法,简化了LS有限差分法在任意各向异性介质中的正演模拟.为了模拟三维大规模复杂各向异性介质的地震响应,提出一个优化的正演模拟计算流程：将模型参数分为模型介质参数和模型构造参数.该计算流程适用于三种有限差分法中的任何一种.使用LS方法实现任意三维各向异性裂缝介质地震响应的三维全波场模拟,通过使用三种不同的有限差分模拟方法进行二维和三维模型数值模拟试验,验证了所提出方法有效.     

卡斯卡迪亚大陆边缘活动冷泉流场定量研究

海底冷泉是一种重要的地质活动,当前对冷泉的研究逐渐从定性研究过渡到定量研究.本文利用遥控深潜器视频资料对美国西海岸卡斯卡迪亚大陆边缘活动冷泉进行研究.观察视频资料发现,研究区活动冷泉的释放主要有两种形式,一是甲烷气泡的缓慢逸散,二是含甲烷流体快速喷发形成羽状流.缓慢逸散的甲烷气泡多为透明椭球形,羽状流喷口处多为锥形,流体浓度高时呈白色不透明状,流体边缘处往往为灰色半透明.羽状流流动方向整体向上,内部往往存在复杂的湍流涡.针对研究区中两种释放形式,分别利用不同方法进行流场的定量研究,其中,冷泉羽流的流场定量分析是通过粒子图像测速技术来获得冷泉羽状流瞬时全场流动定量信息,在此基础上试算了羽流的卷吸系数.本文利用图像测速技术对冷泉流场进行定量研究,可为估算甲烷通量等研究提供新的数据来源,对研究冷泉地质过程及其对全球气候变化的影响都具有重要意义.     

天然气水合物冷泉和气泡羽状流研究进展

近些年来,天然气水合物的调查研究及勘探开发都取得了一定进展,但针对与天然气水合物相关的海底冷泉和气泡羽状流(或羽状水体)的研究报道较少.为了更好地进行天然气水合物基础研究,本文首先简要总结了天然气水合物的发展状况;然后阐述了与天然气水合物相关的冷泉和气泡羽状流的发展状况,天然气水合物羽状流识别技术情况,以及地震识别水合物的发展现状;最后针对气泡羽状流提出天然气水合物的研究方向和科学问题,为水合物的识别和勘查提供新的线索和思路.     

海底电缆双检资料正演模拟及特征分析

海底电缆双检采集是滩浅海地区地震勘探的重要手段.目前,国内外对双检资料的研究主要集中在鸣震压制方面,而对双检资料中有效波和鸣震的传播规律和陷频特征等研究较少.为了理清有效波和鸣震在水检和陆检中不同响应产生的原因,本文针对海底电缆双检资料开展正演模拟研究,利用二维声波方程的一阶应力—速度方程中的压力和垂向速度分量表示水检和陆检进行正演模拟,并采用自由表面的应力镜像法处理气水边界,形成了双检资料针对性正演模拟技术;分析有效波和鸣震在双检中的波场传播特征,以及不同水深双检资料的陷频规律,为双检资料鸣震压制等处理提供理论依据.     

声波介质一次散射波场高斯束Born正演

Born正演是一种常用的地震波场正演模拟方法,也是线性化地震反演的理论基础.在实际应用时,Born正演通常结合常规的地震射线方法进行实现.为了克服常规地震射线方法的弊端,并且保证地震波场的模拟精度和计算效率,本文提出了一种基于高斯束的一阶散射波场Born正演方法.该方法分为两个环节:首先,我们利用高斯束的走时和振幅信息将地下散射点处的反射率映射为地表束中心位置处的局部平面波;然后,我们利用逆倾斜叠加将局部平面波转化为接收点处的时空域散射波场.在具体的实施过程中,我们提出一种以wavelet-bank方式实现的局部平面波合成方法,同现有的算法相比,可以在保持计算精度的同时,大大减少计算时间;此外,我们还利用最速下降法优化了高斯束的迭代循环过程,进一步提高了Born正演的计算效率.两个模型的应用效果证明,本文所提出的高斯束Born正演方法可以精确、高效的实现声波介质一次散射波场的正演模拟,为三维大规模地震波场的正演问题提供了一种切实可行的实现方案.     

Simulation of the effect of methane bubble plumes on vertical mixing in Mono Lake

A one-dimensional vertical mixing model modified for application to hypersaline Mono Lake reproduced mixed layer dynamics well but hypolimnetic heating was underestimated. One possible source of increased hypolimnetic heating is vertical mixing caused by bubble plumes of methane rising from the sediments. Estimates of vertical mixing from methane seepage in Mono Lake were made with the inclusion of a bubble plume algorithm. A methane ebullition rate three hundred times greater than the maximum estimate for Mono Lake was required to simulate the observed hypolimnetic heating.     

Mantle plume,large igneous province and continental breakup——Additionally discussing the Cenozoic and Mesozoic mantle plume problems in East China

Based on the former workers study results such as numerical simulation of fluid mechanics,seismic tomography of the whole earth and igneous rocks,the basic characteristics of mantle plumes are summarized in detail,namely the mantle plume,from the D″layer near the core-mantle bouundary(CMB)of 2900 km deep,is characterized by the spape of large head and thin narrow conduit,by the physical property of high temperature and low viscosity.The LIP(large igneous province)is the best exhibition when the mantle plume ascends to the surface.According to the basic characteristics of the mantle plumes and the LIP,as well as the temporal-spatial relationships between the mantle plume and continental breakup,the detailed research on petrology,geochemistry,temporal-spatial distribution,tectonic background of the Cenozoic-Mesozoic igneous rocks and gravity anomaly distribution in East China has been done.As a result,the Mesozoic igneous rocks in Southeast China should not be regarded as an example of typical LIP related to mantle plumes.for their related characteristics are not consitent with those of the typical LIPs related to mantle plumes.The Cenozoic igneous rocks in Northeast China have no the typical characteristics of mantle plumes and hotspots,so the Cenozoic volcanism in Northeast China might have no the direct relationships with the activity of mantle plumes.     

南海东北部陆坡冷泉系统的浅地层剖面特征与分析

本文通过对南海东北部陆坡区浅地层剖面资料的处理,结合横跨陆坡区的长地震测线,总结了该区浅层气聚集、海底流体运移、泥火山和冷泉活动在浅地层剖面中的特征,并对剖面中出现的声学空白带现象做了分类.研究认为浊反射、增强反射等现象,反映了该区域存在广泛的浅层气聚集.对声学空白带成因的分析,认为窄空白带很可能对应海底流体运移的通道,而宽空白带则可能由沿地层分布的浅层气引起,同时提出了自生碳酸盐岩体散落分布造成声学空白的可能性.我们在南海东北部陆坡区发现了泥火山,并讨论了其运移通道及对BSR的破坏.同时发现研究区水体中存在声学异常体,推断为冷泉羽状流,指示该区存在活动的冷泉.本文依托大量的浅地层剖面数据,对南海东北部陆坡区的冷泉系统有了较全面的认识,对冷泉系统剖面特征的成因进行了初步探讨,可为南海冷泉系统研究的进一步开展提供有益参考.     

莺歌海盆地中央坳陷带海底天然气渗漏系统初探

地球物理资料表明，莺歌海盆地区麻坑众多，浅层气发育，在高分辨率地震剖面上可以识别出到达海底的地震模糊带，在3.5 kHz的浅层剖面上可见穿透到海底的气体聚集带，在声呐图像可以识别出麻坑.莺歌海盆地不仅在其边缘隆起带沉积层较薄的地区存在海底天然气渗漏现象，而且在中央坳陷带沉积层巨厚的地区也存在海底天然气渗漏，海底天然气渗漏在莺歌海盆地是广泛分布的.莺歌海盆地广泛分布的海底天然气渗漏是全球海底天然气渗漏系统的重要组成部分，海底天然气渗漏系统排出的大量甲烷气体的一部分可能穿越水层到达海面.     

Flow rates in the axial hot springs of the East Pacific Rise (21°N): Implications for the heat budget and the formation of massive sulfide deposits

The first in-situ measurements of flow rates in submarine hot springs (temperatures between 275° and 350°C) were made in the hot springs of the East Pacific Rise at 21°N during November, 1981. The flow rates ranged between 0.7 and 2.4 m/s. We estimate that the total rates at which heat and mass are transported from the three hot spring areas—OBS, National Geographic, and Southwest Vents—are ca. 2.2×10⁸ W and 150 kg of fluid/s respectively. The very high rate of heat loss precludes the existence of hot springs as steady state features of the East Pacific Rise at 21°N. We estimate that hot springs are active for a maximum of ca. 40,000 years. Chemical and/or mechanical clogging of the vents can drastically reduce the lifetime of the hot springs.Using a simple model of plume development [1], we calculated that less than 3% of the sulfide particles that are entrained in the hydrothermal plumes at 21°N settle from the plumes before dispersal by a lateral submarine current at a height of ca. 250 m above the seafloor [2]. The time that is required to form the sulfide mounds that surround the hot springs by the accumulation of sulfide particles that settle from the plumes is ca. 70–85 years. The current rate at which sulfide particles settle from the hydrothermal plumes and accumulate near the hot springs appears to be too small to lead to the formation of large massive sulfide deposits such as the Mavrouni deposit on Cyprus [3], within the time limits that are suggested by our thermal energy balance calculations.     

Plume-induced geoid anomalies from 2D axi-symmetric temperature- and pressure-dependent mantle convection models

A collection of numerical simulations of 2D axi-symmetric thermal convection is presented here. The aim is to investigate the shape of geoid anomalies and dynamic topography above a plume. The simulation is based on the Boussinesq approximation and infinite Prandtl number and is carried out in the spherical shell with strongly temperature- and depth-dependent Arrhenius-type viscosity. According to the Arrhenius law, plume models with purely depth-dependent rheology are unphysical and should be taken with care. The strongly coupled temperature- and depth-dependent viscosity enables us to better understand the plume's behavior inside the Earth.The topography and geoid anomalies produced from plumes are sensitive to rheology of the mantle and rheology of the plume; both have effects on shape and amplitude of the geoid anomalies. We determined different categories of the geoid which are related to various rheology. Depth-dependent viscosity models show a geoid with a negative sign above the plume, and temperature-dependent viscosity models depict a bell-shaped geoid. We identified different behaviors in the combined model with temperature- and pressure-dependent viscosity.     

Application of diffracted wave analysis to time‐lapse seismic data for CO2 leakage detection

Time‐lapse seismic analysis is utilized in CO₂ geosequestration to verify the CO₂ containment within a reservoir. A major risk associated with geosequestration is a possible leakage of CO₂ from the storage formation into overlaying formations. To mitigate this risk, the deployment of carbon capture and storage projects requires fast and reliable detection of relatively small volumes of CO₂ outside the storage formation. To do this, it is necessary to predict typical seepage scenarios and improve subsurface seepage detection methods. In this work we present a technique for CO₂ monitoring based on the detection of diffracted waves in time‐lapse seismic data. In the case of CO₂ seepage, the migrating plume might form small secondary accumulations that would produce diffracted, rather than reflected waves. From time‐lapse data analysis, we are able to separate the diffracted waves from the predominant reflections in order to image the small CO₂ plumes. To explore possibilities to detect relatively small amounts of CO₂, we performed synthetic time‐lapse seismic modelling based on the Cooperative Research Centre for Greenhouse Gas Technologies (CO2CRC) Otway project data. The detection method is based on defining the CO₂ location by measuring the coherency of the signal along diffraction offset‐traveltime curves. The technique is applied to a time‐lapse stacked section using a stacking velocity to construct offset‐traveltime curves. Given the amount of noise found in the surface seismic data, the predicted minimum detectable amount of CO₂ is 1000–2000 tonnes. This method was also applied to real data obtained from a time‐lapse seismic physical model. The use of diffractions rather than reflections for monitoring small amounts of CO₂ can enhance the capability of subsurface monitoring in CO₂ geosequestration projects.