一种海域水合物饱和度预测的方法技术

本发明专利技术公开了一种海域水合物饱和度预测的方法，包括以下步骤：1)地震数据预处理，得到观测炮集数据及偏移速度模型；2)基于步骤1)中获得的观测炮集数据，利用地震近道初至提取地震子波；3)基于步骤1)中获得的偏移速度模型，自上而下建立海水及近海底初始速度模型和密度模型；4)利用步骤2)中获得的地震子波以及步骤3)中自上而下建立海水及近海底初始速度模型及密度模型进行全波形反演，得到最终精细速度模型；5)利用步骤4)中获得的精细速度模型识别水合物，并计算水合物饱和度。

A Method for Predicting Hydrate Saturation in Sea Area

The invention discloses a method for predicting hydrate saturation in sea area, which includes the following steps: 1) seismic data preprocessing to obtain observation artillery gather data and migration velocity model; 2) based on observation artillery gather data obtained in step 1, seismic wavelet is extracted from the first arrival of seismic near-track; 3) based on migration velocity model obtained in step 1, the initial velocity of seawater and offshore bottom is established from top to bottom. Degree model and density model; 4) Using the seismic wavelet obtained in step 2 and step 3), the initial velocity model and density model of seawater and offshore seabed are established from top to bottom, and the final fine velocity model is obtained; 5) Using the fine velocity model obtained in step 4, the hydrate is identified and the hydrate saturation is calculated.

【技术实现步骤摘要】
一种海域水合物饱和度预测的方法
本专利技术涉及一种海域水合物饱和度预测的方法，属于油气勘探领域。
技术介绍
天然气水合物是一种由水和甲烷等气体分子在低温高压下形成的似冰的固态物质，由气、液、固三态形成的混合体，具有能源密度高、清洁环保、分布广等特点。据美国能源部预测全球天然气水合物资源量为20万亿桶油当量，是全球煤、石油、天然气中有机碳总量的2倍。所以天然气水合物是21世纪公认的接替煤炭、石油等常规能源的新型绿色能源。近年来世界各国加大了对水合物的勘探开发研究的力度。地震方法是海洋天然气水合物勘探初期的最主要手段。在勘探阶段，利用地震数据识别水合物的位置和估算水合物的饱和度是首要的任务，是设计钻探井位，进而讨论水合物的安全开采方案的基础工作。近些年，国内已经开展了相关方面的研究。运用重磁电技术，景建恩等(琼东南盆地天然气水合物及其成藏模式的海洋可控源电磁研究，地球物理学报，2018)探讨了海洋可控源电磁方法在探测水合物的位置中的应用。张光学等(三维地震与OBS联合勘探揭示的神狐海域含水合物地层声波速度特征，地球物理学报，2014)提出了利用三维地震和OBS联合勘探水合物的声波速度特征。徐华宁等(南中国海神狐海域天然气水合物地震识别及分布特征，地球物理学报，2010)挖掘地震数据和测井数据潜力，针对南中国海神狐海域水合物开展了地震识别工作。此外还有不少从测井出发开始研究水合物的工作，主要集中于水合物饱和度的计算。但是还有三个方面的问题没有得到很好解决。一是水合物的识别还需要进一步精准。在地震剖面上，水合物一般具有比较明显的似海底反射层(BSR)现象。但因为地下介质构造、介质相变等原因也会产生伪BSR现象，为水合物的识别带来困难。二是水合物的研究大多还离不开井，也就是常说的“一孔之见”。水合物的识别及饱和度的计算，要依仗测井资料，但水合物勘探初期测井数量十分有限，即便有井，由于水合物的特殊位置，获得较高质量井资料也比较难。三是一直以来研究重点集中在水合物的勘探，较少关注与水合物伴生的游离气的预测，而是否存在游离气是水合物安全开采的关键因素之一。解决上述问题，还应该脱离井的束缚在地震数据中挖掘潜力。实际上，在深水区普遍存在水合物高速(通常比围岩速度高300-700m/s)，而游离气低速(通常比围岩速度低300-1000m/s)的特征，所以通过速度识别水合物和游离气是一个有效的手段，但由于水合物纵向一般厚度不大、横向不连续、成块状分布，常规层析成像和速度扫描很难刻画出其横向的速度变化，所以需要借助于更高精度的速度建模手段—全波形反演技术。
技术实现思路
针对上述问题，本专利技术的目的是提供一种海域水合物饱和度预测的方法，该方法利用全波形反演技术反演出精准的速度，根据水合物和围岩的速度差异，把水合物和游离气区别开来，再根据水合物速度和饱和度关系定量预测出水合物的饱和度，为水合物的开发井位选择提供可靠的依据。为实现上述目的，本专利技术采用以下技术方案，一种海域水合物饱和度预测的方法，其特征在于，包括以下步骤：1)地震数据预处理，得到观测炮集数据及偏移速度模型；2)基于步骤1)中获得的观测炮集数据，利用地震近道初至提取地震子波；3)基于步骤1)中获得的偏移速度模型，自上而下建立海水及近海底初始速度模型和密度模型，具体过程如下：3.1)确定偏移速度模型中海水层的速度，并确定海水层的密度；3.2)确定偏移速度模型中的精确海底位置；3.3)确定偏移速度模型中海底以下地层的速度和密度；4)利用步骤2)中获得的地震子波以及步骤3)中自上而下建立的海水及近海底初始速度模型及密度模型进行全波形反演，得到最终精细速度模型；5)利用步骤4)中获得的精细速度模型识别水合物，并计算水合物饱和度；依据水合物的速度特征在步骤4)中获得的精细速度模型上识别水合物；将水合物饱和度和有效孔隙度关系式与有效孔隙度与速度的拟合关系式联合起来，计算水合物饱和度。在上述步骤5)中，水合物饱和度和有效孔隙度关系式为：式中，S为水合物饱和度；Rw为含水地层电阻率；Rt为目的层电阻率；φ有效孔隙度；m、n、a1、b1均为经验常数；有效孔隙度与水合物速度的拟合关系式为：φ＝a2v2+b2v+c1式中，a2、b2、c1为拟合系数；v为水合物速度。上述步骤3.1)中，确定偏移速度模型中海水层的速度和海水层的密度，确定海水层的速度，具体过程如下：a)、获得采集时期海水的温度曲线，获取不同深度下的海水温度T；b)、将步骤a)中获得的海水温度T代入乌德公式估算海水速度，获得海水速度曲线图；乌德公式如下：C＝1450+4.21T-0.037T2+1.14(S-35)+0.175P式中，C为海水速度；T为海水温度；S为海水盐度；P为海水压力；确定海水层的密度，具体为：海水的密度取常值，为1.03kg/m3。在上述步骤3.2)中，确定偏移速度模型中的精确海底位置，具体过程如下：i)、基于上述步骤2)获得的地震子波数据、经上述步骤3.1)处理后的具有确定海水层速度和密度的偏移速度模型，采用正演算法获得正演炮集数据；ii)、将步骤i)中获得正演炮集数据中的海底反射旅行时和步骤1)中获得的观测炮集数据的海底反射旅行时进行比对；如果存在较大旅行时差，调整海底位置；iii)、重复步骤i)和ii)，直到步骤i)中获得的正演炮集数据和步骤1)中的观测炮集数据的海底反射旅行时一致，完成偏移速度模型中精确海底位置的确定。在上述步骤3.3)中，确定偏移速度模型中海底以下地层的速度和密度；具体过程如下：①、基于上述步骤2)获得的地震子波数据、经上述步骤3.1)和3.2)处理后的具有确定海水层速度和密度以及精确海底位置的偏移速度模型，利用正演算法获得正演炮集数据；②、比较步骤①中获得的正演炮集数据的海底反射振幅和步骤1)中获得的观测炮集数据的海底反射振幅的大小来确定海底的反射系数；③、基于②获得的反射系数，根据反射系数公式联合Gardner公式确定海底以下地层的速度和密度；反射系数公式如下：式中，R为反射系数；ρ2和V2分别为下部地层密度和速度；ρ1和V1分别为上部地层密度和速度；Gardner公式如下：ρ＝aVc式中，ρ为介质密度；V为介质速度；a、c均为常数。在上述步骤3.3)中，在确定海底以下地层的速度和密度时，还应分别考虑硬海底条件和软海底条件，对于硬海底条件下，上述步骤3.3)中的第③步中，Gardner公式中的常数a为0.31，c为0.25；对于软海底情况下，上述步骤3.3)中的第③步中，Gardner公式中的常数a和c由拟合得到。在上述步骤1)中，地震数据预处理包括地震炮集数据非线性和线性噪音压制、滤波和偏移，同时地震数据通过偏移处理获得偏移速度模型。本专利技术采用以上技术方案，其具有如下优点：1、本专利技术通过对地震数据进行预处理，根据预处理得到的观测炮集数据及偏移速度模型自上而下建立海水及近海底的初始速度模型和密度模型，再基于初始速度密度模型依靠全波形反演技术建立高精度的速度模型；根据水合物的速度特征在速度模型上识别水合物，并将水合物和游离气区别开来，为水合物的开发井位选择提供可靠的依据。2、本专利技术将速度和有效孔隙度的拟合关系式、有效孔隙度和水合物饱和度的关系式联合起来，计算水合物的饱和度，实本文档来自技高网...

【技术保护点】
1.一种海域水合物饱和度预测的方法，其特征在于，包括以下步骤：1)地震数据预处理，得到观测炮集数据及偏移速度模型；2)基于步骤1)中获得的观测炮集数据，利用地震近道初至提取地震子波；3)基于步骤1)中获得的偏移速度模型，自上而下建立海水及近海底初始速度模型和密度模型，具体过程如下：3.1)确定偏移速度模型中海水层的速度，并确定海水层的密度；3.2)确定偏移速度模型中的精确海底位置；3.3)确定偏移速度模型中海底以下地层的速度和密度；4)利用步骤2)中获得的地震子波以及步骤3)中自上而下建立的海水及近海底初始速度模型及密度模型进行全波形反演，得到最终精细速度模型；5)利用步骤4)中获得的精细速度模型识别水合物，并计算水合物饱和度；依据水合物的速度特征在步骤4)中获得的精细速度模型上识别水合物；将水合物饱和度和有效孔隙度关系式与有效孔隙度与速度的拟合关系式联合起来，计算水合物饱和度。

【技术特征摘要】
1.一种海域水合物饱和度预测的方法，其特征在于，包括以下步骤：1)地震数据预处理，得到观测炮集数据及偏移速度模型；2)基于步骤1)中获得的观测炮集数据，利用地震近道初至提取地震子波；3)基于步骤1)中获得的偏移速度模型，自上而下建立海水及近海底初始速度模型和密度模型，具体过程如下：3.1)确定偏移速度模型中海水层的速度，并确定海水层的密度；3.2)确定偏移速度模型中的精确海底位置；3.3)确定偏移速度模型中海底以下地层的速度和密度；4)利用步骤2)中获得的地震子波以及步骤3)中自上而下建立的海水及近海底初始速度模型及密度模型进行全波形反演，得到最终精细速度模型；5)利用步骤4)中获得的精细速度模型识别水合物，并计算水合物饱和度；依据水合物的速度特征在步骤4)中获得的精细速度模型上识别水合物；将水合物饱和度和有效孔隙度关系式与有效孔隙度与速度的拟合关系式联合起来，计算水合物饱和度。2.如权利要求1所述的一种海域水合物饱和度预测的方法，其特征在于，在上述步骤5)中，水合物饱和度和有效孔隙度关系式为：式中，S为水合物饱和度；Rw为含水地层电阻率；Rt为目的层电阻率；φ有效孔隙度；m、n、a1、b1均为经验常数；有效孔隙度与水合物速度的拟合关系式为：φ＝a2v2+b2v+c1式中，a2、b2、c1为拟合系数；v为水合物速度。3.如权利要求1所述的一种海域水合物饱和度预测的方法，其特征在于，上述步骤3.1)中，确定偏移速度模型中海水层的速度和海水层的密度，确定海水层的速度的具体过程如下：a)、获得采集时期海水的温度曲线，获取不同深度下的海水温度T；b)、将步骤a)中获得的海水温度T代入乌德公式估算海水速度，获得海水速度曲线图；乌德公式如下：C＝1450+4.21T-0.037T2+1.14(S-35)+0.175P式中，C为海水速度；T为海水温度；S为海水盐度；P为海水压力；确定海水层的密度，具体为：海水的密度取常值，为1.03kg/m3。4.如权利要求1所述的一种海域水合物饱和度预测的方法，其特征在于，在上述步骤3.2)...

【专利技术属性】
技术研发人员：朱振宇，丁继才，杜向东，于水，张金淼，翁斌，姜秀娣，张益明，张云鹏，王小六，陈剑军，薛东川，
申请(专利权)人：中国海洋石油集团有限公司，中海油研究总院有限责任公司，
类型：发明
国别省市：北京,11