安徽宣城第四纪网纹红土的磁组构特征及其意义

朱宗敏, 杨文强, 林文姣, 刘育燕

朱宗敏, 杨文强, 林文姣, 刘育燕. 安徽宣城第四纪网纹红土的磁组构特征及其意义[J]. 海洋地质与第四纪地质, 2006, 26(4): 105-110.
引用本文: 朱宗敏, 杨文强, 林文姣, 刘育燕. 安徽宣城第四纪网纹红土的磁组构特征及其意义[J]. 海洋地质与第四纪地质, 2006, 26(4): 105-110.
ZHU Zong-min, YANG Wen-qiang, LIN Wen-jiao, LIU Yu-yan. MAGNETIC FABRIC CHARACTERISTICS AND SIGNIFICANCE OF RETICULATE RED EARTH IN XUANCHENG PROFILE, ANHUI PROVINCE[J]. Marine Geology & Quaternary Geology, 2006, 26(4): 105-110.
Citation: ZHU Zong-min, YANG Wen-qiang, LIN Wen-jiao, LIU Yu-yan. MAGNETIC FABRIC CHARACTERISTICS AND SIGNIFICANCE OF RETICULATE RED EARTH IN XUANCHENG PROFILE, ANHUI PROVINCE[J]. Marine Geology & Quaternary Geology, 2006, 26(4): 105-110.

安徽宣城第四纪网纹红土的磁组构特征及其意义

基金项目: 

中国地质大学(武汉)校自然科学基金项目(CUGQNL0526)

国家自然科学基金项目(40502015)

详细信息
    作者简介:

    朱宗敏(1978-),女,助教,硕士,从事岩石磁学研究,E-mail:zhumin@cug.edu.cn

  • 中图分类号: P539.3

MAGNETIC FABRIC CHARACTERISTICS AND SIGNIFICANCE OF RETICULATE RED EARTH IN XUANCHENG PROFILE, ANHUI PROVINCE

  • 摘要: 安徽宣城第四纪网纹红土剖面黄棕色土层和网纹层的KPEFLq等磁组构参数特征显示,i3>60°样品的磁组构特征有效降低了后期风化作用的影响,可近似代表黄棕色土层和网纹层的原生磁组构特征。黄棕色土层的原生磁组构特征显示它形成于沉积动力强度弱而稳定、事件沉积偶尔发生的沉积环境;该层磁组构参数的平均值及组合参数特征均与风成黄土-古土壤沉积的相应特征十分吻合,表明其原始沉积环境与风尘沉积环境一致,从而为黄棕色土层的风积成因提供了沉积动力特征上的证据。网纹层形成于沉积动力强度大、稳定性相对较差的沉积环境;该层的磁组构特征与风成沉积的相应特征差别显著,但与部分水成沉积的相应特征有近似之处,结合该层物质组成上的风成特性,初步推断网纹层可能为早期的风积物经后期水流改造的产物。
    Abstract: Magnetic fabric investigation of yellow-brown layer and reticulate layer of reticulate red earth section in Xuancheng, southern Anhui Province has been carried out and data of magnetic fabric parameters K,P,E,F,L and q of yellow-brown layer and reticulate layer are gained. The characteristics of the above parameters indicate that the original magnetic fabrics of the yellow-brown layer and reticulate layer have been recomposed by weathering in different ways, while the former seemed to be more seriously reworked. The original magnetic fabrics implicate that the yellow-brown layer was deposited in weak but stable force condition with event deposition happening occasionally. However, the reticulated layer was formed in strong and unstable force conditions.
    The characteristics of mean values (K,P,E,F,L and q)and collaborated parameters(F-L,P-q) of yellow-brown layer are coincident with those of eolian loess-paleosol, indicating its eolian deposit formation. However, the characteristics of the corresponding values in the reticulated layer differ significantly from those of eolian deposit, but somewhat close to those of hydraulic deposit. It probably implicate that the reticulated layer has been recomposed by fluviation.
  • 下扬子板块是我国南方海相碳酸盐岩分布的重要地区,从晚震旦世开始到三叠纪结束长达7亿年的海相沉积演化过程中形成了广泛而巨厚的海相沉积[1]。南黄海是下扬子板块向东延伸的主体,北部以苏鲁造山带为界与中朝块体相邻,南部以江绍断裂带为界与华南块体相望,西接下扬子苏北盆地[2]。南黄海盆地经历过加里东、印支、燕山和喜山等一系列的构造运动,现今盆地的构造格局北以千里岩隆起为界,南以勿南沙隆起为界,盆地从北到南可以划分为烟台坳陷、崂山隆起、青岛坳陷3个二级构造单元[3-5](图 1)。

    图  1  南黄海盆地及相邻的大地构造格局
    Figure  1.  Tectonic map of South Yellow Sea basin and its adjacent areas

    对比南黄海盆地和下扬子陆域的钻井及地震资料,中-古生界(以T2界面为界线)从老到新主要发育以下几套主要的碳酸盐岩:上震旦统灯影组(Z2d)潮坪-局限台地相的白云岩,中-上寒武统(C - 2-3)盆地相-台缘斜坡相泥质岩,台地相白云岩及白云质灰岩,奥陶系(O)浅海陆棚相碳酸盐岩,石炭系-下二叠统栖霞组(C-P1q)潮坪-台地相碳酸盐岩和三叠统青龙组(T1q)的碳酸盐岩[6](图 2)。截至目前,下扬子陆域海相地层发现3个油气田:句容油田、朱家墩油气田和黄桥气田,而南黄海油气调查工作已有50余年的历史[7],至今尚未获得油气的重大突破[4, 8, 9]。实际勘探过程中发现碳酸盐岩储层非均质性强、储集性能差以及后期改造较为严重的不利因素,造成地震响应差异较大,因此,目前对碳酸盐岩储层预测是南黄海油气勘探的重点和难点。

    图  2  南黄海地层综合柱状图
    (据文献6)
    Figure  2.  Integrated stratigraphic column of the South Yellow Sea

    南黄海中-古生界海相碳酸盐岩储层主要受沉积、成岩和构造作用影响,前人对南黄海碳酸盐岩储层预测主要利用地震属性变化[3, 10]和速度分析[2, 11-13]方法,在一定程度上能够将碳酸盐岩储层进行有效预测,但是利用叠后地震技术是基于零偏移距道集叠加的地震资料分析,单纯应用叠后地震进行属性分析和叠后波阻抗反演会增加储层预测的多解性以及难以对有效储层进行识别[14],而叠前地震资料反演是基于能够反映振幅随偏移距变化的多个角度叠加的叠前地震数据体,需充分利用测井曲线进行叠前同时反演,得到纵波阻抗、横波阻抗和密度体[15, 16],并通过相关公式得到一系列岩性指示因子和物性指示因子等弹性参数体,才能有效地预测有利储层。本文通过对南黄海碳酸盐岩储层进行岩石物理分析基础上,提取有利储层的敏感弹性参数,采用叠前同时反演方法对研究区碳酸盐岩的有利储层进行有效的预测。

    叠前同时反演是将叠前CRP地震道集中按照远、中、近偏移距进行部分叠加,同时反演出纵波阻抗(ZP)、横波阻抗(ZS)和密度(ρ),根据Zoepprize方程[17]演化而来:

    $$\left[ {\begin{array}{*{20}{c}} {\sin {\alpha _1}}&{\cos {\beta _1}}&{ - \sin {\alpha _1}}&{\cos {\beta _2}}\\ {\cos {\alpha _1}}&{ - \sin {\beta _1}}&{\cos {\alpha _2}}&{\sin {\beta _2}}\\ {\cos 2{\beta _1}}&{\frac{{{\upsilon _{s1}}}}{{{\upsilon _{\rho 1}}}}\sin 2{\beta _1}}&{\frac{{{\rho _2}{\upsilon _{s2}}}}{{{\rho _1}{\upsilon _{\rho 1}}}}\cos 2{\beta _2}}&{\frac{{{\rho _2}{\upsilon _{s2}}}}{{{\rho _1}{\upsilon _{\rho 1}}}}\sin 2{\beta _2}}\\ {\frac{{\upsilon _{s1}^2}}{{{\upsilon _{\rho 1}}}}\sin 2{\alpha _1}}&{{\upsilon _{s1}}\cos 2{\beta _1}}&{\frac{{{\rho _2}2\upsilon {{_s^2}_2}}}{{{\rho _1}{\upsilon _{p2}}}}\sin 2{\alpha _2}}&{\frac{{\rho {\upsilon _{s2}}}}{{{\rho _1}}}\cos 2{\beta _2}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{R_p}}\\ {{R_s}}\\ {{T_p}}\\ {{T_s}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} { - \sin \alpha }\\ {\cos {\alpha _1}}\\ { - \cos 2{\beta _1}}\\ {\frac{{\upsilon _{s1}^2}}{{{\upsilon _{\rho 1}}}}\sin 2{\alpha _1}} \end{array}} \right] $$ (1)

    式中$ {\upsilon _{\rho 1}},{\upsilon _{s1}},{\rho _1}$为上部地层的纵波速度、横波速度和密度;${\upsilon _{\rho 2}},{\upsilon _{s1}},{\rho _2}$为下部地层的纵波速度、横波速度和密度;RP为纵波反射系数,RS为横波反射系数,TP为纵波透射系数,TS为横波透射系数,α1β1α2β2分别为纵波入射角、横波反射角、纵波透射角和横波透射角。

    1980年Aki和Richard对Zoeppritz进行简化,得到Aki-Richard近似表达式[18],得到了当弹性参数的相对变化$ \frac{{\Delta {\upsilon _\rho }}}{{{\upsilon _\rho }}}$、$\frac{{\Delta {\upsilon _s}}}{{{\upsilon _s}}} $、$ \frac{{\Delta \rho }}{\rho }$较小时,利用纵波速度、横波速度和密度的相对变化表示的反射系数公式[19],即:

    $$ {R_p}\left( \theta \right) = \frac{1}{2}\left( {1 - 4\frac{{\upsilon _s^2}}{{\upsilon _p^2}}{{\sin }^2}\theta } \right)\frac{{\Delta \rho }}{\rho } = \frac{{{{\sec }^2}\theta }}{2}\frac{{\Delta {\upsilon _\rho }}}{{{\upsilon _\rho }}} - 4\frac{{\upsilon _s^2}}{{\upsilon _\rho ^2}}\theta {\sin ^2}\theta \frac{{\Delta {\upsilon _s}}}{{{\upsilon _s}}} $$ (2)

    式中vpvsρ分别为界面两侧的纵波速度、横波速度和密度平均值,${\upsilon _R} = \frac{{{\upsilon _{P1}} + {\upsilon _{P1}}}}{2},{\upsilon _s} = \frac{{{\upsilon _{s1}} + {\upsilon _{s1}}}}{2},\rho = \frac{{{\rho _1} + {\rho _2}}}{2} $; ΔvP,Δvs,Δρ分别为界面两侧纵波速度、横波速度和密度的差值,$ \Delta {\upsilon _p} = {\upsilon _{p2}} - {\upsilon _{p1}},\Delta {\upsilon _{s2}} - {\upsilon _{s1}},\Delta \rho = {\rho _2} - {\rho _1}。\theta = \frac{{{\theta _i} + {\theta _t}}}{2}$为入射纵波的入射角和透射纵波的透射角平均值,θi为入射纵波的入射角,θt为透射纵波的透射角,Aki-Richards方程可以写成Fatti方程[20]的形式:

    $$ {R_{pp}}\left( \theta \right) = {c_1}{R_P} + {c_2}{R_p} + {c_3}{R_D} $$ (3)

    其中${c_1} = 1 + {\tan ^2}\theta$,Rp为法向入射时的纵波反射系数,${R_p} = \frac{1}{2}\left[ {\frac{{\Delta {\upsilon _p}}}{{{\upsilon _p}}} + \frac{{\Delta \rho }}{\rho }} \right],{c_2} = - 8{\gamma ^2}\sin \theta$,RS为横波反射系数,$ {R_S} = \frac{1}{2}\left[ {\frac{{\Delta {\upsilon _s}}}{{{\upsilon _s}}} + \frac{{\Delta \rho }}{\rho }} \right],{c_3} = - \frac{1}{2}{\tan ^2}\theta {\rm{ + }}2{\gamma ^2}\sin \theta $,RD为密度反射系数,${R_D} = \frac{{\Delta \rho }}{\rho } $,其中$\gamma = \frac{{{\upsilon _s}}}{{{\upsilon _p}}}$。

    由于方程(3)的各个系数在数量级上是不同的,使得在小角度的情况下求解得到RSRD不稳定,加上ZP、ZSρ两两之间存在着相关性,因此,通过以下方式消除上述问题,假设背景为含水岩层的情况下,存在趋势关系[14, 21-23],由Castagna方程[24]常数,得到:

    $$ {\rm{In}}\left( {Zs} \right) = {\rm{In}}\left( {Zp} \right) = {\rm{In}}\left( \gamma \right) $$ (4)

    由Gardner方程[25],$ \rho = a\upsilon _p^b$,得到:

    $$ {\rm{In}}\left( \rho \right) = \frac{b}{{1 + b}}{\rm{In}}\left( {{Z_p}} \right) + \frac{{{\rm{In}}\left( a \right)}}{{1 + b}} $$ (5)

    由上述关系,可以得到更一般的背景趋势关系:

    $$ \left\{ {\begin{array}{*{20}{c}} {{\rm{In}}\left( {{Z_s}} \right) = k\;{\rm{In}}\left( {{Z_p}} \right) + {k_c} + \Delta {L_S}}\\ {{\rm{In}}\left( \rho \right) = m\;{\rm{In}}\left( {{Z_p}} \right) + {m_c} + \Delta {L_D}} \end{array}} \right. $$ (6)

    km为含水背景趋势的斜率,kcmc为含水背景趋势的截距,ΔLS和ΔLD为烃类响应,当为含水背景时,其值为0。这样,Fatti方程可以写成如下形式:

    $$ \begin{array}{*{20}{l}} {\left\{ \begin{array}{l} T\left( \theta \right) = {{\bar c}_1}W\left( \theta \right)D{L_P} + {{\bar c}_2}W\left( \theta \right)D\Delta {L_D}\\ {{\bar c}_1} = \frac{1}{2}{c_1} + \frac{1}{2}k{c_2} + m{c_3}\\ {{\bar c}_2} = \frac{1}{2}{c_2} \end{array} \right.} \end{array} $$ (7)

    其中W(θ)为入射角θ处的子波;D为导数因子,LP=ln(Zp),T(θ)为地震记录。

    改进后的方程(7)相对方程(3)具有以下优点:(1)各变量之间是相互不关联的;(2)建立了含水背景下各个变量之间的岩石物理特征关系。

    本文叠前同时反演的流程如图 3,具体步骤是在建立井点储层地质模型和岩石物理模型的基础上,利用Xu-white模型通过纵波速度求取横波速度,进行岩石物理弹性参数计算,分析研究区岩性和物性敏感弹性参数,利用叠前地震CRP道集进行精细的AVO处理和分析,寻找不同孔隙储层的响应特征[26, 27],对CRP地震道集进行分角度叠加得到远、中、近的叠加数据体,合成地震记录提取不同角度道集的子波,构建纵波阻抗(ZP)、横波阻抗(ZS)和密度(ρ)低频模型,对远、中、近的叠加数据体以及构建的低频模型代入公式(7)进行叠前同时反演,反演得到纵波阻抗、横波阻抗和密度体,并利用岩石物理关系式得到物性因子和流体因子,从而对碳酸盐岩有利储层进行解释[27]

    图  3  叠前同时反演流程图
    Figure  3.  Processing chart of pre-stack simultaneous inversion

    南黄海碳酸盐岩按沉积环境划分为颗粒灰岩、泥晶灰岩和白云岩类。颗粒灰岩是南黄海地区发育规模较大的浅水滩相沉积,在震旦系、寒武系、奥陶系和下二叠统都有发育,根据沉积相分析石炭系和下三叠统也发育小型台地点滩和潮间滩,按照不同填隙物将颗粒灰岩细分为亮晶颗粒灰岩和泥晶颗粒灰岩,其中亮晶颗粒灰岩发育种类较多,有亮晶内碎屑灰岩、细晶-微晶生物碎屑灰岩、亮晶球粒灰岩,生物碎屑灰岩。微晶灰岩的晶粒大小约为0.004~0.10mm,为细晶-微晶状,见少量晶粒>0.10mm或<0.004mm的方解石不均匀分布,晶粒间结合紧密,见重结晶现象,部分方解石呈似砂屑状,微裂缝发育,缝宽约0.01mm,泥晶方解石充填(图 4a)。灰质白云岩具细晶-粉晶结构,块状构造,矿物成分主要为白云石,粒状,粒径0.06~0.2mm,见雾心亮边结构,雾心可能是白云石交代方解石的残余物,裂缝发育被方解石、铁泥质充填(图 4b)。总的来说,下扬子地区灰岩孔隙不太发育,主要以构造缝和晶间孔为主;白云岩在震旦系、寒武系、奥陶系、上石炭统船山组和下三叠统青龙组都有发现,且溶孔发育,孔隙条件较好[4]。不仅在灰岩,还有在白云岩层段,都能看到溶蚀孔,特别在构造活跃层段,能见到较多的角砾,沿着角砾发生溶蚀作用,角砾间溶蚀扩大孔较为发育,其次为少量的晶间溶蚀扩大孔(图 4c),另外还发育构造缝,缝中半充填,且有微弱的溶蚀扩大,晶间孔以白云石或方解石晶间孔为主(图 4d)。

    图  4  下扬子中-古生界海相碳酸盐岩铸体薄片显微照片
    (a)石炭系船山组, 含生屑细晶G微晶灰岩, 细晶G微晶状, 生物碎屑主要为介屑、棘屑等, 矿物成分为方解石, 不均匀分布, 晶粒大小约为0.004~0.16mm, CZ12G1G1井, 2074?m; (b)寒武系冷泉王组, 灰质白云岩, 孔隙不发育, 白云石矿物粒径0.06~0.2mm, 在脉体中偶见方解石溶蚀孔, 安徽巢湖汤山; (c)震旦系灯影组白云岩中的溶蚀孔, 粒屑结构, 块状构造, 粒径0.3~3.5mm, 边缘不清晰, 围绕团块具栉壳结构, 原生孔隙中中晶、细晶方解石充填, 巢湖青苔山; (d)下奥陶统伦山组, 灰岩中的构造缝, 缝宽约0.01~0.02mm, 部分被方解石充填
    Figure  4.  Micrographs of Mesozoic-Paleozoic carbonates in the Lower Yangtze areas

    据南黄海井点物性条件统计,下三叠统上青龙组灰岩孔隙度为2.7%~5.6%,平均为4.01%,渗透率<0.01mD;下青龙组灰岩孔隙度为0.8%~1.6%,平均为1.12%,渗透率平均0.04mD;下二叠统栖霞组灰岩孔隙度在0.46%~2.05%之间,平均为1.24%,渗透率<0.01mD;上石炭统船山组藻灰岩孔隙度为0.86%~3.55%,平均为1.81%,渗透率<0.01mD,总体上碳酸盐岩储层为低孔低渗储层。但局部存在高孔高渗层段,如石炭系顶部和下二叠统顶部的风化壳和淋滤带为白云岩发育段,CZ12-1-1井船山组顶部风化剥蚀面孔隙度高达14%,CZ35-2-1井栖霞组发育孔隙为6.6%~10.4%,WX13-3-1井栖霞组顶部孔隙度为3.33%~20.19%,WX5-ST1井下青龙组灰岩中发育25m厚高孔渗白云岩,孔隙度为6%~8%,CZ12-1-1井和州组白云质灰岩孔隙度在4.1%~7.8%。另外在孔、洞、缝发育的碳酸盐岩地层,其储层的孔隙条件较好,易构成较有利的储层。

    一般情况下,白云岩比灰岩的声波阻抗大,当白云岩中有孔隙发育时,其阻抗会降低,从而与灰岩的阻抗相近;另外碳酸盐风化壳一般具有较高的速度和密度,由于孔、洞、缝的存在造成速度和密度的下降。基于以上碳酸盐岩的非均质性,使得储层与非储层的纵波阻抗频率直方图相互叠置,从而导致利用常规纵波阻抗反演方法难以区分储层和非储层。

    通常情况下,不同的岩石具有不同的物理特性,岩石物理特性与地震特征有直接的关系[23, 28],岩石物理分析是叠前岩性与物性反演的基础,也是连接地震-测井-储层物性参数之间的桥梁和纽带,为了更好地开展碳酸盐岩储层预测工作,首先要对不同储层类型进行岩石物理分析工作,优选有利储层敏感的弹性参数,为有利储层预测工作指明方向[29]

    横波资料是岩石物性研究的基础,叠前地震资料反演重要参数,而南黄海测井曲线缺少横波,因此,需要对横波进行预测,一般是通过岩石物理建模、流体替代以及经验公式模型等多种技术手段预测横波,获得各种弹性模量和弹性参数,为正演模拟、AVO属性分析,叠前同时反演提供基础参数和信息。以Xu-White岩石物理模型为基础,根据骨架和泥质不同的纵横直径比以及其他有关参数来估算岩石的体积模量和剪切模量,并按照孔隙中流体的性质来估计流体对岩石模量的影响,再此由密度测井的密度测井曲线和Xu-White模型得出的体积模量和剪切模量,便可导出岩石的纵、横波速度[30]

    图 5中第4道蓝色线为实测的纵波曲线,红色线为预测的纵波曲线,可以看出实测和预测的纵波曲线吻合较好,两者的相对误差整体较小,为0~0.1%,在浅部地层两者的相对误差最小,仅为0~0.05%,而随着深度的加大,误差也增大,达到0~0.1%,但误差都在可容范围之内[30]图 5第6道为最终求得横波速度,其范围值为2000~4000m/s,符合陆上苏北三叠系地层碳酸盐岩地层的纵、横波速度分布,说明通过Xu-White岩石物理模型进行横波预测具有一定的岩石物理规律,预测横波速度相对可靠,能够满足叠前反演需求,为后续的AVO分析和叠前同时反演奠定了基础。

    图  5  Xu-White模型预测CZ35-2-1井的纵横波速度曲线
    Figure  5.  The predicted P&s wave velocity logging of CZ35-2-1well by Xu-White model

    通过对横波速度进行预测,根据岩石物理方程,得到泊松比、拉梅系数、体积模量、剪切模量、杨氏模量等系数,另外通过对纵横波阻抗合适的坐标旋转以及转换,构建物性指示因子,关系式如下:

    $$ Y = {Z_p} \times \cos \left( \theta \right) + {Z_S} \times \sin \left( \theta \right) $$ (8)

    研究发现,新的指示因子随着θ角度的改变而变化,在不同的角度上可以与传统的储层或弹性参数具有较高的相关性,由以上方法我们可以得到物性指示因子(RPI):与孔隙度参数相关最大的角度所计算出的阻抗组合。此种方法具有至少具有3个方面的先进性:首先是增大了划分岩性和物性性质的区分度和识别率;其次是把多参数交会识别有利区降维至单参数识别;最后是一定程度上克服了背景趋势给识别有利储层带来的干扰。

    根据CZ35-2-1井碳酸盐岩和碎屑岩的测井曲线(图 6),可以看出碳酸盐岩相对碎屑岩具有高密度、高泊松比、高拉梅系数、高剪切模量和杨氏模量、高纵横波速度、高纵横波速度比和低自然伽马,在碳酸盐岩有效储层段为低密度、物性指示因子高。

    图  6  CZ35-2-1井碳酸盐岩和碎屑岩的测井曲线
    Figure  6.  Well logging of carbonate and clastic of well CZ35-2-1

    对碳酸盐岩物性的敏感弹性进行分析,通过密度计算孔隙度,发现密度越低,孔隙度越大的特点,但由于本区目的层埋深大,还有次生孔隙缘故,密度与孔隙度之间的线性关系不明显。将物性指示因子与剪切模量进行交汇图(图 7)分析,图中以孔隙度为色标,可以看出物性指示因子越高,孔隙度越大,因此,物性指示因子为研究区块物性敏感弹性参数,利用物性指示因子能够圈定碳酸盐岩有利储层。

    图  7  CZ35-2-1井碳酸盐岩物性敏感弹性参数交汇图
    Figure  7.  Crossplot of carbonate physical property sensitive elastic parameters of well CZ35-2-1

    在叠前时间偏移处理后的CRP道集上,可以看到很多地震同相轴没有拉平的现象,如图 7a所示,造成这种现象的原因很多,地层横向速度变化是其中的一个因素。这些现象会给叠前反演后属性分析等带来较大影响,因此,需要进行剩余时差校正,开展道集拉平处理。具体方法是:(1)选择随机噪音小的、靠近小偏移距的一组道进行部分叠加形成参考道;(2)选择一个对于所有偏移距至少包含一个波组的时窗;将该时窗内的所有道与参考道进行互相关,求得每道最大相关系数所对应的时移量;(3)根据相关系数的门槛值和时移量门槛值剔除掉异常时移值后,得到各道的校正时移量;(4)将校正时移量应用于该时窗数据;(5)时窗沿时间方向滑动至道尾完成一个道集的处理工作,逐道集处理完成整个工区道集拉平处理工作[31]。通过一系列的优化处理,如图 7b充分挖掘叠前地震资料信息的潜力,恢复叠前地震信息,尽可能提高信噪比,使道集质量获得明显的提高,能量更加集中,振幅随偏移距变化关系更加明显,这样的道集更有利于AVO分析和叠前同时反演[29]

    图  8  原始CRP地震道集(a)和优化后CRP地震道集(b)
    Figure  8.  Orignial CRP seismic gather(a)and optimized CRP seismic gather(b)

    子波是参与反演运算一个重要参数,进行叠前反演时,通常隐含一个假设,即子波不随偏移距变化(包括振幅、相位和带宽),然而实际上子波是随着偏移距变化的,因此,在进行叠前反演是对每个部分角度叠加地震道集都用它自身合成地震记录的子波来转化。进行子波提取计算的时窗主要在目的层内,子波主频与目的层地震主频一致,子波长度一半为地震采样间隔的整数倍,根据子波的形态、振幅谱及相位谱对提取的子波进行质量分析,提取的子波具有波形稳定且频域内为单峰的特征[32]。本次采用了4个部分叠加地震道集方案(图 9a),并且对井进行分角度提取子波,各个子波形态如图 9b

    图  9  分角度叠加合成地震记录(a)及提取的子波(b)
    Figure  9.  Angle stacking synthetic seismogram(a) and the extraction of wavelet(b)

    由于大地低通滤波作用,实际采集到的地震数据频带较窄,不仅难以获得较高的频率成分,同时也缺少甚低频数据,低频模型的建立也是叠前同时反演的一个重要环节。如何建立低频模型,主要由两个关键步骤组成:构建构造格架和采用插值方式对格架进行属性填充[33]。其中建立构造格架是分析地层的沉积模式和层位接触关系,按照地层沉积的顺序从下向上逐层依次搭建,由于南黄海存在多期的逆冲断层构造,上下盘地层叠置严重,断层倾角变化较大而严重影响了阻抗反演的精度,本文采用了层位分开建模的方法,并通过将断层转换为层位与地层相融合来建立初始的地层格架模型[34]。同时由于目前勘探的钻井最深钻到志留系地层,对于往下的地层属性无法确定的情况下,对CRP地震道集进行层析速度分析的基础上提取虚拟井,对以上横波预测结果与纵波进行关系式拟合,求取虚拟井的横波速度,并采用Gardenr公式将纵波速度转化为密度,将虚拟井进行合成地震记录后,对地震进行带限阻抗反演,经过多次迭代,最终得到的纵波阻抗低频模型(图 10a)、横波阻抗低频模型(图 10b)和密度低频模型(图 10c),可以看出,建立的低频模型符合本区碳酸盐岩和碎屑岩的纵波阻抗、横波阻抗以及密度的分布。

    图  10  层析速度分析和虚拟井带限阻抗反演结合建立的低频模型
    (a)低频纵波阻抗模型;(b)低频横波阻抗模型;(c)低频密度模型
    Figure  10.  Low frequency model combined with tomographic velocity analysis and the virtual

    利用分角度叠加地震道集资料,以及纵波阻抗、横波阻抗和密度低频模型的约束下进行叠前同时反演,得到纵波阻抗、横波阻抗和密度体,根据公式(8)转换得到物性指示因子(RPI)剖面(图 11),红色代表物性好,蓝色代表物性差,黄色的物性处于中等。由图中可以看出总体上碳酸盐岩的物性较好,尤其是位于古高地的局部发育风化壳、淋滤带、白云岩,其孔、洞、缝比较发育,能够提高碳酸盐岩储层的孔隙条件,易形成有利的储层。上震旦统灯影组(Z2d)白云岩处于蒸发潮坪相,以溶蚀孔、粒间孔和粒内孔的储集空间发育,物性较好;中-上寒武统—奥陶系(C - 2-3—O)灰岩除了发育风化壳,还有溶蚀孔和洞缝型储层,以及后期构造裂缝为主要储集空间,缝洞呈层状分布;石炭系-下二叠统(C-P1q)碳酸盐岩储层风化壳和滩相都较为发育,发育有瘤状灰岩和生物碎屑灰岩,石炭系和栖霞组顶面由于海西运动遭受风化剥蚀,存在岩溶风化壳,具有一定的孔隙发育。

    图  11  叠前反演结果的物性剖面
    Figure  11.  The physical property section from pre-stack inversion

    总的来说,叠前同时反演结果反映的物性剖面具有横向变化分布特征,由于受古构造等地质因素影响,位于古高地为局限台地相沉积,处于高潮面和最大风暴面之间,大部分时间暴露于海平面之上,受气候影响比较大,海水蒸发引起盐度高,并且水循环受限制,灰岩地层容易发生交代作用形成间歇性白云岩,后期由于淡水淋滤作用,易形成岩溶型储层,提高孔隙空间,因此古高地的孔隙度相对较大,位于凹陷带主要沉积了泥晶灰岩或者泥质灰岩,水体能量低,可能为安静的海相或泻湖相沉积,孔隙度相对较小。碳酸盐地层由于后期的构造运动形成或裂缝也会改造储层的孔隙空间,使得孔隙度也有所提高。

    (1) 根据岩石物理分析结果,发现碳酸盐岩相对碎屑岩具有高阻抗和高泊松比,物性指示因子对碳酸盐岩孔隙度比较敏感,物性指示因子越高,孔隙度越大,因此,物性指示因子可作为碳酸盐岩物性定量判别标准;

    (2) 通过对南黄海碳酸盐岩叠前同时反演,其结果对碳酸盐岩物性具有很好的指示作用,物性较好的局部地区说明发育有风化壳、林滤带和白云岩地层,由于孔、洞、缝比较发育,碳酸盐岩的孔隙条件发育较好,易形成有利储层;

    (3) 碳酸盐岩的分布受沉积、成岩和构造控制作用,位于古高地的局限台地相碳酸盐岩由于交代作用形成白云岩,由于淡水淋滤作用,易形成岩溶型储层,加上后期的构造运动,裂缝的发育一定程度上提高孔隙空间。

    (4) 叠前同时反演方法依赖于精确的速度,因研究区地震资料品质相对较差,深部钻井较少,特别是缺少震旦-奥陶系的钻井资料,深部储层反演存在一定误差。若相关资料不断完善、品质不断提升,可以进一步提升反演精度。

  • [1] 刘东生. 黄土与环境[M]. 北京:科学出版社,1985:191-400.[LIU Tung-sheng. Loess and Environment[M]. Beijing:Science Press, 1985:191

    -400.]

    [2]

    Guo Z T,Ruddiman W F,Hao Q Z,et al. Onset of Asian desertification by 22Myr ago inferred from loess deposit s in China[J]. Nature,2002,416:159-163.

    [3] 李四光. 冰期之庐山[M]//李四光主编.前国立研究院地质研究所专刊(乙种第2号).南京地质研究所,1947:28-33.[LI Si-guang. Quaternary glaciations in the Lushan area[M]//LI Si-guang ed. Contributions from the National Institute of Geology Academia Sinica. Nanking:Institute of Geology,Academia Sinica,1947:B (2

    ):28-33.]

    [4] 李驭亚.华南第四纪网纹红土虫状白斑的成因探讨[J]. 地质论评,1965,23(2):144-145.

    [LI Yu-ya. Study on the formation of worm-like white mottles in the Quaternary reticulate red clay in South China[J]. Geological Reviews, 1965,23(2):144-145.]

    [5] 孙殿卿. 中国第四纪亚冰期划分之商讨[J]. 第四纪研究,1989,9(1):3-23.

    [SUN Dian-qing. The discussion of sub-ice age division in the fourth age in China[J]. Quaternary Science, 1989,9(1):3-23.]

    [6] 于振江,黄多成. 安徽省沿江地区网纹红土和下属土的形成环境及其年龄[J]. 安徽地质,1996,6(3):48-55.

    [YU Zhen-jiang, HUANG Duo-cheng. Formation environment of net-veined laterite and Xiashu loess and their ages in the area along the Yangtze River, Anhui Province[J]. Geology of Anhui, 1996, 6(3):48-55.]

    [7] 朱诚. 对庐山东麓第四系沉积物物源及新构造运动的新认识[J].地理学报, 1995,50(1):41-50.

    [ZHU Cheng. The new recognitions to genesis Quaternary sediments and neotectonics in the east piedmont of the Lushan mountain[J]. Acta Geographica Sinica, 1995, 50(1):41-50.]

    [8] 朱显谟.中国南方的红土与红色风化壳[J]. 第四纪研究,1993,13(1):75-84.

    [ZHU Xian-mo. Red clay and red residuum in South China[J]. Quaternary Sciences, 1993,13(1):75-84.]

    [9] 杨达源,韩辉友,周旅复,等. 安徽宣城地区中晚更新世风成堆积与环境变迁[J]. 海洋地质与第四纪地质,1991,11(2):97-104.

    [YANG Da-yuan, HAN Hui-you, ZHOU Lü-fu, et al. Eolian deposit and environmental change of middle-late Pleistocene in Xuancheng, Anhui Province[J]. Marine Geologgy and Quaternary Geology, 1991, 11(2):97-104.]

    [10] 杨达源. 中国东部的第四纪风尘堆积与季风变迁[J].第四纪研究,1991,11(4):354-360.

    [YANG Da-yuan. The Quaternary dust-fall accumulation and monsoon variability in Eastern China[J]. Quaternary Sciences, 1991,11(4):354-360.]

    [11]

    Xiong S F, Sun D H, Ding Z L. Aeolian origin of the red earth in southeast China[J]. Journal of Quaternary Science,2002,17(2):181-191.

    [12]

    Hu X F, Yuan G D, Gong Z T. Origin of Quaternary red clay of southen Anhui Province[J]. Pedoshere,1998,8(3):267-272.

    [13] 胡雪峰,龚子同.江西九江泰和第四纪红土成因的比较研究[J].土壤学报,2001,38(1):1-9.

    [HU Xue-feng, GONG Zi-tong. Comparative study on the origin of Quaternary red earth in Jiujiang and Taihe,Jiangxi Province[J]. Acta Pedologica Sinica, 2001,38(1):1-9.]

    [14] 胡雪峰,朱煜,沈铭能. 南方网纹红土多元成因的粒度证据[J]. 科学通报,2005,50(9):918-925.

    [HU Xue-feng, ZHU Yu,SHEN Ming-neng. Grain-size evidence for the multi-genesis of reticulated clay in Southern China[J]. Chinese Science Bulletin,2005,50(9):918-925.]

    [15] 夏应菲,杨浩. 安徽宣城第四纪红土剖面石英颗粒扫描电镜观察[J]. 南京师大学报自然科学版,1998,21(1):120-124.

    [XIA Ying-fei, YANG Hao. SEM Scanning of Quartz of the Quaternary Red Earth in Xuancheng,Anhui[J].Journal of Nanjing University(Natural Science), 1998,21(1):120-124.]

    [16]

    LIU X, XU T,LIU T. The Chinese loess in Xifeng,Ⅱ, A study of anisotropy of magnetic susceptibility of loess from Xifeng[J]. Geophys.Journal, 1988, 92:349-353.

    [17] 张家强,丛友滋,李从先. 水成沉积与风成沉积及古土壤的磁组构特征[J]. 海洋地质与第四纪地质,1999,19(2):85-94.

    [ZHANG Jia-qiang, CONG You-zi,LI Cong-xian. Magnetic fabric characteristics of hydraulic deposit, eolian deposit and Paleosoil[J]. Marine Geology and Quaternary Geology, 1999, 19(2):85-94.]

    [18] 丛友滋,韩国忠,张家强,等. 标志沉积层磁组构特征及其古环境意义研究[J]. 海洋地质与第四纪地质,1997,17(3):31-36.

    [CONG You-zi, HAN Guo-zhong, ZHANG Jia-qiang, et al. Studies on the magnetic fabric characteristics of index bed and the paleoenvironmental significance[J]. Marine Geology and Quaternary Geology, 1997, 17(3):31-36.]

    [19] 陈庆强,李从先,丛友滋. 沉积物磁组构与其动力沉积特征对应关系研究[J]. 科学通报. 1998,43(10):1106-1109.

    [CHEN Qing-qiang, LI Cong-xian, CONG You-zi. Relationship between magnetic fabric parameters and dynamic characteristics of deposits in different environments[J]. Chinese Science Bulletin, 1998,43(22):1917-1921.]

    [20] 阎桂林. 岩石磁化率各向异性在地学中的运用[M]. 武汉:中国地质大学出版社,1996:8-89.[YAN Gui-lin. Anisotropy of magnetic susceptibility of rocks and its application in geology[M]. Wuhan:China University of Geosciences Press, 1996:8

    -89.]

    [21] 刘宝柱,李从先,业治铮,等. 沉积磁组构在长江三角洲晚第四纪古土壤研究中的运用[J]. 海洋地质与第四纪地质,1994,14(2):55-62.

    [LIU Bao-zhu, LI Cong-xian, YE Zhi-zheng, et al. Applications of magnetic fabrics to the studies of Late Quaternary paleosoils in the Changjiang River delta area[J]. Marine Geology and Quaternary Geology, 1994,14(2):55-62.]

    [22] 刘良梧,龚子同.宣城第四纪红色粘土剖面的发育特征[J].第四纪研究,2000,20(5):465-468.

    [LIU Liang-wu, GONG Zi-tong. Characteristics of development of Quaternary red clay in Xuancheng, Anhui Province[J]. Quaternary Sciences, 2000,20(5):465-468.]

    [23] 李徐生,杨达源,鹿化煜.皖南风尘堆积序列氧化物地球化学特征与古气候记录[J].海洋地质与第四纪地质,1999,19(4):75-82.

    [LI Xu-sheng, YANG Da-yuan, LU Hua-yu. Oxide-geochemistry features and paleoclimatic record of the aeolian-dust depositional sequence in southern Anhui[J]. Marine Geology and Quaternary Geology,1999,19(4):75-82.]

    [24] 黄孝刚,孙继敏.末次间冰期以来黄土-古土壤序列的磁组构特征及其指示的古风向[J].第四纪研究,2005,25(4):516-522.

    [HUANG Xiao-gang, SUN Ji-min.Study of the magnetic fabrics in Chinese loess-paleosols since the last interglacial:implication of the paleowind direction[J].Quaternary Sciences, 2005,25(4):516-522.]

    [25] 许同春.磁组构(下)[J]. 地震地磁观测与研究,1989,10(3):72-80.

    [XU Tong-chun. Magnetic fabrics(Ⅱ)[J].Seismological and Geomagnetic Observation and Research, 1989, 10(3):72-80.]

    [26] 吴海斌,陈发虎,王建民,等.干旱区湖泊沉积物磁组构参数量值特征与环境变化研究[J]. 海洋地质与第四纪地质,1999,19(2):95-102.

    [WU Hai-bin, CHEN Fa-hu, WANG Jian-min, et al. Study on magnetic fabric parameters of lake sediments in arid area and environmental changes[J]. Marine Geology and Quaternary Geology,1999,19(2):95-102.]

    [27] 陈庆强,李从先. 长江三角洲地区晚更新世硬粘土层成因研究[J]. 地理科学,1998,18(1):53-57.

    [CHEN Qi-qiang, LI Cong-xian. Study on origin of the late Pleistocene stiff clays at the Yangtz delta[J]. Science Geographica Sinica, 1998, 18(1):53-57.]

  • 期刊类型引用(7)

    1. 胡雪莹,龚锐,高顺莉,张建明,刘玉柱,耿建华. 南黄海中部隆起带折射波层析成像. 地球物理学进展. 2023(02): 803-813 . 百度学术
    2. 陈建文,杨长清,张莉,钟广见,王建强,吴飘,梁杰,张银国,蓝天宇,薛路. 中国海域前新生代地层分布及其油气勘查方向. 海洋地质与第四纪地质. 2022(01): 1-25 . 本站查看
    3. 张然,德勒恰提·加娜塔依. 基于卷积神经网络的油气地震储层预测研究. 计算机仿真. 2022(01): 471-475 . 百度学术
    4. 王亚,易远元,王成泉,王孟华,杨洲鹏,张红文,王盛亮,贾敬. 叠后反演技术在杨税务潜山裂缝孔隙型储层预测中的应用. 物探与化探. 2020(05): 1208-1214 . 百度学术
    5. 吴淑玉,刘俊,陈建文,梁杰,张银国,袁勇,许明. 南黄海崂山隆起石炭系—下二叠统孔隙型碳酸盐岩储层预测. 海洋地质与第四纪地质. 2020(05): 136-148 . 本站查看
    6. 陈建文,梁杰,张银国,杨长清,袁勇,许明,王建强,雷宝华,李刚,杨艳秋,杨传胜,孙晶. 中国海域油气资源潜力分析与黄东海海域油气资源调查进展. 海洋地质与第四纪地质. 2019(06): 1-29 . 本站查看
    7. Jian-wen Chen,Ming Xu,Bao-hua Lei,Jie Liang,Yin-guo Zhang,Shu-yu Wu,Jian Shi,Yong Yuan,Jian-qiang Wang,Yu-xi Zhang,Gang Li,Wen-juan Wang. Prospective prediction and exploration situation of marine Mesozoic-Paleozoic oil and gas in the South Yellow Sea. China Geology. 2019(01): 67-84 . 必应学术

    其他类型引用(2)

计量
  • 文章访问数:  2110
  • HTML全文浏览量:  153
  • PDF下载量:  11
  • 被引次数: 9
出版历程
  • 收稿日期:  2005-11-29
  • 修回日期:  2006-02-20

目录

/

返回文章
返回