CN108647461A - 盐底劈构造带水合物稳定底界的二维数值模拟方法及系统 - Google Patents
盐底劈构造带水合物稳定底界的二维数值模拟方法及系统 Download PDFInfo
- Publication number
- CN108647461A CN108647461A CN201810476185.1A CN201810476185A CN108647461A CN 108647461 A CN108647461 A CN 108647461A CN 201810476185 A CN201810476185 A CN 201810476185A CN 108647461 A CN108647461 A CN 108647461A
- Authority
- CN
- China
- Prior art keywords
- bsr
- hydrate
- salt
- split
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于天然气水合物勘探技术领域,公开了一种盐底劈构造带水合物稳定底界的二维数值模拟方法及系统,分析地温梯度、盐和沉积物不同的热导率以及盐底劈构造带的盐度高异常对水合物稳定底界(地震资料上表现为BSR)的影响,并将其反映在水合物相平衡条件中对BSR进行数值模拟;系统包括BSR特征选择模块、双程走时确定模块、天然气水合物的相平衡稳定曲线选择模块、BSR发育深度处的压力值确定模块、BSR的深度值计算模块、地震剖面上海底各点的温度值确定模块、地震剖面的温度场分布计算模块、水合物稳定底界移向模块、对比分析模块。本发明可用于任何边界明显的盐底劈构造带来预测其水合物稳定底界的位置。
Description
技术领域
本发明属于天然气水合物勘探技术领域,尤其涉及一种盐底劈构造带水合物稳定底界的二维数值模拟方法及系统。具体地说,涉及一种利用地震解释和数值模拟方法,来描述和预测盐底劈构造带水合物稳定底界的方法。
背景技术
目前,业内常用的现有技术是这样的:
海域天然气水合物稳定区的底界(BHSZ:Base of Gas hydrate stability zone)在地震资料上一般表现为似海底反射层(Bottom simulating reflection:BSR),代表了含水合物地层与底部含游离气地层之间的相界面,是天然气水合物重要的地震反射标志。一般认为海底到BHSZ之间的区域为水合物稳定区(GHSZ:Gas hydrate stability zone)。由于天然气水合物稳定性较差,主要受温度、压力、盐度和气体组分等因素影响,因此水合物稳定底界BHSZ常受沉积和构造作用、海底热流等事件的影响而发生垂向迁移,在地震资料上表现为BSR向上或向下迁移。
对水合物相平衡曲线的分析表明,地温梯度、海底的水深和温度、气体组分和孔隙水盐度都会对BSR深度造成影响,但其中地温梯度和海底温度对其影响最大。
综上所述,现有技术存在的问题是:
目前,对盐底劈构造带水合物稳定底界BSR的分析通常是定性的,认为BSR受盐底劈构造带热流的影响而发生向上迁移,即BSR变浅。另外,孔隙水盐度的增大会抑制水合物的形成,从而使水合物相平衡曲线向左迁移,也会造成BSR的向上迁移。目前针对盐底辟构造的热效应分析也有一定进展,很多学者对盐底劈构造带的地温场进行了定量表征。由于盐的导热率较高,为普通沉积岩的2-3倍,因此盐底劈构造带会影响地温场,造成盐上地层的高温异常和盐下地层的低温异常。
解决上述技术问题的难度和意义:
本发明针对盐底劈构造带水合物稳定底界BSR缺少定量深入分析的现状,将盐底辟构造带具体的、量化的热效应和盐度效应体现到对BSR的数值模拟中去,充分考虑盐岩和普通沉积物的导热率差别,以及盐底劈构造带的盐度异常对水合物稳定条件的影响,利用数值模拟软件MATLAB来实现对不同控制条件下BSR位置的模拟。盐底辟构造在海洋环境较为常见,本发明提供了一种定量分析盐底辟构造带的热效应和盐度效应对水合物稳定底界的影响,为分析其他盐底辟构造地区的水合物稳定底界的发育位置、水合物稳定底界的主要控制因素等提供了依据,为将来盐底辟构造带水合物藏的开发和环境影响评价等工作打下基础。
发明内容
针对现有技术及存在的问题,本发明提供了一种盐底劈构造带水合物稳定底界的二维数值模拟方法及系统。本发明针对盐底劈构造带的天然气水合物稳定底界的分布位置,提供一种二维数值模拟方法,来预测盐底劈构造带的水合物稳定底界;通过和观测的BSR进行拟合,从而认识盐底劈构造的热效应和盐度效应对水合物稳定底界的影响。
本发明是这样实现的,一种盐底劈构造带水合物稳定底界的二维数值模拟方法,包括以下步骤:
步骤①:在盐底劈构造带,选择水合物发育且地震资料上BSR特征明显的地震剖面;
步骤②:确定地震剖面上海底界面的双程走时Tsb、观测BSR的双程走时Tbsr和盐底劈构造带顶部边界的双程走时Tsalt;
步骤③:根据分析区的气体组分信息,选择合适的水合物相平衡曲线。本发明选择Moridis(2003)的纯甲烷水合物的相平衡稳定曲线,相平衡公式如下:
In(P)=a+bT+cT2+dT3+fT4+gT5 (1)
式中,P和T分别为甲烷水合物的稳定压力条件和稳定温度条件,a、b、c、d、f、g为经验常数,分别为a=-1.94138504464560×105,b=3.31018213397926×103,c=-2.25540264493806×101,d=7.67559117787059×10-2,f=-1.30465829788791×10-4,g=8.8606531668757×10-8;
步骤④:由于水合物稳定底界的深度较浅,故假设地震剖面上BSR发育深度处的压力值为静水压力,压力可通过以下公式求取:
Pbsr=ρsw g Hbsr (2)
式中,ρsw为海水密度,为1028kg/m3,g是重力加速度,为9.81m/s2,Hbsr是BSR的深度值;
步骤⑤:BSR的深度值可通过下面公式求得:
Hbsr=Vsw×Tbsr/2 (3)
式中,Vsw为海水速度,为1500m/s;Tbsr为BSR位置的双程走时,单位是s,可通过地震剖面读取;
步骤⑥:通过World Ocean Database等数据库查询,确定地震剖面上海底各点的温度值Tsb;
步骤⑦:通过二维稳态热传导方程,计算地震剖面的温度场分布:
式中,T为温度(℃),x为横向距离(km),z为垂向距离(km),kx为横向导热率(W m- 1K-1),kz为垂向导热率(W m-1K-1)。一般认为沉积物各向同性且均质性较强,因此kx=kz,但沉积物和盐的热导率差别较大,分别设置为2.5W m-1K-1和5.9W m-1K-1。
步骤⑧:盐底劈构造带的盐度异常会抑制水合物的形成,将水合物稳定底界移向较低的温度条件,如下:
式中,m为盐度对水合物稳定条件的影响参数,计算公式如下:
式中,Sw为盐度值;
步骤⑨:通过设置不同的地温梯度,可利用公式(1)(2)(5)来模拟对应的BSR深度,将模拟的BSR位置与地震剖面上的观测BSR进行对比,拟合程度较高的代表了较符合实际的地温梯度条件。
进一步,步骤①选择盐底劈构造带BSR特征明显的地震剖面是为了验证该模拟方法的可靠性,模拟BSR与观测BSR拟合程度好则说明设定的地温梯度、盐度和热导率等参数条件合理,可用于其他BSR特征不明显的盐底劈构造带的BSR深度预测。
进一步,步骤②确定海底、观测BSR和盐底劈构造带顶部边界双程走时,分别是为了在MATLAB软件中模拟时可以设定边界、绘制出观测BSR的位置以及在进行二维稳态热传导公式计算时针对盐底劈构造的分布设定不同于周围沉积物的热导率。
进一步,步骤⑦中的二维稳态热传导方程的上部边界条件是将海底温度设定为海水在海底深度处的温度值,下部边界是通过将地温梯度在在海平面以下5km的位置设置为均值,即区域地温梯度值,可通过设定不同的地温梯度值得到不同条件下的模拟BSR位置。
进一步,步骤⑧考虑到了盐底劈构造带的盐度高异常会抑制水合物的形成,造成水合物稳定条件向较低的温度偏移,本方法未将盐度设置为较高的35‰。
本发明的另一目的在于提供一种实现所述盐底劈构造带水合物稳定底界的二维数值模拟方法的计算机程序。
本发明的另一目的在于提供一种实现所述盐底劈构造带水合物稳定底界的二维数值模拟方法的信息数据处理终端。
本发明的另一目的在于提供一种计算机可读存储介质,包括指令,当其在计算机上运行时,使得计算机执行所述盐底劈构造带水合物稳定底界的二维数值模拟方法。
本发明的另一目的在于提供一种所述盐底劈构造带水合物稳定底界的二维数值模拟方法的盐底劈构造带水合物稳定底界的二维数值模拟系统,其特征在于,所述盐底劈构造带水合物稳定底界的二维数值模拟系统包括:
BSR特征选择模块,用于在水合物发育的盐底劈构造带,选择BSR特征明显的地震剖面;
双程走时确定模块,用于确定地震剖面上海底界面的双程走时Tsb、观测BSR的双程走时Tbsr和盐底劈构造带顶部边界的双程走时Tsalt;
天然气水合物的相平衡稳定曲线选择模块,用于选择天然气水合物的相平衡稳定曲线,纯甲烷水合物相平衡公式如下:
In(P)=a+bT+cT2+dT3+fT4+gT5
式中,P和T分别为甲烷水合物的稳定压力条件和稳定温度条件,a、b、c、d、f、g为经验常数,分别为a=-1.94138504464560×105,b=3.31018213397926×103,c=-2.25540264493806×101,d=7.67559117787059×10-2,f=-1.30465829788791×10-4,g=8.8606531668757×10-8;
BSR发育深度处的压力值确定模块,地震剖面上BSR发育深度处的压力值为静水压力,压力通过以下公式求取:
Pbsr=ρsw g Hbsr
式中,ρsw为海水密度,为1028kg/m3,g是重力加速度,为9.81m/s2,Hbsr是BSR的深度值;
BSR的深度值计算模块,用于BSR的深度值通过下面公式求得:
Hbsr=Vsw×Tbsr/2
式中,Vsw为海水速度,为1500m/s;Tbsr为BSR位置的双程走时,单位是s,通过地震剖面读取;
地震剖面上海底各点的温度值确定模块,通过数据库查询,确定地震剖面上海底各点的温度值Tsb;
地震剖面的温度场分布计算模块,通过二维稳态热传导方程,计算地震剖面的温度场分布:
式中,T为温度(℃),x为横向距离(km),z为垂向距离(km),kx为横向导热率(W m- 1K-1),kz为垂向导热率(W m-1K-1)。一般认为沉积物各向同性且均质性较强,因此kx=kz,但沉积物和盐的热导率差别较大,分别设置为2.5W m-1K-1和5.9W m-1K-1;
水合物稳定底界移向模块,将水合物稳定底界移向低的温度条件,
式中,m为盐度对水合物稳定条件的影响参数,计算公式如下:
式中,Sw为盐度值;
对比分析模块,通过设置不同的地温梯度,模拟对应的BSR深度,将模拟的BSR位置与地震剖面上的观测BSR进行对比。
本发明的另一目的在于提供一种搭载有所述盐底劈构造带水合物稳定底界的二维数值模拟系统的信息数据处理终端。
综上所述,本发明的优点及积极效果为:
针对盐底劈构造带水合物稳定底界BSR缺少定量深入分析的现状,将盐底辟构造的具体、量化的热效应和盐度效应体现到对BSR的数值模拟中去,充分考虑盐岩和普通沉积物的导热率差别,绘制受盐岩影响的温度等值线,结合盐底劈构造带的盐度异常对水合物稳定条件的影响,利用数值模拟软件MATLAB来实现对不同控制条件下BSR位置的模拟。本发明能够用于模拟其他盐底劈构造地区水合物稳定底界的发育位置,为分析水合物稳定底界的主控因素提供依据,也为将来盐底辟构造带水合物藏的开发和环境影响评价等工作打下基础。
本发明通过选取合适的水合物相平衡公式,对BSR进行数值,其目的是通过设定不同的参数值(通常是地温梯度)来模拟不同参数条件下的BSR位置,再将模拟BSR位置与观察到的BSR位置进行比较,则拟合程度较高的参数值代表了较符合实际的参数条件。在此参数基础上,可对BSR不存在的地区进行数值模拟,分析这些地区潜在的水合物稳定底界的位置。
本发明可在其他水合物发育但BSR地震反射特征不明显的盐底劈构造区域进行应用,预测盐底劈构造带附近的BSR界面;
本发明在对BSR进行数值模拟时考虑了盐和沉积物的不同导热率的影响,以及盐底劈构造带的盐度高异常的影响,模拟条件参数考虑较全面,因此,模拟BSR和观测BSR拟合程度较高的条件可能较符合实际的地温梯度、导热率和盐度条件。
附图说明
图1是本发明实施例提供的盐底劈构造带水合物稳定底界的二维数值模拟方法流程图。
图2是本发明实施例提供的天然气水合物相平衡曲线示意图,反映了天然气水合物稳定底界BSR随盐度和热流(地温梯度)的增大而发生的向上偏移。
图3是本发明实施例提供的西非地区盐底劈构造带的水合物BSR界面的模拟结果,当地温梯度设置为32℃ km-1,盐和沉积物的热导率分别为5.9W m-1K-1和2.5W m-1K-1,盐度设置为35‰时,模拟BSR和观测BSR拟合效果最好。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明综合考虑了地温梯度、盐和沉积物不同的热导率以及盐底劈构造带的盐度高异常对水合物稳定底界(地震资料上表现为BSR)的影响,并将其反映在水合物相平衡条件中对BSR进行数值模拟。模拟BSR和观测BSR拟合程度高则代表了所选地温梯度、热导率、盐度等参数条件较符合实际。确定了该模拟方法可行后,还可在其他BSR特征不明显的盐底劈构造带进行模拟,预测BSR位置。总之,该方法适用性较强,可用于任何边界明显的盐底劈构造带来预测其水合物稳定底界的位置。
如图1、图2、图3所示,本发明实施例提供的盐底劈构造带水合物稳定底界的二维数值模拟方法。步骤如下:
步骤①:在盐底劈构造带,选择水合物发育且地震资料上BSR特征明显的地震剖面;
步骤②:确定地震剖面上海底界面的双程走时Tsb、观测BSR的双程走时Tbsr和盐底劈构造带顶部边界的双程走时Tsalt;
步骤③:根据分析区的气体组分信息,选择合适的水合物相平衡曲线。本发明选择Moridis(2003)的纯甲烷水合物的相平衡稳定曲线,相平衡公式如下:
In(P)=a+bT+cT2+dT3+fT4+gT5 (1)
式中,P和T分别为甲烷水合物的稳定压力条件和稳定温度条件,a、b、c、d、f、g为经验常数,分别为a=-1.94138504464560×105,b=3.31018213397926×103,c=-2.25540264493806×101,d=7.67559117787059×10-2,f=-1.30465829788791×10-4,g=8.8606531668757×10-8;
步骤④:由于水合物稳定底界的深度较浅,故假设地震剖面上BSR发育深度处的压力值为静水压力,压力可通过以下公式求取:
Pbsr=ρsw g Hbsr (2)
式中,ρsw为海水密度,为1028kg/m3,g是重力加速度,为9.81m/s2,Hbsr是BSR的深度值;
步骤⑤:BSR的深度值可通过下面得
Hbsr=Vsw×Tbsr/2 (3)
式中,Vsw为海水速度,为1500m/s;Tbsr为BSR位置的双程走时,单位是s,可通过地震剖面读取;
步骤⑥:通过World Ocean Database等数据库查询,确定地震剖面上海底各点的温度值Tsb;
步骤⑦:通过二维稳态热传导方程,计算地震剖面的温度场分布:
式中,T为温度(℃),x为横向距离(km),z为垂向距离(km),kx为横向导热率(W m- 1K-1),kz为垂向导热率(W m-1K-1)。一般认为沉积物各向同性且均质性较强,因此kx=kz,但沉积物和盐的热导率差别较大,分别设置为2.5W m-1K-1和5.9W m-1K-1。
步骤⑧:盐底劈构造带的盐度异常会抑制水合物的形成,将水合物稳定底界移向较低的温度条件,如下:
式中,m为盐度对水合物稳定条件的影响参数,计算公式如下:
式中,Sw为盐度值;
步骤⑨:通过设置不同的地温梯度,可利用公式(1)(2)(5)来模拟对应的BSR深度,将模拟的BSR位置与地震剖面上的观测BSR进行对比,拟合程度较高的代表了较符合实际的地温梯度条件。
本发明实施例提供的盐底劈构造带水合物稳定底界的二维数值模拟系统包括:
BSR特征选择模块,用于在盐底劈构造带,选择水合物发育且地震资料上BSR特征明显的地震剖面;
双程走时确定模块,用于确定地震剖面上海底界面的双程走时Tsb、观测BSR的双程走时Tbsr和盐底劈构造带顶部边界的双程走时Tsalt;
天然气水合物的相平衡稳定曲线选择模块,用于选择纯天然气合物的相平衡稳定曲线,纯甲烷水合物的相平衡公式如下:
In(P)=a+bT+cT2+dT3+fT4+gT5
式中,P和T分别为甲烷水合物的稳定压力条件和稳定温度条件,a、b、c、d、f、g为经验常数,分别为a=-1.94138504464560×105,b=3.31018213397926×103,c=-2.25540264493806×101,d=7.67559117787059×10-2,f=-1.30465829788791×10-4,g=8.8606531668757×10-8;
BSR发育深度处的压力值确定模块,地震剖面上BSR发育深度处的压力值为静水压力,压力通过以下公式求取:
Pbsr=ρsw g Hbsr
式中,ρsw为海水密度,为1028kg/m3,g是重力加速度,为9.81m/s2,Hbsr是BSR的深度值;
BSR的深度值计算模块,用于BSR的深度值通过下面公式求得:
Hbsr=Vsw×Tbsr/2
式中,Vsw为海水速度,为1500m/s;Tbsr为BSR位置的双程走时,单位是s,通过地震剖面读取;
地震剖面上海底各点的温度值确定模块,通过数据库查询,确定地震剖面上海底各点的温度值Tsb;
地震剖面的温度场分布计算模块,通过二维稳态热传导方程,计算地震剖面的温度场分布:
式中,T为温度(℃),x为横向距离(km),z为垂向距离(km),kx为横向导热率(W m- 1K-1),kz为垂向导热率(W m-1K-1)。一般认为沉积物各向同性且均质性较强,因此kx=kz,但沉积物和盐的热导率差别较大,分别设置为2.5W m-1K-1和5.9W m-1K-1;
水合物稳定底界移向模块,将水合物稳定底界移向低的温度条件,
式中,m为盐度对水合物稳定条件的影响参数,计算公式如下:
式中,Sw为盐度值;
对比分析模块,通过设置不同的地温梯度,模拟对应的BSR深度,将模拟的BSR位置与地震剖面上的观测BSR进行对比。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种盐底劈构造带水合物稳定底界的二维数值模拟方法,其特征在于,所述盐底劈构造带水合物稳定底界的二维数值模拟方法包括:
综合分析地温梯度、盐和沉积物不同的热导率以及盐底劈构造带的盐度高异常对水合物稳定底界的影响,通过选取合适的水合物相平衡公式,对BSR进行数值模拟;通过设定不同的地温梯度、热导率、盐度参数值模拟不同参数条件下的BSR位置;再将模拟BSR位置与观察到的BSR位置进行拟合程度的比较。
2.如权利要求1所述的盐底劈构造带水合物稳定底界的二维数值模拟方法,其特征在于,所述盐底劈构造带水合物稳定底界的二维数值模拟方法具体包括:
步骤一:在盐底劈构造带,选择水合物发育且地震资料上BSR特征明显的地震剖面;
步骤二:确定地震剖面上海底界面的双程走时Tsb、观测BSR的双程走时Tbsr和盐底劈构造带顶部边界的双程走时Tsalt;
步骤三:根据分析区天然气水合物的气体组分信息,选择合适的水合物相平衡曲线。如纯甲烷水合物的相平衡公式如下:
In(P)=a+bT+cT2+dT3+fT4+gT5
式中,P和T分别为甲烷水合物的稳定压力条件和稳定温度条件,a、b、c、d、f、g为经验常数,分别为a=-1.94138504464560×105,b=3.31018213397926×103,c=-2.25540264493806×101,d=7.67559117787059×10-2,f=-1.30465829788791×10-4,g=8.8606531668757×10-8;
步骤四:地震剖面上BSR发育深度处的压力值为静水压力,压力通过以下公式求取:
Pbsr=ρsw g Hbsr
式中,ρsw为海水密度,为1028kg/m3,g是重力加速度,为9.81m/s2,Hbsr是BSR的深度值;
步骤五:BSR的深度值通过下面公式求得:
Hbsr=Vsw×Tbsr/2
式中,Vsw为海水速度,为1500m/s;Tbsr为BSR位置的双程走时,单位是s,通过地震剖面读取;
步骤六:通过数据库查询,确定地震剖面上海底各点的温度值Tsb;
步骤七:通过二维稳态热传导方程,计算地震剖面的温度场分布:
式中,T为温度℃,x为横向距离km,z为垂向距离km,kx为横向导热率W m-1K-1,kz为垂向导热率W m-1K-1;一般认为沉积物各向同性且均质性较强,因此kx=kz,但沉积物和盐的热导率差别较大,分别设置为2.5W m-1K-1和5.9W m-1K-1;
步骤八:将水合物稳定底界移向低的温度条件,
式中,m为盐度对水合物稳定条件的影响参数,计算公式如下:
式中,Sw为盐度值;
步骤九:通过设置不同的地温梯度,模拟对应的BSR深度,将模拟的BSR位置与地震剖面上的观测BSR进行对比。
3.一种实现权利要求1~2任意一项所述盐底劈构造带水合物稳定底界的二维数值模拟方法的计算机程序。
4.一种实现权利要求1~2任意一项所述盐底劈构造带水合物稳定底界的二维数值模拟方法的信息数据处理终端。
5.一种计算机可读存储介质,包括指令,当其在计算机上运行时,使得计算机执行如权利要求1-2任意一项所述盐底劈构造带水合物稳定底界的二维数值模拟方法。
6.一种如权利要求1所述盐底劈构造带水合物稳定底界的二维数值模拟方法的盐底劈构造带水合物稳定底界的二维数值模拟系统,其特征在于,所述盐底劈构造带水合物稳定底界的二维数值模拟系统包括:
BSR特征选择模块,用于在水合物发育的盐底劈构造带,选择BSR特征明显的地震剖面;
双程走时确定模块,用于确定地震剖面上海底界面的双程走时Tsb、观测BSR的双程走时Tbsr和盐底劈构造带顶部边界的双程走时Tsalt;
天然气水合物的相平衡稳定曲线选择模块,用于选择水合物的相平衡稳定曲线,纯甲烷水合物的相平衡公式如下:
In(P)=a+bT+cT2+dT3+fT4+gT5
式中,P和T分别为甲烷水合物的稳定压力条件和稳定温度条件,a、b、c、d、f、g为经验常数,分别为a=-1.94138504464560×105,b=3.31018213397926×103,c=-2.25540264493806×101,d=7.67559117787059×10-2,f=-1.30465829788791×10-4,g=8.8606531668757×10-8;
BSR发育深度处的压力值确定模块,地震剖面上BSR发育深度处的压力值为静水压力,压力通过以下公式求取:
Pbsr=ρsw g Hbsr
式中,ρsw为海水密度,为1028kg/m3,g是重力加速度,为9.81m/s2,Hbsr是BSR的深度值;
BSR的深度值计算模块,用于BSR的深度值通过下面公式求得:
Hbsr=Vsw×Tbsr/2
式中,Vsw为海水速度,为1500m/s;Tbsr为BSR位置的双程走时,单位是s,通过地震剖面读取;
地震剖面上海底各点的温度值确定模块,通过数据库查询,确定地震剖面上海底各点的温度值Tsb;
地震剖面的温度场分布计算模块,通过二维稳态热传导方程,计算地震剖面的温度场分布:
式中,T为温度℃,x为横向距离km,z为垂向距离km,kx为横向导热率W m-1K-1,kz为垂向导热率W m-1K-1;
水合物稳定底界移向模块,将水合物稳定底界移向低的温度条件,
式中,m为盐度对水合物稳定条件的影响参数,计算公式如下:
式中,Sw为盐度值;
对比分析模块,通过设置不同的地温梯度,模拟对应的BSR深度,将模拟的BSR位置与地震剖面上的观测BSR进行对比。
7.一种搭载有权利要求6所述盐底劈构造带水合物稳定底界的二维数值模拟系统的信息数据处理终端。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810476185.1A CN108647461B (zh) | 2018-05-17 | 2018-05-17 | 盐底劈构造带水合物稳定底界的二维数值模拟方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810476185.1A CN108647461B (zh) | 2018-05-17 | 2018-05-17 | 盐底劈构造带水合物稳定底界的二维数值模拟方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108647461A true CN108647461A (zh) | 2018-10-12 |
CN108647461B CN108647461B (zh) | 2022-04-19 |
Family
ID=63756764
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810476185.1A Active CN108647461B (zh) | 2018-05-17 | 2018-05-17 | 盐底劈构造带水合物稳定底界的二维数值模拟方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108647461B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109598019A (zh) * | 2018-10-26 | 2019-04-09 | 中国石油天然气股份有限公司 | 超低渗透油气藏微纳米孔隙泡点压力数值模拟方法及装置 |
CN110456024A (zh) * | 2019-08-23 | 2019-11-15 | 中国石油大学(华东) | 一种分析天然气水合物稳定区边界处碳循环过程的方法及系统 |
CN110954960A (zh) * | 2019-11-20 | 2020-04-03 | 中国石油大学(华东) | 一种海底冷泉地震响应数值模拟方法及系统 |
CN117408186A (zh) * | 2023-12-14 | 2024-01-16 | 青岛海洋地质研究所 | 地震剖面中天然气水合物似海底反射层的判断方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007290986A (ja) * | 2006-04-21 | 2007-11-08 | National Institute Of Advanced Industrial & Technology | ハイドレートの3相平衡条件の予測方法 |
WO2014078624A2 (en) * | 2012-11-16 | 2014-05-22 | Chevron U.S.A. Inc. | Methods and systems for determining minimum porosity for presence of clathrates in sediment |
CN105224704A (zh) * | 2014-06-25 | 2016-01-06 | 中国科学院海洋研究所 | 一种海底不稳定性的评价方法 |
CN107861158A (zh) * | 2017-11-07 | 2018-03-30 | 中国石油大学(华东) | 一种确定天然气水合物bsr界面与海底相交位置的研究方法 |
-
2018
- 2018-05-17 CN CN201810476185.1A patent/CN108647461B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007290986A (ja) * | 2006-04-21 | 2007-11-08 | National Institute Of Advanced Industrial & Technology | ハイドレートの3相平衡条件の予測方法 |
WO2014078624A2 (en) * | 2012-11-16 | 2014-05-22 | Chevron U.S.A. Inc. | Methods and systems for determining minimum porosity for presence of clathrates in sediment |
CN105224704A (zh) * | 2014-06-25 | 2016-01-06 | 中国科学院海洋研究所 | 一种海底不稳定性的评价方法 |
CN107861158A (zh) * | 2017-11-07 | 2018-03-30 | 中国石油大学(华东) | 一种确定天然气水合物bsr界面与海底相交位置的研究方法 |
Non-Patent Citations (2)
Title |
---|
徐浩等: "海域天然气水合物稳定带厚度与其影响因素的半定量研究", 《现代地质》 * |
毕海波: "台西南海域天然气水合物含量估算及地球化学特征分析", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109598019A (zh) * | 2018-10-26 | 2019-04-09 | 中国石油天然气股份有限公司 | 超低渗透油气藏微纳米孔隙泡点压力数值模拟方法及装置 |
CN109598019B (zh) * | 2018-10-26 | 2021-09-28 | 中国石油天然气股份有限公司 | 超低渗透油气藏微纳米孔隙泡点压力数值模拟方法及装置 |
CN110456024A (zh) * | 2019-08-23 | 2019-11-15 | 中国石油大学(华东) | 一种分析天然气水合物稳定区边界处碳循环过程的方法及系统 |
CN110954960A (zh) * | 2019-11-20 | 2020-04-03 | 中国石油大学(华东) | 一种海底冷泉地震响应数值模拟方法及系统 |
CN117408186A (zh) * | 2023-12-14 | 2024-01-16 | 青岛海洋地质研究所 | 地震剖面中天然气水合物似海底反射层的判断方法 |
CN117408186B (zh) * | 2023-12-14 | 2024-03-12 | 青岛海洋地质研究所 | 地震剖面中天然气水合物似海底反射层的判断方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108647461B (zh) | 2022-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Thompson et al. | The Antarctic Slope Current in a changing climate | |
CN108647461A (zh) | 盐底劈构造带水合物稳定底界的二维数值模拟方法及系统 | |
Moeck | Catalog of geothermal play types based on geologic controls | |
Gan et al. | A modeling study of shelf circulation off northern California in the region of the Coastal Ocean Dynamics Experiment: Response to relaxation of upwelling winds | |
Li et al. | The difficulty for subducted oceanic crust to accumulate at the Earth's core‐mantle boundary | |
Waples | Modeling of Sedimentary Basins and Petroleum Systems: Chapter 18: Part IV. Identification and Characterization | |
Qiu et al. | Deformation of a warm eddy in the northern South China Sea | |
Pham et al. | Assessment of CO2 injection into the south Utsira-Skade aquifer, the North Sea, Norway | |
CN109885962A (zh) | 海域天然气水合物分解诱发海底滑坡的数值模拟预测方法 | |
Cedeño et al. | The impact of salt tectonics on the thermal evolution and the petroleum system of confined rift basins: insights from basin modeling of the Nordkapp Basin, Norwegian Barents Sea | |
Barr et al. | Pre-development fracture modelling in the Clair field, west of Shetland | |
Pppelreiter et al. | Structural control on sweet-spot distribution in a carbonate reservoir: Concepts and 3-D models (Cogollo Group, Lower Cretaceous, Venezuela) | |
Aldenius | Subsurface characterization of the Lund Sandstone–3D model of the sandstone reservoir and evaluation of the geoenergy storage potential, SW Skåne, south Sweden-merged | |
Feng et al. | A high‐resolution 3‐D circulation model in a complex archipelago on the coastal Scotian shelf | |
CN110456024A (zh) | 一种分析天然气水合物稳定区边界处碳循环过程的方法及系统 | |
Brink | The evolution of the North German Basin and the metamorphism of the lower crust | |
Surrette et al. | Quantifying heterogeneity in variably fractured sedimentary rock using a hydrostructural domain | |
Ciftci et al. | Outcrop-based three-dimensional modeling of the Tensleep sandstone at Alkali Creek, Bighorn Basin, Wyoming | |
CN106869915A (zh) | 一种水平井井间隔夹层预测方法及装置 | |
Li et al. | Optimizing injection well trajectory to maximize storage security and minimize geomechanical risk | |
Pollack et al. | Stochastic Structural Modeling of a Geothermal Field: Patua Geothermal Field Case Study | |
Lescoffit et al. | Quantifying the impact of fault modeling parameters on production forecasting for clastic reservoirs | |
Harding et al. | The siliciclastic Upper Cretaceous play of eastern Mississippi Canyon | |
DeGraff et al. | Fracture network modeling and dual-permeability simulation of carbonate reservoirs | |
Grammer et al. | Carbonate platforms: Exploration-and production-scale insight from modern analogs in the Bahamas |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |