Environmental reconstruction for the paleo-lake of ZhuⅠdepression and the depositional model for high-quality source rocks
-
摘要: 珠一坳陷古近纪陆相湖盆发育优质烃源岩，通过对文昌组、恩平组泥岩系统取样，进行常量、微量元素、同位素含量等分析，优选出对环境敏感的元素与元素比值，系统分析了珠一坳陷烃源岩形成时的古环境，包括古盐度、古水深、氧化还原条件、古气候及古水体的封闭性，并建立了有机碳与环境敏感元素比值的定量关系。开展优质烃源岩发育的水化学性质研究，明确研究区优质烃源岩形成演化所具有的古环境特征，最终探讨了文昌组优质烃源岩的发育模式。研究结果表明，珠一坳陷（除韩江凹陷无样本点）文昌期温暖略干气候条件下，湖盆水体主要为淡水—微咸水，发育中—深湖相，整体表现为开放—半封闭、还原—强还原的环境，极利于有机质的保存和优质烃源岩的发育，各凹陷古环境差异不大。烃源岩质量与咸度、水深、还原条件呈正相关，共同影响了优质烃源岩的发育。总之，温暖略干、微咸、深水及强还原、开放—半封闭环境共同控制了珠一坳陷优质烃源岩的发育。Abstract: High-quality source rocks are well occurred in the Paleogene lacustrine deposits of the ZhuⅠdepression. Element geochemistry and isotope analysis for mudstone cuttings from the Wenchang Formation and Enping Formation is carried out by this paper in order to optimize the study of element content and elements ratio which are sensitive to paleoenvironment reconstruction. Paleoenvironmental parameters for the source rocks in ZhuⅠdepression are systematically analyzed, which include paleosalinity, paleowater depth, redox conditions, paleoclimate, degree of environmental closure and so on. Upon the basis, we revealed the quantitative relationship between the organic matter content and the ratio of environmental sensitive elements, and established the depositional model for the high-quality source rocks in the Wenchang Formation. The results show that under the warm and slightly dry climate conditions during the Wenchang period, the ZhuⅠdepression was a medium to deep and open to semiclosed lake dominated by fresh water and brackish water under a reductive or strongly reductive environment, which was conducive to the preservation of organic matter and the formation of high-quality source rocks. The quality of source rocks is positively related to salinity, water depth and reduction conditions. In conclusion, in a warm and slightly dry lake with slightly salty deep water as the ZhuⅠdepression, the strong reduction and open-semiclosed environment jointly control the development of high-quality source rocks in ZhuⅠdepression.
-
随着大气中CO2浓度不断攀升，全球平均温度较1850—1900年间已升高近1.5℃，由此导致了极端气候频发、粮食产量减少、海平面上升等问题[1]。为了应对此类现象，全球各国陆续公布了碳中和目标。CO2捕集、利用和封存技术（Carbon Capture, Utilization and Storage简称CCUS），被认为是实现减排目标的关键技术[2-3]。其中，地质封存是实现CO2永久封存最有效的方式[4]。地质封存的主要场所是海底深部的咸水层或枯竭的油气藏，其中深部咸水层因其巨大的封存潜力以及不具备其他经济价值而被视为最理想的CO2地质封存场所[5-6]，据统计，咸水层CO2封存潜力占总封存潜力的98%以上[4]。国外目前正在运行的比较成功的咸水层封存项目包括挪威北海的Sleipner项目、Snøhvit项目，加拿大的Gorgon项目等，其中Sleipner作为世界首个海上咸水层CO2封存示范项目，以运行时间长、封存量大的特点著称，并部署了各种封存监测手段。自1996年起开始注入，至今已封存超过
2400 万吨CO2，是研究CO2咸水层封存的示范案例[7-8]。在国内，鄂尔多斯盆地的神华咸水层封存项目正在持续运行，中国海油建设的中国首个海上咸水层CO2封存项目恩平15-1油田也已正式开始CO2注入，但目前这些项目仍未部署CO2监测。针对CO2注入后的运移和泄露状况的动态监测问题[9]，地球物理被认为是一种有效的监测手段，用于确认地质封存的安全性、可靠性[10]。其中，与重力、测井、微地震等其他方法相比，四维地震在监测CO2体积和波及范围方面具有独特优势，经过先进方法[11-12]处理后能够提供随时间变化的地下动态[13]、高分辨率成像[14]。Huang等[15]对德国的咸水层CO2封存项目通过四维地震的振幅差异监测注入的CO2波及范围的变化，Chadwick和Noy[16]对Sleipner四维地震数据定性解释来匹配CO2羽流迁移的历史，Roach和White[17]使用了类似的方法来描述加拿大的Aquistore封存点的CO2羽流迁移，Fawad[18]为提高监测准确性利用四维地震和电磁测量相结合的方法来描绘CO2羽流。
为了建立CO2在海底咸水层封存过程中四维地震属性对CO2饱和度变化的定量描述关系，提高监测CO2封存波及范围的准确性。本研究结合Sleipner咸水层CO2封存项目所采集的四维地震监测资料和测井资料，通过Gassmann流体替换方程和Xu-White模型建立不同CO2饱和度下的速度模型，应用井控地震属性分析技术研究CO2咸水层封存过程中CO2-盐水两相介质变化引起的各向异性地震响应特征，提取并分析对CO2饱和度变化敏感的地震属性。同时，对Sleipner采集的7次时延地震数据进行解释，提取沿层地震属性，监测CO2封存波及范围的时空变化。
1. 研究区概况
1.1 Sleipner CO2封存项目概况
Sleipner是Equinor（原挪威国家石油公司）运营的一个重要天然气田，位于挪威北海。Sleipner CCS项目是其商业化的CCS站点，项目将Sleipner East气田的天然气生产过程中所产生的CO2压缩并通过约2.3 km的斜井注入高孔隙度Utsira咸水层，以防止CO2释放到大气中[19-21]（图1）。自1996年9月15日开始，Sleipner平台通过注入15/9-A-16井开始注入CO2，注入点位于海平面以下
1012 m的咸水层Utsira底部，距离储层顶部约200 m，注入射孔长度为38 m[8]。为了达到所需的井口压力，CO2需要先经过4个压缩阶段，在注入井井口CO2处于气液两相流动的相变阶段。最初几年的年注入率约为0.9 Mt，计划注入200年[22]。由于Sleipner East气田的天然气流量减少，后期略有下降，截至2010年，累计注入量和注入速度统计见图2。![]()
![]()
图 2 CO2日均注入体积及累计注入量Figure 2. CO2 Average daily injection volume and cumulative injection迄今为止，Sleipner共进行了10次三维地震调查和4次重力调查，通过这些调查，深入研究了储层中CO2的运移行为。当CO2进入储层时，通常处于超临界状态，与初始的咸水层形成了强烈的声速对比，为地震监测提供了有利条件。Sleipner的地震监测分别于1994年（Base line）、1999年、2001年、2002年、2004年、2006年、2008年、2010年、2013年和2016年进行了重复三维拖曳地震勘测，其中1994年的地震监测获得了注入CO2之前的基线三维地震数据[23-24]。
1.2 地质概况
Sleipner CO2封存地点位于Utsira砂层[25]，这是一个深度约800～
1100 m、厚度约200～300 m、面积约26000 km2的咸水层，南北长400多千米、东西长50～100 km[26]，主要由中新世晚期至上新世早期的砂质单元和部分泥质储层组成[27]（图3）。通过伽玛测井曲线以及强地震反射（图4）可以有效地识别Utsira储层中的砂泥岩。声速和密度分别在地层单元顶部和底部急剧下降和增加，构成了明显的阻抗差异，并引起强烈的清晰可识别的地震反射，根据15/9-13的测井曲线，对砂泥岩进行速度分析（图5），结果显示随着泥质含量的提高，声波速度整体呈上升的趋势，泥岩的速度整体要大于砂岩的速度。
![]()
图 4 15/9-13井GR测井曲线图（a）及2010年Inline 187地震剖面（b）数据来自https://CO2datashare.org。Figure 4. Well 15/9-13 GR logging profile (a), and the Inline 187 seismic profile in 2010 (b)Data from https://CO2datashare.org.根据测井曲线和三维地震资料，将目标地层分为9个主要的砂层（图4），Utsira储层之上的Nordland页岩层厚度约200～300 m，构成主要的盖层，从测井曲线可以看出，在8层和9层之间还有一层比较厚、横向分布比较广泛的泥岩，厚度约为5 m（图4a）。储层中的薄泥岩层平均厚度为1～1.5 m，构成了储层砂层内重要的渗透性屏障，并对CO2在储层内的运移和圈闭产生影响。储层单元的含砂量一般为0.7～1.0，15/9-A-16井获得的不同深度Utsira储层岩芯样品的主要岩石和碎屑成分分析结果见表1。
表 1 15/9-A-16井不同深度岩芯样本主要岩石和碎屑成分百分比Table 1. Percentage of major rock and debris compositions in core samples from Wells 15/9-16 at different depths碎屑含量/% 850～860 m 890～900 m 1000 ～1010 m石英 50.7 66.7 76.7 长石 7.3 3.7 2.7 方解石 18 17 7.7 页岩 4.3 1 4.7 根据挪威国家石油公司数据，Sleipner地区在
1038 m附近的地下温度测量值为37℃。地下温度梯度为33℃/km，因此，Sleipner地区Utsira储层的温度范围为储层顶部的28℃到储层底部的41℃[27]，具体温度、压力以及CO2密度关系如图6所示。![]()
根据表1矿物成分以及Sleipner地区的实际地质条件、储层温度和压力状况得到岩石基质矿物成分和储层内流体的相关模量和密度。其中在储层条件下，盐水的体积模量为2.3 GPa[28]，在储层温度和压力下，盐水的密度假设为1.03 g/cm3[29]。大部分CO2以超临界状态封存，考虑初始孔隙压力为8 MPa[30]，储层温度范围为28℃（Utsira上部地层）至41℃（注入点附近），选择CO2参考密度为700 kg/m3[31]。CO2的体积模量变化范围为0.02～0.075 GPa[32]，根据现场压力和温度条件，选择使用0.075 GPa的CO2体积模量[31]，具体参量统计见表2。
表 2 储层中矿物成分和流体的弹性模量及密度Table 2. Elastic modulus and density of mineral components and fluids in reservoirs体积模量/GPa 剪切模量/GPa 密度/(g/cm3) 石英 37.00 44.00 2.65 长石 37.50 15.00 2.70 方解石 76.80 32.00 2.71 盐水 2.30 0 1.03 CO2 0.075 0 0.70 Utsira储层主要是高孔隙度砂岩，通过岩芯样品的液体侵入测量和密度测井得出的孔隙度值范围为36.0%～40.1%，平均38.0%[27]。从15/9-A-16井的中子孔隙度测井曲线得到孔隙度为35%～39%，平均37%。Williams[33]对渗透率和电缆测井数据的重新评估认为东西方向渗透率约为1.975 μm2，南北方向渗透率高于东西方向，约为7.896 μm2。
从地震剖面可以看出，CO2通过几个相对高渗透性“烟囱”垂直向上运移，穿过储层及其地层内的泥岩[21]（图7）。“烟囱1”和“烟囱3”对应于地层中微妙的不连续面，位于或接近储层顶部（图7a）。“烟囱1”在重复测量中特别突出（图7b），并且可以作为一个独特的垂直圆形反射，半径约为30 m，它被解释为CO2通过羽流向上运移的主要垂直管道[34]。
![]()
图 7 四维地震剖面对比显示3个可能的CO2羽流烟囱a：1994年基准地震剖面，b：2010年时移地震剖面。Figure 7. Comparison of 4D seismic profiles showing three possible CO2 plume chimneysa: The 1994 reference seismic profile, b: 2010 time-shifted seismic profile.2. CO2封存地震响应特征正演模拟
为了研究CO2封存过程中随CO2饱和度变化的四维地震响应特征，需通过测井数据建立不同CO2饱和度的正演模型，具体包括以下3个方面的内容：① 采用流体替换方法预测不同CO2饱和度条件下的纵波和横波速度；② 基于流体替换前后的速度和密度曲线构建水平层状模型，正演模拟得到自激自收地震记录；③ 通过对不同CO2饱和度下的地震记录提取相应的地震属性，分析对于CO2饱和度变化敏感的属性以用于地震资料解释和CO2波及范围监测。
2.1 不同CO2饱和度地层的速度模型
在进行CO2封存的过程中，CO2注入导致储层内原有流体被置换，同时岩石框架和饱和岩石体积模量等岩石物理参数发生改变[35]。通过分析CO2注入期间地震资料的变化，能够监测到不同阶段CO2分布状况。因此需要建立CO2饱和度与地震速度之间的关系，即对注入CO2后的地震速度进行估算。
Sleipner咸水层CO2封存项目的储层为砂泥岩互层，因此本文采用Xu-White模型[36]求取纵横波速度。该模型是在Gassmann流体替换方程的基础之上，综合考虑了黏土影响和孔隙纵横比，结合Kuster-Toksoz模型[37]和差分有效介质（DEM）理论，提出的一种依据孔隙度和泥质含量来估算泥质砂岩纵波与横波速度的方法。其中纵横波速度计算公式如下：
$$ \begin{array}{c}{V}_{\mathrm{P}}=\sqrt{\frac{{K}_{\mathrm{s}\mathrm{a}\mathrm{t}}+\dfrac{4}{3}{\mu }_{\mathrm{s}\mathrm{a}\mathrm{t}}}{\rho }}\end{array} $$ (1) $$ {V}_{\mathrm{S}}=\sqrt{\frac{{\mu }_{\mathrm{s}\mathrm{a}\mathrm{t}}}{\rho }} $$ (2) 式（1）、（2）中，VS为横波速度，VP为纵波速度，Ksat为饱和岩石的体积模量，μsat为饱和岩石的剪切模量，ρ为饱和岩石密度。
Xu-White速度模型的公式如下：
$$ \phi ={\phi }_{\mathrm{s}}+{\phi }_{\mathrm{c}} $$ (3) $$ {\phi }_{\mathrm{s}}={f}_{\mathrm{s}}\frac{\phi }{1-\phi } $$ (4) $$ {\phi }_{\mathrm{c}}={f}_{\mathrm{c}}\frac{\phi }{1-\phi } $$ (5) $$\begin{split} K_{\mathrm{d}\mathrm{r}\mathrm{y}}-K_{\mathrm{m}\mathrm{a}}=\;&\frac{1}{3}\left(K_{\mathrm{f}}-K_{\mathrm{m}\mathrm{a}}\right)\frac{3K_{\mathrm{d}\mathrm{r}\mathrm{y}}+4K_{\mathrm{m}\mathrm{a}}}{3K_{\mathrm{m}\mathrm{a}}+4\mu_{\mathrm{m}\mathrm{a}}}\times\\&\sum_{l=s,c}^{ }\phi_lT_{iijj}\left(\alpha_l\right) \end{split}$$ (6) $$ \begin{split} & \mu_{\mathrm{d}\mathrm{r}\mathrm{y}}-\mu_{\mathrm{m}\mathrm{a}}=\frac{\mu_{\mathrm{f}}-\mu_{\mathrm{m}\mathrm{a}}}{5}\times \\ &\quad\frac{6\mu_{\mathrm{d}\mathrm{r}\mathrm{y}}\left(K_{\mathrm{m}\mathrm{a}}+2\mu_{\mathrm{m}\mathrm{a}}\right)+\mu_{\mathrm{m}\mathrm{a}}\left(9K_{\mathrm{m}\mathrm{a}}+8\mu_{\mathrm{m}\mathrm{a}}\right)}{5\mu_{\mathrm{m}\mathrm{a}}\left(3K_{\mathrm{m}\mathrm{a}}+4\mu_{\mathrm{m}\mathrm{a}}\right)} \times\\ &\quad\sum_{l=s,c}^{ }\phi_lF\left(\mathrm{\alpha}_l\right) \end{split} $$ (7) $$ F\left(\mathrm{\alpha }\right)={T}_{ijij}\left(\mathrm{\alpha }\right)-\frac{{T}_{iijj}\left(\mathrm{\alpha }\right)}{3} $$ (8) 式中，ϕ为总孔隙度；ϕs为砂岩孔隙度；ϕc为泥岩孔隙度；fs为砂的体积分数；fc为黏土的体积分数；Kdry、Kma、Kf分别为干岩石体积模量、组成岩石的固体矿物基质的体积模量以及孔隙流体体积模量。μdry、μma、μf分别为对应的剪切模量。F(α)、Tijij(α)为关于孔隙纵横比α的函数。公式（6）和（7）为Kuster-Toksoz方程。
Xu-White模型将岩石视为砂质和泥质的复合体。其假定岩石骨架的泊松比大致为恒定值，强调了孔隙形态对地震波速的显著影响。通过引入孔隙的纵横比以刻画孔隙的几何特性，并将其与孔隙率、岩石基质的剪切及体积模量相结合，从而推导出干岩石体积模量Kdry与剪切模量μdry的简化表达式：
$$ \begin{array}{c}{K}_{\mathrm{d}\mathrm{r}\mathrm{y}}={K}_{\mathrm{m}\mathrm{a}}(1-\phi {)}^{p}\end{array} $$ (9) $$ {\mu }_{\mathrm{d}\mathrm{r}\mathrm{y}}={\mu }_{\mathrm{m}\mathrm{a}}(1-\phi {)}^{q} $$ (10) 式中，参数$ p $与$ q $的表达式分别为：
$$ p=\frac{1}{3}\mathop\sum\nolimits _{l=\mathrm{1,2}}{f}_{l}{T}_{ij}\left({\mathrm{\alpha }}_{l}\right) $$ (11) $$ q=\frac{1}{5}\mathop\sum\nolimits _{l=\mathrm{1,2}} {f}_{l}{F}_{ij}\left({\mathrm{\alpha }}_{l}\right) $$ (12) 同时，引入了Voigt-Reuss-Hill（VRH）[38-39]平均理论计算岩石基质模量，公式如下:
$$ {M}_{\mathrm{V}\mathrm{R}\mathrm{H}}=\frac{{M}_{\mathrm{\nu }}+{M}_{\mathrm{R}}}{2} $$ (13) 式中，MVRH为岩石基质的弹性模量；Mν为Voigt有效弹性模量；MR为Reuss等效弹性模量。
引入Brie[40]经验模型计算不同CO2饱和度下的混合流体模量，公式为：
$$ {K}_{\mathrm{f}}=\left({K}_{\mathrm{w}}-{K}_{\mathrm{g}}\right)\cdot {S}_{\mathrm{w}}^{e}+{K}_{\mathrm{g}} $$ (14) 式中，Kf为混合流体的体积模量，Kw为液相体积模量，Kg为气相体积模量，Sw为液相饱和度，e为经验指数因子。
基于上述模型编程计算得到储层在不同CO2饱和度状态下的速度曲线，图8分别对比了从0开始注入CO2，在储层CO2饱和度为0、20%、40%、60%、80%时的流体替换反射纵横波速度曲线。可以看出，随着储层CO2饱和度的提高，纵横波速度均持续下降，这是由于随CO2饱和度的增加，Xu-White模型计算得到的岩石体积模量Ksat和剪切模量μsat均减小，因此流体替换得到的速度减小。
![]()
图 8 不同CO2饱和度的纵波（a）及横波（b）速度曲线对比Figure 8. Comparison of primary wave velocity profiles (a) and shear wave velocity profiles (b) for different CO2 saturations图9对比了不同饱和度下的均方根速度，可以看出，随着CO2饱和度的升高，纵横波均方根均下降，注入CO2后到CO2饱和度为60%之前，横波速度下降幅度略小于纵波速度下降幅度，并且幅度随饱和度增加而减小，当CO2饱和度大于60%之后，横波速度减小的幅度将超过纵波速度。
![]()
图 9 随CO2饱和度升高纵横波均方根速度变化趋势Figure 9. Trend of root-mean-square velocity of longitudinal and transverse waves with increasing CO2 saturation2.2 地震正演模拟与敏感属性分析
基于不同CO2饱和度替换得到的速度和密度曲线建立水平层状模型。对其进行自激自收的地震正演模拟得到储层CO2饱和度分别为0、20%、40%、60%、80%的地震记录，如图10所示。对比发现，振幅变化最为明显，随CO2饱和度的增加，盖层振幅保持不变，储层振幅提高。在储盖层交界处，相位发生变化。
在不同饱和度的地震记录中沿剖面提取振幅、频率、相位、吸收衰减地震属性。
（1） 均方根振幅属性
其中均方根（RMS）振幅属性是计算振幅平方和的平方根，并将其除以特定时间窗内的样本数量。使用RMS振幅属性是为了突出显示异常的高振幅，包括负值和正值。RMS振幅属性是一个很好的声学对比指标，可以指示岩性和油气的变化[41-42]，计算公式如下：
$$ {A}_{\mathrm{R}\mathrm{M}\mathrm{S}}=\sqrt{\frac{1}{\mathrm{N}}\mathop\sum\nolimits _{i=1}^{\mathrm{N}} {a}_{i}^{2}} $$ (15) 式中，N为采样点总数，ai为各采样点振幅值。
图11a对比了不同CO2饱和度的均方根振幅属性，随着饱和度的增加，盖层均方根振幅值不变，储层均方根振幅值变大。通过将不同饱和度的均方根振幅属性做差（图11b）可以看出，开始注入时，振幅变化最明显，随着CO2饱和度增加，振幅增大的幅度减弱，饱和度越大振幅变化越不明显，但仍能看出变化趋势。
![]()
图 11 不同CO2饱和度均方根振幅属性（a）及均方根振幅属性差（b）Figure 11. RMS amplitude attributes (a) and difference in RMS amplitude attributes (b) for different CO2 saturations（2） 瞬时频率属性
瞬时频率作为相位随时间变化率的量度，能够反映地层的岩性变化。受地层影响，其变化幅度较大，如当地震波穿过某些特殊的储层时，高频成分将会显著衰减，其中含油气储层经常引起高频衰减[41]。因此瞬时频率属性可反映地震波衰减特征，以此估计地层的周期性，其计算公式如下：
$$ \omega \left(t\right)=\frac{\mathrm{d}\phi \left(t\right)}{\mathrm{d}t}$$ (16) 图12a为不同CO2饱和度的瞬时频率属性，当注入CO2之后，储盖层的分界面处以及上下两侧频率发生变化，将不同饱和度的瞬时频率属性做差（图12b）可以发现，在储盖层交界处瞬时频率升高，在上下两侧瞬时频率降低，但随着CO2饱和度的增大，瞬时频率属性的变化敏感性减弱。饱和度60%以上，频率变化几乎为0。
![]()
图 12 不同CO2饱和度瞬时频率属性（a）及瞬时频率属性差（b）Figure 12. Instantaneous frequency attributes (a) and difference in instantaneous frequency attributes (a) for different CO2 saturations（3） 瞬时相位属性
瞬时相位属性一般是用来反映同相轴的连续性，不考虑振幅强度，地震波穿越不同的地层会引起地震波的相位变化，因此通常是用来识别断层、尖灭点、河道等[43]，计算公式如下：
$$ \phi \left(t\right)={\mathrm{tan}}^{-1}\left[\frac{g\left(t\right)}{f\left(t\right)}\right] $$ (17) 图13a对比了不同饱和度的相位属性，CO2注入前后几乎看不出变化。将不同饱和度的瞬时相位属性做差可以发现（图13b），当开始注入CO2时，在储层和盖层交界处出现相位变化，当CO2饱和度增大到40%以上时，瞬时相位属性基本未发生改变。
![]()
图 13 不同CO2饱和度瞬时相位属性（a）及 瞬时相位属性差（b）Figure 13. Instantaneous phase attributes (a) and the difference in instantaneous phase attributes (b) for different CO2 saturations（4） 瞬时Q值属性
瞬时Q值属性也叫做瞬时品质因子，是一种吸收衰减属性，可通过低频振幅与高频振幅的比值得到。其主要用于描述地震波在介质中传播时的衰减特性。瞬时Q值越大，说明低频能量与高频能量的比值越大，高频的衰减越剧烈。流体饱和度的变化会影响介质的衰减特性，因此可以用来评估地层中的流体饱和度。
图14a为不同饱和度CO2的瞬时Q值属性图。当CO2注入后，在储盖层交界处周围发生了较明显的衰减变化。将不同饱和度的瞬时Q值属性做差可以发现（图14b），注入CO2后，在储层周围出现正负相间的衰减变化，随CO2饱和度增长，衰减变化值变小，但仍具有敏感性。
![]()
图 14 不同CO2饱和度瞬时Q值属性（a）及瞬时相位属性差（b）Figure 14. Instantaneous Q attributes (a) and the difference in instantaneous phase attributes (b) for different CO2 saturations3. CO2咸水层封存波及范围地震监测
为了研究CO2羽流在平面以及垂向上的发育情况，沿着泥岩层反射分别解释了上（层9）、中（层5）、下(层1)3个层位(图15)，网格间距为16×16，解释结果如图16所示。根据正演结果选择均方根振幅属性来描述CO2波及范围。沿解释层位，创建时窗并提取均方根振幅属性。图17显示，以中间层第5层1999年、2006年、2010年3次时延地震资料为例，图中绿色和黄色为振幅高值，蓝和紫色为振幅低值，通过同一层位几次不同时间的属性图对比可以明显看出整体CO2波及范围随着时间增长，此外，随CO2饱和度增加，均方根振幅值升高。周围储层的均方根振幅值为0.1左右，CO2饱和度最高值处振幅在2左右，CO2羽流波及前缘处振幅值大于0.5，对波及范围边界的识别效果较好，同时与地震正演结果相符。通过均方根振幅属性开展不同层位、不同时间的CO2波及范围对比，分上（第1）、中（第5）、下（第9层）3个层位提取2010—1999年6次四维地震数据的均方根振幅属性（图18），对比不同时间，上、中、下三层平面上CO2羽流的纵向扩散速度均大于横向扩散速度，整体均呈SSW-NNE方向展布。
![]()
图 17 Sleipner咸水层封存第5层四维地震均方根振幅属性Figure 17. 4-D seismic RMS amplitude attributes for Layer 5 of the Sleipner saline aquifer storage![]()
图 18 挪威Sleipner咸水层封存项目CO2平面波及范围四维地震属性预测结果Figure 18. Predicted 4D seismic attributes of CO2 planar spread range for the Sleipner saline aquifer storage project in Norway图18显示了位于注入点上方的第1层平面上的CO2波及范围，通过其随时间的变化可以看出，开始注入时CO2呈一个高饱和度的团状羽流，随着时间的推移，CO2持续注入，CO2在平面上不断扩散，并在2004年达到最大横向饱和度，这时储层中注入的CO2量为6.9 Mt，随后沿纵向继续扩散，整体纵向扩散大于横向，随CO2羽流扩散，中心的饱和度不断减小，表明CO2正在不断向上层砂岩运移。
图18展示了中部第5层的CO2波及范围变化情况，在2008年时第5层达到了最大横向扩散，此时储层中注入的CO2量为10.2 Mt，在2010年注入12 Mt后整体饱和度达到最大，中心振幅小于两端，表明CO2在波及范围两端的构造高点处形成高饱和度积聚，中心通过高渗透率“烟囱”持续向上运移。对比发现，1999年时第1层和第5层的CO2波及范围大小相近，随后CO2波及范围在第5层的增长速度大于第1层的增长速度，最终的CO2波及范围第1层也大于第5层。
图18同时显示了储层最顶部砂岩层第9层CO2波及范围变化情况，可以看出1999年已经有CO2羽流到达最顶层，这进一步证明了存在CO2高渗透率“烟囱”。2010年CO2波及范围达到最大，除主体羽流发育以外，CO2沿着右侧一条通道向北快速运移，到达一个构造高点后积聚，继续向不同的方向扩散，整体上东侧振幅高于西测，通道处振幅值最高，CO2到达第9层后，先是平面扩散，随后沿通道运移，并在通道中聚集使饱和度提高。第9层整体上CO2波及范围小于第5层而大于第1层。
4. 结论与建议
本文基于Sleipner咸水层CO2封存项目的测井资料、四维地震资料，开展了针对CO2海底咸水层封存的流体替换和地震正演模拟，以及CO2波及范围地震监测，得到如下认识：
（1）均方根振幅属性对CO2饱和度变化最为敏感。随CO2饱和度的增加，饱和岩石的体积模量、体积密度、纵波速度和横波速度均有所下降。相反，地震记录中总体振幅升高，且随着CO2注入量的增加，振幅变化幅度减小。
（2）利用均方根振幅属性对Sleipner海底咸水层CO2波及范围进行描述，在注入期间，CO2在层内主要沿SSW-NNE方向运移，并在构造高部位聚集。垂向上，CO2从注入点向上层运移，下层达到最大波及范围的时间早于上层。结合构造解释结果和储层物性发现，CO2在储层内的波及范围主要受各项异性渗透率和构造高低控制。由于井15/9-A-16缺少实测速度值，本文在速度预测中使用的是固定的孔隙纵横比的Xu-White模型，在今后的研究中，可以采用由实测速度进行矫正的岩石物理模型，提高速度预测精度。CO2在咸水层封存过程中会改变储层的孔隙度等储层参数，在今后岩石物理建模过程中，可以将CO2对储层孔渗条件的影响考虑在内，建立动态岩石物理模型。
-
![]()
图 1 珠一坳陷地质背景纲要图
a. 珠一坳陷构造区划图，b. 珠一坳陷地层柱状图，c. 珠一坳陷纹层状烃源岩背散射图，d. 珠一坳陷分散状烃源岩背散射图。黄色箭头所指为有机碳。
Figure 1. Outline of geological background of ZhuⅠdepression
a. Tectonic map， b. stratigraphic column，c. laminated source rock，d. dispersed source rock. Note: the yellow arrow indicates organic carbon.
![]()
图 2 珠一坳陷文昌期湖盆古盐度恢复参数（Rb/K，Sr/Ba）箱线图
Figure 2. Box diagram of paleosalinity restoration for Wenchang Formation in ZhuⅠdepression using element parameters Rb/K and Sr/Ba
![]()
图 3 珠一坳陷文昌期湖盆古水深恢复参数（Mn/Fe，m值）箱线图
Figure 3. Box diagram of ancient water depth restoration for Wenchang Formation in ZhuⅠdepression using element parameters Mn/Fe and m value
![]()
图 4 珠一坳陷文昌期湖盆古氧相恢复、元素参数（Th/U，Cu/Zn）箱线图及烃源岩镜下小粒径黄铁矿特征
Figure 4. Box diagram of paleooxygen facies restoration for Wenchang Formation in ZhuⅠdepression using element parameters Th/U and Cu/Zn, and microscopic characteristics of fine pyrite in source rocks
![]()
图 5 珠一坳陷文昌期湖盆古气候恢复参数（Sr/Cu，Rb/Sr）箱线图
Figure 5. Box diagram of paleoclimate restoration for Wenchang Formation in ZhuⅠdepression using element parameters Sr/Cu and Rb/Sr
![]()
图 6 珠一坳陷文昌组、恩平组及东营凹陷沙河街组（数据引自文献33-34）无机碳、氧同位素δ值交汇图
Figure 6. Intersection diagram of inorganic Carbon and oxygen isotopes δ Value for Wenchang Formation, Enping Formation in ZhuⅠdepression and Shahejie Formation in Dongying depression (data from reference 33 and 34)
![]()
图 7 珠一坳陷文昌组地层特征及凝灰质泥岩显微照片
Figure 7. Geochemical element stratigraphic column and micrograph of tuffaceous mudstone for Wenchang Formation in ZhuⅠdepression
![]()
图 8 珠一坳陷古环境与烃源岩质量相关性分析
Figure 8. Correlation analysis between paleoenvironment and source rock quality in ZhuⅠdepression
![]()
图 9 珠一坳陷优质烃源岩形成环境模式图
Figure 9. Environment model for High-quality source rocks formation in Zhu Ⅰ depression
表 1 珠一坳陷文昌期湖盆古盐度重建指标及其指示意义
Table 1 Reconstruction index of paleosalinity for Wenchang Formation in ZhuⅠdepression and its implications
判别指标 淡水 微—半咸水 咸水 B/Ga ＜1.5 5～6 ＞7 Rb/K ＜0.004 0.004～0.006 ＞0.006 Sr/Ba ＜0.6为陆相，＞1为海相 Z值 ＜120 − ＞120 
下载: 导出CSV
表 2 珠一坳陷文昌期湖盆古氧相重建指标及其指示意义
Table 2 Reconstruction index of paleooxygen facies for Wenchang Formation in ZhuⅠdepression and its indicative significance
判别指标 氧化环境 弱还原—氧化 还原环境 强还原环境 Th/U ＞30 10～30 4～10 ＜4 V/(V+Ni) ＜0.46 0.46～0.60 ＞0.60 − Cu/Zn ＞0.63 0.38～0.63 0.21～0.38 ＜0.21
下载: 导出CSV
表 3 珠一坳陷文昌期湖盆古气候重建指标及其指示意义
Table 3 Reconstruction index of paleoclimate in ZhuⅠdepression for Wenchang Formation and its indicative significance
判别指标 低值 高值 Sr/Cu 1.3～5.0，温湿 ＞5，干热 Sr/Ca 温湿 干热 Mg/Ca 潮湿 干旱 Rb/Sr 干燥 湿润 Al2O3/MgO 干旱 湿润
下载: 导出CSV
-
[1] 朱伟林. 中国近海新生代含油气盆地古湖泊学与烃源条件[M]. 北京: 地质出版社, 2009. ZHU Weilin. Paleolimnology and Source Rock Studies of Cenozoic Hydrocarbon-bearing Offshore Basins in China[M]. Beijing: Geological Publishing House, 2009.
[2] Mello M R, Maxwell J R. Organic geochemical and biological marker characterization of source rocks and oils derived from lacustrine environments in the Brazilian continental margin[C]//Lacustrine Basin Exploration: Case Studies and Modern Analogs. Tulsa: AAPG, 1999, 50: 77-97.
[3] Rosendahl B R, Reynolds D J, Lorber P M, et al. Structural expressions of rifting: lessons from Lake Tanganyika, Africa [J]. Geological Society, London, Special Publications, 1986, 25(1): 29-43. doi: 10.1144/GSL.SP.1986.025.01.04
[4] 孟庆涛, 刘招君, 胡菲, 等. 桦甸盆地始新世古湖泊生产力与有机质富集机制[J]. 中国石油大学学报:自然科学版, 2012, 36(5):38-44 MENG Qingtao, LIU Zhaojun, HU Fei, et al. Productivity of Eocene ancient lake and enrichment mechanism of organic matter in Huadian Basin [J]. Journal of China University of Petroleum, 2012, 36(5): 38-44.
[5] 吴克强, 姜雪. 孙和风. 近海富生油凹陷湖相烃源岩发育模式: 以黄河口凹陷古近系为例[J]. 地质科技情报, 2015, 34(2):63-70 WU Keqiang, JIANG Xue, SUN Hefeng. Model of lacustrine source rocks in offshore oil kitchen sag: a case study of Paleogene in Huanghekou sag [J]. Geological Science and Technology Information, 2015, 34(2): 63-70.
[6] 马小祥, 姚素平, 张柏林, 等. 渤海湾盆地东濮凹陷古近系古湖盆氧化还原条件及其优质烃源岩的发育模式[J]. 高校地质学报, 2019, 25(6):801-812 MA Xiaxiang, YAO Suping, ZHANG Bolin, et al. Redox conditions of Paleogene paleolake and development models of high-quality source rocks in the Dongpu sag, Bohai Bay Basin [J]. Geological Journal of China Universities, 2019, 25(6): 801-812.
[7] 王书荣, 宋到福. 何登发. 三塘湖盆地火山灰对沉积有机质的富集效应及凝灰质烃源岩发育模式[J]. 石油学报, 2013, 34(6):1077-1087 WANG Shurong, SONG Daofu, HE Dengfa. The enrichment effect of organic materials by volcanic ash in sediments of the Santanghu Basin and the evolutionary pattern of tuffaceous source rocks [J]. Acta Petrolei Sinica, 2013, 34(6): 1077-1087.
[8] 施和生, 舒誉, 杜家元, 等. 珠江口盆地古近系石油地质[M]. 北京: 地质出版社, 2017 SHI Hesheng, SHU Yu, DU Jiayuan, et al. Petroleum Geology of Paleogene in Pearl River Mouth Basin[M]. Beijing: Geological Publishing House, 2017.
[9] 刘海伦. 珠江口盆地珠一坳陷裂陷结构: 基底属性与区域应力联合制约[D]. 中国地质大学博士学位论文, 2018 LIU Hailun. Rift style controlled by basement attribute and regional stress in ZhuⅠ depression, Pearl River Mouth basin[D]. Doctor Dissertation of China University of Geosciences, 2018.
[10] 吴克强, 刘志峰, 王升兰, 等. 珠一坳陷北部洼陷带始新统半深—深湖相烃源岩综合判识[J]. 中国海上油气, 2015, 27(3):10-15, 24 WU Keqiang, LIU Zhifeng, WANG Shenglan, et al. Composite recognition of eocene semi-deep and deep lacustrine facies source rocks in northern subsags belt of ZhuⅠ depression, Pearl River Mouth basin [J]. China Offshore Oil and Gas, 2015, 27(3): 10-15, 24.
[11] 高阳东, 林鹤鸣, 汪旭东, 等. 幕式裂陷控洼背景下的烃源岩分布及岩浆改造——以珠一坳陷番禺4洼为例[J]. 2021, 41(4): 151-160. GAO Yangdong, LIN Heming, WANG Xudong, et al. Source rock distribution pattern in an episodic rifting sag and later stage magmatiic reformation: A case from Panyu 4 sag, ZhuⅠDepression[J]. 2021, 41(4): 151-160.
[12] 朱明, 张向涛, 黄玉平, 等. 珠江口盆地烃源岩特征及资源潜力[J]. 石油学报, 2019, 40(S1):53-68 ZHU Ming, ZHANG Xiangtao, HUANG Yuping, et al. Source rock characteristics and resource potential in Pearl River Mouth Basin [J]. Acta Petrolei Sinica, 2019, 40(S1): 53-68.
[13] 毛光周, 刘晓通, 安鹏瑞, 等. 无机地球化学指标在古盐度恢复中的应用及展望[J]. 山东科技大学学报:自然科学版, 2018, 37(1):92-102, 118 MAO Guangzhou, LIU Xiaotong, AN Pengrui, et al. Application and outlook of inorganic geochemical indexes in reconstruction of palaeosalinity [J]. Journal of Shandong University of Science and Technology:Natural Science, 2018, 37(1): 92-102, 118.
[14] 王益友, 郭文莹, 张国栋. 几种地球化学标志在金湖凹陷阜宁群沉积环境中的应用[J]. 同济大学学报:自然科学版, 1979(2):51-60 WANG Yiyou, GUO Wenyin, ZHANG Guodong. Application of some geochemical indications in determining of sedimentary environment of the Funing Group (Paleogene), Jin-Hu Depression, Kiangsu Province [J]. Journal of Tongji University:Natural Science, 1979(2): 51-60.
[15] 李进龙, 陈东敬. 古盐度定量研究方法综述[J]. 油气地质与采收率, 2003, 10(5):1-3 LI Jinlong, CHEN Dongjing. Summary of quantified research method on paleosalinity [J]. Petroleum Geology and Recovery Efficiency, 2003, 10(5): 1-3.
[16] 张金亮, 张鑫. 塔里木盆地志留系古海洋沉积环境的元素地球化学特征[J]. 中国海洋大学学报, 2006, 36(2):200-208 ZHANG Jinliang, ZHANG Xin. The element geochemical features of ancient oceanic sedimentary environments in the Silurian period in the Tarim Basin [J]. Periodical of Ocean University of China, 2006, 36(2): 200-208.
[17] Keith M L, Weber J N. Carbon and oxygen isotopic composition of selected limestones and fossils [J]. Geochimica et Cosmochimica Acta, 1964, 28(10-11): 1787-1816. doi: 10.1016/0016-7037(64)90022-5
[18] 庞军刚, 杨友运, 郝磊. 湖盆古水深恢复研究现状综述[J]. 长江大学学报:自然科学版 理工, 2012, 9(9):42-45 PANG Jungang, YANG Youyun, HAO Lei. Review on the research status of ancient water depth restoration in the lake basin [J]. Journal of Yangtze University:Natural Science Edition, Science and Engineering, 2012, 9(9): 42-45.
[19] 逄淑伊, 操应长, 梁超. 渤海湾盆地东营凹陷沙四上亚段—沙三下亚段岩相特征及沉积环境——以樊页1井为例[J]. 石油与天然气地质, 2019, 40(4):799-809 PANG Shuyi, CAO Yingchang, LIANG Chao. Lithofacies characteristics and sedimentary environment of Es4U and Es3L: A case study of Well FY1 in Dongying sag, Bohai Bay Basin [J]. Oil & Gas Geology, 2019, 40(4): 799-809.
[20] 许中杰, 程日辉, 张莉, 等. 华南陆缘晚三叠-早、中侏罗世海平面相对升降与古气候演化的地球化学记录[J]. 地球科学:中国地质大学学报, 2012, 37(1):113-124 XU Zhongjie, CHENG Rihui, ZHANG Li, et al. The geochemistry records of sea-level relative movement and paleoclimatic evolution of the South China continental margin in Late Triassic-Early-Middle Jurassic [J]. Earth Science:Journal of China University of Geosciences, 2012, 37(1): 113-124.
[21] 李明龙, 陈林, 田景春, 等. 鄂西走马地区南华纪古城期-南沱早期古气候和古氧相演化: 来自细碎屑岩元素地球化学的证据[J]. 地质学报, 2019, 93(9):2158-2170 doi: 10.3969/j.issn.0001-5717.2019.09.004 LI Minglong, CHEN Lin, TIAN Jingchun, et al. Paleoclimate and paleo-oxygen evolution during the Gucheng Period-early Nantuo Period of Nanhua system in the Zouma area, West Hubei: Evidence from elemental geochemistry of fine Clastic rocks [J]. Acta Geologica Sinica, 2019, 93(9): 2158-2170. doi: 10.3969/j.issn.0001-5717.2019.09.004
[22] Jones G E. An ecological survey of open ocean and estuarine microbial populations. 1. The importance of trace metal ions to microorganisms in the sea[M]//Stevenson L H, Colwell R R. Estuarine Microbial Ecology. Columbia: University of South Carolina Press, 1973: 233-241.
[23] Emerson S R, Huested S S. Ocean anoxia and the concentrations of molybdenum and vanadium in seawater [J]. Marine Chemistry, 1991, 34(3-4): 177-196. doi: 10.1016/0304-4203(91)90002-E
[24] 梅水泉. 岩石化学在湖南前震旦系沉积环境及铀来源研究中的应用[J]. 湖南地质, 1988, 7(3):25-31, 49 MEI Shuiquan. Application of rock chemistry in the study of presinian sedimentary environment and the source of uranium mineralization in Hunan province [J]. Hunan Geology, 1988, 7(3): 25-31, 49.
[25] 田涛, 付德亮, 周世新, 等. 米仓山-汉南隆起区牛蹄塘组页岩古氧相及其与有机质富集的关系[J]. 兰州大学学报:自然科学版, 2020, 56(1):37-47,55 TIAN Tao, FU Deliang, ZHOU Shixin, et al. The paleo-redox conditions of the shale in Niutitang formation and its effects on organic matter enrichment of the Micangshan-Hannan Uplift [J]. Journal of Lanzhou University:Natural Sciences, 2020, 56(1): 37-47,55.
[26] 刘莹, 展翅飞. 探究影响Pr/Ph值的地质因素及其生烃演化[J]. 中国石油和化工标准与质量, 2014(5):203-203 LIU Ying, ZHAN Chifei. Explore the geological gactors affecting Pr/Ph and its hydrocarbon generation evolution [J]. China Petroleum and Chemical Standard and Quality, 2014(5): 203-203.
[27] Lerman A, Imboden D M, Gat J R. Physics and Chemistry of Lakes[M]. 2nd ed. Berlin, Heidelberg: Springer, 1995.
[28] 石军, 邹艳荣, 余江, 等. 准噶尔盆地阜康凹陷芦草沟组高有机碳页岩发育的古环境[J]. 天然气地球科学, 2018, 29(8):1138-1150 doi: 10.11764/j.issn.1672-1926.2018.06.016 SHI Jun, ZOU Yanrong, YU Jiang, et al. Paleoenvironment of organic-rich shale from the Lucaogou Fromation in the Fukang Sag, Junggar Basin, China [J]. Natural Gas Geoscience, 2018, 29(8): 1138-1150. doi: 10.11764/j.issn.1672-1926.2018.06.016
[29] 戴贤铎, 杜远生, 马千里, 等. 秭归盆地峡口剖面中侏罗统千佛崖组古气候演化的地球化学记录[J]. 地质科技情报, 2018, 37(2):50-58 DAI Xianduo, DU Yuansheng, MA Qianli, et al. Palaeoclimate change of the Middle Jurassic Qianfoya formation revealed by geochemical records of Xiakou section in Zigui Basin [J]. Geological Science and Technology Information, 2018, 37(2): 50-58.
[30] 余素华, 文启忠, 张士三, 等. 中国西北地区晚第四纪黄土中镁铝地球化学与古气候意义[J]. 沉积学报, 1994, 3(12):112-116 YU Shuhua, WEN Qizhong, ZHANG Shisan, et al. The geochemistry and paleoclimate significance of magnesium and aluminium in loess of late Quaternary in Northwestern China [J]. Acta Sedimentologica Sinica, 1994, 3(12): 112-116.
[31] 吴国瑄. 珠江口盆地东部古近系沉积环境及烃源研究[R]. 同济大学内部科技报告, 2007. WU Guoxuan. Study on Paleogene Sedimentary Environment and Hydrocarbon Sources in the Eastern Pearl River Basin[R]. Tongji University, Internal technology report, 2007.
[32] Talbot M R. A review of the palaeohydrological interpretation of carbon and oxygen isotopic ratios in primary lacustrine carbonates [J]. Chemical Geology:Isotope Geoscience Section, 1990, 80(4): 261-279. doi: 10.1016/0168-9622(90)90009-2
[33] 刘庆. 渤海湾盆地东营凹陷烃源岩碳氧同位素组成及地质意义[J]. 石油实验地质, 2017, 39(2):247-252 doi: 10.11781/sysydz201702247 LIU Qing. Composition and geologic significance of carbon and oxygen isotopes in hydrocarbon source rocks, Dongying Sag, Bohai Bay Basin [J]. Petroleum Geology and Experiment, 2017, 39(2): 247-252. doi: 10.11781/sysydz201702247
[34] 刘传联. 东营凹陷沙河街组湖相碳酸盐岩碳氧同位素组分及其古湖泊学意义[J]. 沉积学报, 1998, 16(3):109-114 LIU Chuanlian. Carbon and oxygen isotopic compositions of lacustrine carbonates of the Shahejie formation in the Dongying depression and their paleolimnological significance [J]. Acta Sedimentologica Sinica, 1998, 16(3): 109-114.
[35] 孙莎莎, 姚艳斌, 吝文. 鄂尔多斯盆地南缘铜川地区油页岩元素地球化学特征及古湖泊水体环境[J]. 矿物岩石地球化学通报, 2015, 34(3):642-645 SUN Shasha, YAO Yanbin, LIN Wen. Elemental geochemical characteristics of the oil shale and the paleo-lake environment of the Tongchuan Area, Southern Ordos Basin [J]. Bulletin of Mineralogy, Petrology and Geochemistry, 2015, 34(3): 642-645.
[36] 金强, 朱光有, 王娟. 咸化湖盆优质烃源岩的形成与分布[J]. 中国石油大学学报:自然科学版, 2008, 32(4):19-23 JIN Qiang, ZHU Guangyou, WANG Juan. Deposition and distribution of high-potential source rocks in saline lacustrine environments [J]. Journal of China University of Petroleum, 2008, 32(4): 19-23.
[37] Katz B J. Factors controlling the development of lacustrine petroleum source rocks-An update [J]. AAPG Studies in Geology, 1995, 40: 61-79.
[38] 秦建中. 中国烃源岩[M]. 北京: 科学出版社, 2005. QIN Jianzhong. Source Rocks in China[M]. Beijing: Science Press, 2005.
[39] 王建, 王权, 钟雪梅, 等. 二连盆地优质烃源岩发育特征及成藏贡献[J]. 石油实验地质, 2015, 37(5):641-647 doi: 10.11781/sysydz201505641 WANG Jian, WANG Quan, ZHONG Xuemei, et al. Characteristics of high-quality hydrocarbon source rocks and their contributions to reservoirs in the Erlian Basin [J]. Petroleum Geology and Experiment, 2015, 37(5): 641-647. doi: 10.11781/sysydz201505641
计量
- 文章访问数: 2238
- HTML全文浏览量: 589
- PDF下载量: 66
