研究区位于普光气田主体的中部，西翼被普光–东岳寨断层所控制并通过该断层与土主–大湾背斜构造带相接触，受断层影响，倾角较陡，深度较大；构造东翼以浅、缓的凹陷与老君庙构造相接，构造西北侧受东岳寨–普光断裂及普光7断层遮挡，西南部的东岳寨–清溪场–老君山之间存在一条北西西向延伸的长兴组–飞仙关组沉积相变带，为一典型的构造–岩性复合圈闭。

本次研究区西起回注1井，东至老君3井，北到普光6井~普光202井一线，南至东岳井1井~清溪2井，面积100 km²。研究区构造处于北低南高的斜坡带。

浅层从须家河至地表整体表现为向西、南抬升的单斜，具西高东低、南高北低特征。区内勘探面积320 km²。普光东洼总体构造变形强度较弱，沿普光302-1-普光9-普光3井以东一线，构造相对较低。在普光1-普光3-普光7井一线以西，地层很陡，而在普光3井以东构造非常平缓。因受燕山、喜山运动影响，断裂非常发育，但规模较小。平面上分为两组，即：北西向断层和北东向断层。其中，北东向断层数量少，仅在西部高部位有少量发育，北西向断层较为发育但断距通常较小，断距一般为50~100 m，断层一般消失于千佛崖组，区内无大型应力构造，仅西部高部位靠近毛坝断层发育大湾3块、回注1块等4个小型的断块构造，面积5.07~16.29 km²，闭合幅度400~700 m，高点埋深1300~1850 m。

须家河组岩心以灰色为主，体现出水下还原沉积环境，砂岩中多发育块状层理、交错层理及平行层理。须家河组砂岩整体展现出低伽马、高电阻的电测曲线特征，测井曲线多以箱形、钟形及漏斗形为主，岩性剖面显示以砂泥岩互层为主，沉积物在纵向上显示为河道、河口坝、席状砂砂体与支流间湾泥岩互层( 图2 )。研究区内整体发育三角洲前缘亚相，在沉积微相分析的基础上结合连井相分析，研究区内横向上砂体广泛发育，垂向上砂体相互叠置，河道间发育分流间湾。须家河组沉积环境为三角洲相。具体可分为三角洲平原亚相、三角洲前缘亚相和前三角洲亚相。其中，须二段主要为三角洲前缘沉积，须四段北部为三角洲平原沉积，东部及西部以三角洲前缘沉积为主，须六段为三角洲前缘沉积( 图3 )。

须家河组为特低孔超低渗致密砂岩气藏。通过铸体薄片资料显示，研究区孔隙类型主要以粒间孔为主。对孔隙度及渗透率实测资料进行统计，储层孔隙度在0.14%~7.84%之间，平均为2.10%；渗透率在0.01~15.83 × 10⁻³μm²之间，平均为0.089 × 10⁻³μm²( 图4 、 图5 )。

其中，须二段物性相对较好，孔隙度平均值为2.56%，渗透率平均值为0.0446 × 10⁻³μm²；须四段物性次之，孔隙度平均值为1.97%，渗透率平均值为0.0306 × 10⁻³μm²；须六段物性最差，孔隙度平均值为1.52%，渗透率平均值为0.0343 × 10⁻³μm²；由下而上须家河组储层物性存在变差趋势。

选择研究区内测井质量较高，资料较全的普陆1井作为标准井。基于地层厚度基本相同、岩性无明显变化的特点，采用模式匹配的方法，对多井目标地层测井采取序列的自然伽马、声波、中子、密度等曲线进行一致性检查和校正处理。确保曲线校正后依然满足实际地质特征在横向上的变化规律，保留不同岩性间的相应差异。对多井目标地层测井采集序列的自然伽马、声波、中子、密度等曲线进行一致性检查和校正处理。

在建立理论岩石模型基础上，应用流体和矿物构成以及岩石结构等基础信息，获得有效的岩石弹性属性，并根据对实测的弹性波曲线与合成曲线的对比结果来确定模型参数。首先选取测井地层评价处理后的各矿物含量、孔隙度、饱和度等参数( 表1 )对普光气田12口目标井进行正演。

地震反演技术对于预测储层具有重要意义 [14] ，地震反演是研究地层特征的关键技术，该技术以地震资料为基础，并精确显示出地质特征 [15] 。首先进行井震标定( 图6 )，为后续储层预测奠定基础。再建立低频模型进行边界约束( 图7 )，以得到可信的纵横波阻抗。继而建立纵波阻抗低频模型( 图8 )，反映出井点出低频的变化，也反映出大套地层弹性参数的变化。经过多次调试，确定叠后反演参数( 表2 )，以此为基础进行叠后确定性反演处理，得到纵波阻抗的反演结果。

根据叠后反演结果提取出纵波阻抗井曲线(红)与实际井点的纵波阻抗曲线(蓝)进行对比( 图9 )，发现两者纵向变化趋势一致，阻抗变化范围相同，为1.0e10⁷~1.7e10⁷，说明叠后反演的结果是可信的。

结合叠后连井剖面对比发现反演结果与单井纵波阻抗一致( 图10 )，且反演结果井间阻抗横向变化清楚，说明反演的参数选取合理。

约束稀疏脉冲反演通过调整反射系数序列的稀疏性来拓宽输入地震数据的有效频带宽度，并得到一个弹性参数模型。叠后反演的结果是非常稳定的，通过叠后反演可以进行一些关键参数及子波的测试，检测反演中所适用的低频、子波及参数是否合适，为后续的叠前反演提供参考数据。

研究区的叠后确定性阻抗反演结果稳定，对反演结果与岩性结合分析( 图11 )，在普光5-普光6井区，须六段砂岩储层表现为高纵波阻抗特征，而其他井区及其他层段砂岩的阻抗特征不明显。普光5-普光6井须六段的部分干砂岩为在交会图中为高纵波阻抗，可以与泥岩区分，分析表明，只有普光5、普光6井须六段的部分干砂岩具高纵波阻抗特征，其他井均不具备，所以，在普光5、普光6井区可以利用高纵波阻抗的特征将部分干砂岩的识别出来，而对于其他井区，该方法的不确定性较大。对于有效储层，则表现为相对较低的阻抗特征，与泥岩的阻抗特征叠置，纵波阻抗参数难以识别。

首先建立纵波阻抗、横波阻抗(VPVS)及密度的低频模型( 图12 )。通过调试和选取各项反演参数( 表3 )，以此为基础进行精确的叠前反演。

本区泥岩对应中高VPVS，砂岩对应中低VPVS，叠前同时反演得到的VPVS对岩性有一定的识别能力，但在VPVS叠置范围会有一些误判；而纵波阻抗在砂泥岩储层中是叠置的，无法有效识别岩性。

反演得到的纵波阻抗联井剖面，反演结果与井点阻抗较一致，并且叠前反演的纵波阻抗与叠后反演的纵波阻抗是一致的，说明反演结果是稳定、可信的( 图13 )。

叠前反演VPVS连井剖面显示( 图14 )，砂岩表现为较低VPVS，其截止值为1.75，图中亮色范围为小于1.75的VPVS分布，与井中GR曲线对比表明，低VPVS可以反映厚砂层及大套的砂层组，与岩石物理分析结果是一致的。说明，叠前反演利用VPVS可以对厚砂层及砂层组进行描述及刻画。但是，小于1.75的部分也包括了部分泥岩，所以利用单一的VPVS以1.75为截止进行分析时，不可避免地会将部分泥岩判断为砂岩，带来岩性解释的误判，增加反演解释结果的不确定性。

叠前确定性同时反演处理表明，针对普光地区的储层预测研究，利用叠前反演可以解决纵波阻抗叠置从而不能识别岩性的问题。叠前反演处理得到的VPVS可以有效刻画其纵向及横向的空间展布。

根据叠前反演VPVS及FFP结果砂岩概率体剖面发现( 图15 )，图中VPVS色标红–黄–绿的亮色为VPVS较低区域，蓝色为高VPVS区域，根据岩石物理量板分析，VPVS较低区域基本为砂岩分布区，VPVS较高区域为泥岩展布区。井中曲线的颜色表示不同的岩性，橙黄色为砂岩，灰色为泥岩。下图中的黑色曲线充填为GR，为砂岩分布层段；剖面中的红色为砂岩概率大于70%的展布区域，黄色为60%概率，绿色为50%概率。对比分析可知，对于厚砂岩储层及大套的砂层组，叠前反演结果均有响应。

图15. PG9-PG5-PG8-PL1-QX2联井反演剖面(上图为VPVS剖面，下图为砂岩概率体剖面)

普陆1井须六段1砂组钻井砂岩厚度46米，2砂组厚度33.5米，在反演结果中均有反映( 图16 )。须四段1砂组及3砂组，分别厚达54米及38米，在反演结果中表现为低VPVS及较高的砂岩概率特征，而虽然2砂组的累计厚度达到28米，但是由于是多套单砂层的叠加，造成在地震资料中反映不明显。须2段三套砂组均较厚，厚度达到50米、38米及84米，同时单套砂层也较厚，所以反演结果的识别较好。1砂组及2砂组之间的泥岩隔层也能够准确识别。

过普陆1井的有效储层概率剖面显示( 图17 )，图中井点曲线为岩性曲线，红色为I+II类有效储层，绿色为干砂岩，灰色为泥岩；剖面中的色标表示有效储层的概率，由蓝色–红色表明其概率增加。

井点分析表明须六~须二段的有效储层单层厚度均较薄，在3~10米之间，但是其累积厚度相对较厚，须六段的1砂组累积厚度为18米，2砂组累积厚度为8米；须四段1砂组累积厚度15米，2砂组累积厚度9米，3砂组累积厚度1.5米；须二段1砂组累积厚度20米，2砂组累积厚度13米，3砂组累积厚度23米。所以，须六段1砂组、须四段1砂组、及须二段1、3砂组的累积厚度均达到了地震反演可识别的范围，在有效储层概率剖面上均能够很好识别。但是，从岩石物理模板分析，叠前反演对有效储层有一定的识别能力；但一方面由于有效储层的弹性参数特征与干砂层的弹性参数特征有部分叠置，另一方面，有效储层厚度均较薄，单层厚度均超出地震分辨率的范围，故对有效储层的预测只能反映有效储层的叠合效果，如其累积厚度达到地震分辨能力，在反演结果会存在相应的响应特征，同时，受地震资料的品质等多种因素的制约，精确的有效储层解释存在风险。