储层物性参数反演方法、装置、设备及存储介质.pdf

(19)国家知识产权局(12)发明专利申请(10)申请公布号CN116027398A(43)申请公布日2023.04.28(21)申请号202111241483.0(22)申请日2021.10.25(71)申请人中国石油化工股份有限公司地址100728北京市朝阳区朝阳门北大街22号申请人中国石油化工股份有限公司石油物探技术研究院(72)发明人袁媛刘卫华索重辉韩德超(74)专利代理机构北京聿宏知识产权代理有限公司11372专利代理师陈超德吴昊(51)Int.Cl.G01V1/28(2006.01)G01V1/30(2006.01)G01V1/36(2006.01)权利要求书3页说明书22页附图3页(54)发明名称储层物性参数反演方法、装置、设备及存储介质(57)摘要本申请提供的一种储层物性参数反演方法、装置、设备及存储介质，包括：分析浊积砂岩储层的基本岩石物理特征，得到工区浊积砂岩的性质；根据研究工区浊积砂岩的性质，对软砂岩模型进行校正；根据Gassmann流体替换方程和校正后的软砂岩模型，求得纵波速度、横波速度和密度的表达式，构成岩石物理模型；根据建立的岩石物理模型进行储层物性参数敏感度分析；基于敏感度分析结果对储层物性参数进行反演；将反演结果与实测曲线对比分析，验证结果的准确性，使反演结果更稳定且提高了计算效率。CN116027398ACN116027398A权利要求书1/3页1.一种储层物性参数反演方法，其特征在于，包括：S1：分析浊积砂岩储层的基本岩石物理特征，得到工区浊积砂岩的性质；S2：根据工区浊积砂岩的性质，对软砂岩模型进行校正；S3：根据Gassmann流体替换方程和校正后的软砂岩模型，求得纵波速度、横波速度和密度的表达式，构成岩石物理模型；S4：根据建立的岩石物理模型进行储层物性参数敏感度分析；S5：基于敏感度分析结果对储层物性参数进行反演；S6：将反演结果与实测曲线对比分析，验证结果的准确性。2.根据权利要求1所述的方法，其特征在于，根据建立的岩石物理模型进行储层物性参数敏感度分析之前，所述方法还包括：以工区内测井数据为基础对建立的岩石物理模型进行验证；所述验证的具体方法包括：通过计算将储层物性参数转换为纵波速度、横波速度和密度，再与实际测井曲线进行对比分析，如果曲线拟合较好，说明建立的岩石物理模型适用于该工区，反之则需要对岩石物理模型的参数进行调整，重新进行检验，直到模型正演得到的纵波速度、横波速度和密度与实际测井的纵波速度、横波速度和密度达到较好的拟合。3.根据权利要求1所述的方法，其特征在于，所述软砂岩模型的具体公式包括：式中，Kdry和μdry分别表示干岩石骨架的等效体积模量和剪切模量；φC为临界孔隙度，n为每个颗粒接触面的个数，R为颗粒的等效半径；Fn和Fτ分别表示干燥条件下颗粒在一定围压P作用下产生的法向和切向刚度,μ表示颗粒集合体的剪切模量，υ表示泊松比，为接触面的平均曲率半径；KHM：等径颗粒堆积物的等效体积模量μHM：等径颗粒堆积物的等效剪切模量φ：孔隙度Kg：固相颗粒的等效体积模量μg：固相颗粒的等效剪切模量。2CN116027398A权利要求书2/3页4.根据权利要求3所述的方法，其特征在于，所述对软砂岩模型进行校正的具体方法包括：引入了切向刚度校正因子K，接触颗粒面的局部平均曲率半径则被定义为：式中R1和R2分别为两个接触颗粒的半径；式中K的取值范围在0～1之间，反映了颗粒接触面的特征和松弛作用的影响；的值大于或小于1，当岩石是等径大小的球状颗粒时该值等于1。5.根据权利要求1所述的方法，其特征在于，所述根据Gassmann流体替换方程和校正后的软砂岩模型，求得纵波速度、横波速度和密度的表达式的具体方法包括：根据Gassmann流体替换方程得到饱含流体状态下岩石的体积模量和剪切模量，从而求得纵波、横波速度以及密度的表达式，μsat＝μdryρ＝φρfl+(1‑φ)ρg(10)式中，Kma：岩石组成基质的体积模量；Kfl：介质孔隙流体的体积模量；ρfl：介质孔隙流体密度；ρg：基质矿物颗粒等效密度；ksat：饱含流体状态下岩石的体积模量；μsat：饱含流体状态下岩石的剪切模量；vp：纵波速度；vs：横波速度；ρ：密度。3CN116027398A权利要求书3/3页6.根据权利要求1所述的方法，其特征在于，所述根据建立的岩石物理模型进行储层物性参数敏感度分析的具体方法包括：通过校正后的软砂岩模型和Gassmann流体替换方程将储层物性参数转换为岩石弹性参数vp,vs,ρ，再由精确Zoeppritz方程计算得到反射系数，进而分析储层物性参数变化对反射系数的影响。7.根据权利要求6所述的方法，其特征在于，基于敏感度分析结果对孔隙度、地层有效压力和含水饱和度参数进行反演的具体方法包括：利用最小二乘