钻井液侵入海洋含水合物地层的一维数值模拟研究

本文以墨西哥湾水合物区域为背景,利用数值模拟方法研究了过平衡钻井条件下,当钻井液温度高于地层中水合物稳定温度时,水基钻井液侵入海洋含水合物地层的动态过程及其一般性规律.与侵入常规油气地层相比,耦合水合物分解和再形成是钻井液侵入海洋含水合物地层的主要特征.模拟结果表明,钻井液密度、温度和盐度都对侵入过程有影响.在一定条件下,钻井液密度越大,温度和含盐量越高,则钻井液侵入程度越深,热量传递越远,水合物分解程度越大.分解的水气在合适条件下又会重新形成水合物,影响了钻井液进一步侵入.而重新形成的水合物的饱和度甚至可能高于原位水合物饱和度,在井周形成一个"高饱水合物"环带.这一现象归因于钻井液侵入的驱替推挤、水合物分解的吸热以及地层传热的滞后等因素共同作用.在地层物性一定的条件下,高饱水合物环带的出现与否主要受钻井液温度和盐度控制.水合物分解以及高饱水合物环带的出现对井壁稳定和电阻率测井解释有很大影响.因此,为维护井壁稳定、确保测井准确和减少水合物储层伤害,就必须对钻井液密度、温度和滤失量进行严格控制,防止地层中的水合物大量分解.最好采用控制压力钻井(MPD)和深侧向测井方式,同时尽量选用低矿化度的含水合物动力学抑制剂的钻井液体系,采取低温快速循环方式.     

Thermal conductivity measurement of the synthetic samples of bottom sediments containing methane hydrates

The experimental setup is described, which makes it possible to simulate the quartz sand samples, containing methane hydrates, and to measure their thermal conductivity, using a needle probe of constant power. The method and results of measurements at different temperatures and pressures are considered. It is established that under the P-T-conditions close to the equilibrium for methane hydrate, the measurements result in the essential overestimation of the thermal conductivity the samples, i.e., to an anomalous increase in its calculated values. This is because of the dissociation (with the heat consumption) of the part of hydrates near the needle probe under the action of its heater. It is possible to conclude that this feature (the anomalous increase in the calculated values of thermal conductivity) is certain evidence for the presence of a noticeable quantity of hydrates in the sediments. This observation offers a new possibility of utilization of the geothermal method for prospecting the subsea gas hydrate accumulations. Our conclusions are confirmed by the results of measurements in situ of the thermal conductivity of the bottom sediment of the Black Sea [Kutas et al., 2005].     

南海北部神狐海域天然气水合物成藏模式研究

南海北部陆坡神狐海域是我国海洋天然气水合物勘探开发研究的重点靶区,独特的水合物成藏特征,难以利用当前观测到的沉积速率和流体流动条件对其成藏机理进行解释和量化说明,对其形成演化模式和控制因素尚不明确.本文构建了海洋天然气水合物形成演化过程的动力学模型,模型的主控参量为海底沉积速率和水流速率,以此计算了神狐海域天然气水合物聚集演化过程,并与饱和度的盐度测试值进行对比.最后,在研究神狐海域地质构造活动和水合物成藏动力学基础上建立了神狐天然气水合物形成演化模式.认为神狐海域当前的天然气水合物是在上新世末—更新世早期断裂体系水合物基础上继承演化而来的,神狐海域天然气水合物形成演化具有典型的二元模式.第一阶段水合物形成发生在距今1.5 Ma之前构造活动形成的断裂体系中,高达50 m/ka的孔隙水流动携带了大量的甲烷进入水合物稳定带,导致了水合物的快速生成,在4万年内形成了饱和度达20%的甲烷水合物;第二阶段发生在1.5 Ma以来,泥质粉砂沉积使沉积体渗透率骤减,0.7 m/ka的低速率水流使甲烷供给不足,在海底浅层新沉积体中无法生成水合物,仅在水合物稳定带底部有缓慢的水合物继承增长,并因此形成了神狐海域当前观测到的水合物产出特征,而且水合物资源量仍在减少.     

Numerical simulation of accumulation of gas hydrates during sedimentation and compaction of sediments under subaqueous conditions

The mutual influence of compaction of depositing sediments and formation of gas hydrates in a porous medium is examined. Within the framework of the general mathematical model developed in this study, these interrelated processes are described by a nonlinear system of eight partial differential equations and are analyzed using numerical simulation in terms of the model including compaction of depositing sediments, movement of pore fluids, and formation of gas hydrates from dissolved methane accumulated during sedimentation. On the basis of model examples calculated with the use of representative parameters of the system, it is shown that the hydrate formation rate depends on the sedimentation rate and thermal conditions. Gas hydrate formation is shown to influence the pore fluid velocity in the near-surface zone of sediments.     

海洋天然气水合物的地球物理研究(III):似海底反射

对天然气水合物研究中与似海底反射有关的一些观点进行讨论 ,以推动天然气水合物地震研究的认识 .30年的似海底反射研究表明 ,似海底反射仍然是指示天然气水合物沉积存在的最好手段之一 .有似海底反射 ,是可以认为存在天然气水合物的 .虽然存在“游离气带速度模型”与“水合物楔速度模型” ,但似海底反射主要由天然气水合物稳定带底界下方的游离气引起 .BSR上方的天然气水合物、下方的游离气与天然气再循环和含甲烷的流体流有关 .由于天然气水合物稳定带计算控制因素难以准确确定等因素 ,似海底反射与天然气水合物稳定带底界只是近似的对应关系 .需从动态的观点考虑天然气水合物 天然气体系及其与似海底反射的关系 .     

The Oxygen Cycle in Lake Baikal

The oxygen cycle and the balance of free oxygen and oxygen bound in CO₂in Lake Baikal are discussed based on new data on the gas exchange between the lake and atmosphere, O₂consumption in bottom sediments, and the rate of aerobic decomposition in the abyssal zone of the lake.     

琼东南盆地天然气水合物及其成藏模式的海洋可控源电磁研究

南海琼东南盆地是天然气水合物重要远景区之一.由于盆地大部分地区海底地形平缓、地层近于水平,增加了利用地震反射剖面识别似海底反射（BSR,bottom simulating reflector）的难度,从而影响了对水合物的评价.为了进一步开展琼东南盆地水合物调查研究,本文在研究海域进行了海洋可控源电磁探测试验,将自主研发的10台接收机以500 m的间距,投放至水深约为1360 m的海底,完成了一条4.5 km剖面的电磁数据采集.通过对采集的数据进行处理与二维（2D）反演,获得了研究剖面海底的电阻率断面图.反演结果显示,研究区海底60~330 mbsf（meter bottom of seafloor）的地层中,存在多个横向不连续分布的高阻异常体,电阻率介于2~10 Ωm之间;在海底330 mbsf之下,横向上发育了电阻率为2~4 Ωm的3个高阻体.根据研究区热力学条件,本文估算了生物成因气与热成因天然气的水合物稳定带（GHSZ,gas hydrate stability zone）厚度,结合高阻体的分布特征推断了地震剖面上BSR的位置.在此基础上,对反演的电阻率断面进行解释,推断了研究区水合物的分布及游离气运移通道.研究表明,勘探区具有形成天然气水合物矿藏的地质与地球物理条件,其成藏模式可能属于"断层、裂隙输导的下生上储型",水合物的气源为生物成因气.     

Dynamic-stochastical modeling of long-term fluctuations in Lake Baykal levels and Angara River runoff

The results of modeling of the long-term fluctuations in Lake Baikal levels and Angara River runoff are presented. The dynamic-stochastical model was developed and used to obtain analytical relationships between the statistic characteristics of the Baikal levels fluctuations and the parameters of forcing processes, including water inflow into the lake and the effective evaporation from its surface. The parameter λ determining the inertia of the Baikal levels fluctuations, decreased by about three times after Irkutsk Reservoir filling, resulting in an increase in the variance of the lake levels fluctuations. The dependence of variances of the Baikal levels fluctuations and Angara runoff on parameter λ was studied. It is shown that the decreasing in the variance of the Baikal levels fluctuations corresponds to increasing in the variance of Angara river runoff and vice versa. The results of numerical modeling of Baikal levels fluctuations show that water withdrawal from inflow into the lake has an adverse effect on Angara runoff regime.     

水平井开采南海神狐海域天然气水合物数值模拟

2007年5月,南海北部神狐海域的实地钻探结果表明该区域海底存在大量天然气水合物,其作为未来我国潜在的可开发能源的调查和资源评价工作正在展开.本文以SH7站位的钻探、测井数据为基础,建立了实际水合物藏分层地质模型,主要包括上盖层、水合物层和下盖层,其中上下盖层均为可渗透的沉积物.本文利用新型的开采井设计方式,进行了单一...     

Gas hydrates (clathrates) causing pore-water freshening and oxygen isotope fractionation in deep-water sedimentary sections of terrigenous continental margins

The occurrence of gas hydrates in deep-water sections of the continental margins predicted from anomalous acoustic reflectors on seismic profiles has been confirmed by recent deep-sea drilling results. On the Pacific continental slope off Guatemala gas hydrates were brought up for the first time from two holes (497, 498A) drilled during Leg 67 of the DSDP in water depths of 2360 and 5500 m, respectively. The hydrates occur in organic matter-rich Pleistocene to Miocene terrigenous sediments. In the hydrate-bearing zone a marked decrease in interstitial water chlorinities was observed starting at about 10–20 m subbottom depth. Pore waters at the bottom of the holes (near 400 m subbottom) have as little as half the chlorinity of seawater (i.e. 9‰). Similar, but less pronounced, trends were observed during previous legs of the DSDP in other hydrate-prone segments of the continental margins where recharge of fresh water from the continent can be excluded (e.g. Leg 11). The crystallization of hydrates, like ice, excludes salt ions from the crystal structure. During burial the dissolved salts are separated from the solids. Subsidence results in a downward motion of the solids (including hydrates) relative to the pore fluids. Thawing of hydrates during recovery releases fresh water which is remixed with the pore fluid not involved in hydrate formation. The volume of the latter decreases downhole thus causing downward decreasing salinity (chlorinity). Hydrate formation is responsible for oxygen isotope fractionation with¹⁸O-enrichment in the hydrate explaining increasingly more positive δ¹⁸O values in the pore fluids recovered (after hydrate dissociation) with depth.     

时移地震监测天然气水合物开采可行性分析

时移地震储藏监测是一种提高油气采收率的有效技术,但该技术应用需要满足一定的储藏条件,因此时移地震的可行性研究非常有必要.天然气水合物的开采过程,主要是通过有关措施(加温、降压、注入化学材料等)使沉积物中的固态水合物变成气体开采出来,这将改变水合物储层的物理性质和地震反射特性,原理上为进行时移地震监测提供了条件.本文通过数值模拟分析发现:扩散型水合物藏在开采过程中有明显的地震信号变化,可以用时移地震对这类水合物藏的开采进行监测;渗漏型水合物藏的时移地震监测是否可行与其所含气体含量大小有关,气体含量低的可以用时移地震进行检测,气体含量高的是否能用时移地震进行检测要视具体情况而定.     

南海北部神狐海域天然气水合物形成聚集的数值模拟研究

利用2D数值方法对南海北部陆坡神狐海域水合物形成聚集过程进行了模拟,对气烟囱、泥底辟与水合物成藏间的关系进行研究.模拟结果表明,来自深部的甲烷热解气在向上运移过程中以垂向运动为主,且局限在某一狭窄的范围内,故在地震剖面上显示为气烟囱及顶部BSR.只有当其越过水合物稳定带底界,才能形成水合物,此时BSR等于水合物稳定带底...     

AVO study of a gas-hydrate deposit, offshore Costa Rica

Amplitude variation with offset (AVO) analysis and waveform inversion are techniques used to determine qualitative or quantitative information on gas hydrates and free gas in sediments. However, the quantitative contribution of gas hydrates to the acoustic impedance contrast observed at the bottom‐simulating reflector and the reliability of quantitative AVO analyses are still topics of discussion. In this study, common‐midpoint gathers from multichannel wide‐angle reflection seismic data, acquired offshore Costa Rica, have been processed to preserve true amplitude information at the bottom‐simulating reflector for a quantitative AVO analysis incorporating angles of incidence of up to 60°. Corrections were applied for effects that significantly alter the observed amplitudes, such as the source directivity. AVO and rock‐physics modelling indicate that free gas immediately beneath the gas‐hydrate stability zone can be detected and low concentrations can be quantified from AVO analysis, whereas the offset‐dependent reflectivity is not sensitive to gas‐hydrate concentrations of less than about 10% at the base of the gas‐hydrate stability zone. Bulk free‐gas saturations up to 5% have been determined from the reflection seismic data assuming a homogeneous distribution of free gas in the sediment. Assuming a patchy distribution of free gas increases the estimated concentrations up to 14%. There is a patchy occurrence of bottom‐simulating reflectors south‐east of the Nicoya Peninsula on the continental margin, offshore Costa Rica. AVO analysis indicates that this phenomenon is related to the local presence of free gas beneath the gas‐hydrate stability zone, probably related to a focused vertical fluid flow. In areas without bottom‐simulating reflectors, the results indicate that no free gas is present.     

南海北部陆坡块体搬运沉积体系的地震响应与成因机制

块体搬运沉积体系是一种重力作用下形成的深水沉积体系,在全球大陆边缘沉积体系中扮演着重要角色.块体搬运沉积作用不仅造成危害极大的深水地质灾害,而且与海洋天然气水合物和深水油气的形成过程与富集有密切联系.根据新获得的海底地形地貌资料和地震资料,揭示了块体搬运沉积体系识别特征,明确了其分布范围.这些块体搬运沉积体系具有明显的...     

青藏高原冻土带天然气水合物的形成条件与分布预测

冻土带是天然气水合物发育的两个重要地质环境之一.青藏高原平均海拔在4000m以上,多年冻土面积约1.4×106km2.本文根据青藏高原冻土层厚度和地温梯度特征,运用天然气水合物的热力学稳定域预测方法,确定中低纬度高海拔区冻土带天然气水合物的产出特征.青藏高原多年冻土带热成因天然气水合物形成的热力学相平衡反映,水合物顶界埋深约27~560m,底界埋深约77~2070m.初步计算表明,青藏高原冻土带水合物天然气资源约1.2×1011~2.4×1014m3.在冻土层越厚、冻土层及冻土层之下沉积层的地温梯度越小的地区,最有利于天然气水合物的发育.气温的季节性变化对天然气水合物影响不大.在全球气温快速上升的背景下,青藏高原天然气水合物将处于失稳状态,天然气水合物顶界下降、底界上升,与冻土带的退化相似,分布区逐渐缩小,最终将完全消失.     

Ultrasonic attenuation of pure tetrahydrofuran hydrate

Improved estimates of the amount of subsurface gas hydrates are needed for natural resource, geohazard, and climate impact assessments. To evaluate gas hydrate saturation from seismic methods, the properties of pure gas hydrates need to be known. Whereas the properties of sediments, specifically sands, and hydrate‐bearing sediments are well studied, the properties of pure hydrates are largely unknown. Hence, we present laboratory ultrasonic P‐wave velocity and attenuation measurements on pure tetrahydrofuran hydrates as they form with reducing temperatures from 25°C to 1°C under atmospheric pressure conditions. Tetrahydrofuran hydrates, with structure II symmetry, are considered as proxies for the structure I methane hydrates because both have similar effects on elastic properties of hydrate‐bearing sediments. We find that although velocity increased, the waveform frequency content and amplitude decreased after the hydrate formation reaction was complete, indicating an increase in P‐wave attenuation after hydrate formation. When the tetrahydrofuran hydrate was cooled below the freezing point of water, velocity and quality factor increased. Nuclear Magnetic Resonance results indicate the presence of water in the “pure hydrate” samples above the water freezing point, but none below. The presence of liquid water between hydrate grains most likely causes heightened attenuation in tetrahydrofuran hydrates above the freezing point of water. In naturally occurring hydrates, a similarly high attenuation might relate to the presence of water.     

南海北部陆坡水合物分解引起海底不稳定性的定量分析

分析了天然气水合物分解引起海底斜坡不稳定性的原因,探讨了水合物分解引起沉积物层孔隙压力变化规律.以南海北部神狐海域水合物样品钻获区为研究对象,应用极限平衡法探讨了海底无限斜坡的稳定性问题,计算了设定的不同滑动面的安全系数,建立了目标不稳定区相关参数条件下的水合物分解与海底斜坡不稳定的模式图,模拟计算了不同程度的水合物分解对海底斜坡失稳的影响作用.结果表明：在假设条件下,南海神狐海域当水深为1200 m,沉积物层厚度为200 m,对于20°的较大坡角,沉积物层中5%水合物分解将引起海底斜坡失稳;斜坡坡角为5°时,沉积物层中15%水合物分解后,将会引起海底斜坡失稳;当斜坡坡角为3°或更小时,沉积物层中25%水合物分解后,将会引起海底斜坡失稳.     

哈纳斯湖的水下地貌与表层沉积物的若干特征

测量表明,哈纳斯湖最大水深达188.5m,平均水深达120.1m,优质淡水容量达53.78×10~8m~3,均在我国湖泊中居突出地位。横断面为罕见的典型箱形,断陷作用仍在强烈进行之中。湖底表层沉积物粒度参数之分布主要受纵向搬运作用控制,局部受横向搬运作用影响。沉积物的矿物学特征与元素分布特征明显受沉积环境控制。沉积物中营养元素含量总体偏低。~(210)Pb法测得平均沉积速率为0.52mm/a,难以补偿断陷作用,湖盆将继续加深。     

海底天然气水合物层界面反射AVO数值模拟

本文采用AVO数值模拟方法,共选取水合物系统分层结构6个模型,对水合物、游离气和饱水沉积物接触界面的反射特征进行了数值模拟,研究了BSR及双BSR的存在条件与水合物体系垂向分布的关系,对一些现象从理论上进行了阐明.主要结论是：（1）强的似海底反射界面BSR与游离气体的存在密不可分,实际地震剖面中的“BSR”可能不对应水合物而只对应气体,无明显BSR的地方可能有水合物.（2）水合物顶部有可能存在游离气体,它可以在正常BSR（BSR1）之上形成另一个具有正极性的似海底反射界面BSR2.（3）正常BSR之下的双BSR（BSR0）其弹性机理有两种可能,一是水合物之下游离气上升迁移遇到某种阻隔层或不同气体组分的自然分层所造成的气体垂向分布的梯度性差异;二是当水合物之下的游离气体中存在残存的水合物时,同样会形成一定强度的似海底反射,在这种情形下BSR0的极性比较难以判断,取决于残存水合物上下游离气的饱和度和残存水合物的厚度.