CN110456024B - 一种分析天然气水合物稳定区边界处碳循环过程的方法及系统 - Google Patents
一种分析天然气水合物稳定区边界处碳循环过程的方法及系统 Download PDFInfo
- Publication number
- CN110456024B CN110456024B CN201910784394.7A CN201910784394A CN110456024B CN 110456024 B CN110456024 B CN 110456024B CN 201910784394 A CN201910784394 A CN 201910784394A CN 110456024 B CN110456024 B CN 110456024B
- Authority
- CN
- China
- Prior art keywords
- hydrate
- boundary
- carbon
- natural gas
- coral reef
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/22—Fuels, explosives
- G01N33/222—Solid fuels, e.g. coal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/22—Fuels, explosives
- G01N33/225—Gaseous fuels, e.g. natural gas
Abstract
本发明属于天然气水合物碳循环技术领域,公开了一种分析天然气水合物稳定区边界处碳循环过程的方法及系统,基于地震解释、数值模拟、水合物分解量计算、珊瑚礁固碳量计算等方法,定量评估水合物稳定区边界处的碳循环过程。本发明通过计算水合物分解释放的游离气量和冷水珊瑚礁的固碳量,来估算进入海水中的甲烷量,从而充分认识水合物稳定区边界处,发生在岩石圈‑生物圈‑水圈‑大气圈中的碳循环过程,从而定量、准确评价天然气水合物系统在全球碳循环系统中的作用,帮助人们正确认识水合物对环境、气候变化等的影响。
Description
技术领域
本发明属于天然气水合物碳循环技术领域,尤其涉及一种分析天然气水合物稳定区边界处碳循环过程的方法及系统。
背景技术
目前,最接近的现有技术:
天然气水合物存储有大量的甲烷等烃类气体,其估算地质储量超过了全球陆地已知天然气地质储量总和。海域天然气水合物位于沉积物浅层,受海底洋流、沉积作用、底辟作用等影响天然气水合物稳定底界(Base of Gas Hydrate Stability Zone:BHSZ)易发生垂向迁移,导致天然气水合物分解。一般认为海底陆坡的天然气水合物稳定区GHSZ是一个向海增厚的楔状体。天然气水合物稳定区GHSZ向陆一侧分布边界(Landward Limit of GasHydrate Stability Zone:LLGHSZ)在地震剖面上表现为BSR与海底的相交处,该位置对海底的温度变化事件最敏感,天然气水合物易分解,且深部的水合物分解而释放的游离气还可沿水合物底界向上移动至此处,因此LLGHSZ位置处可能会存在较大的甲烷释放通量,释放的甲烷等气体若进入水体、甚至大气中,就会发生岩石圈-水圈-大气圈之间的碳循环,会造成一定的环境和气候影响。
甲烷作为一种大气痕量气体,其温室效应是二氧化碳的20多倍,其大量排放能够影响全球气候。前人研究通常用甲烷和气候变化的对比曲线来作为与过去气候变化有联系的甲烷释放的证据,认为很多环境和气候事件都可能与水合物分解有关,如新元古带雪球事件(600-800Ma)、P-T界线附近古环境与古生物灾变事件(251Ma)、侏罗纪Toarcian期的海洋缺氧事件(180Ma),以及古新世末的LPTM事件(55Ma)。水合物分解后,一部分甲烷在浅层沉积物中可进行甲烷厌氧氧化反应(Anaerobic Methane Oxidation:AMO),另外的甲烷进入海水中。前者AMO通常发育于海底沉积物表面至其下150cm的狭窄范围内,该带内CH4与硫酸盐浓度急剧下降,在贫/缺氧环境下以微生物为媒介的CH4氧化作用,通常与硫酸盐的还原作用(Sulfate Reduction:SR)同时发生。研究表明,水合物分解释放的甲烷90%以上可在沉积物中被AMO消耗,阻止甲烷继续向上运移。另外,进入海水中的甲烷还能在海水中与溶解氧发生反应,生成二氧化碳。这一反应非常强烈且反应速率很快,有学者认为进入海水中的甲烷也会消耗95%以上。
综上,目前学术界对水合物与气候环境之间的关系存在不同的认识,均已定性研究为主。有人认为水合物分解会造成环境和气候影响,还有人认为水合物分解释放的大多数甲烷会在抵达大气圈前被消耗。而本次申请者认为,水合物的存在对全球气候和环境具有双重作用。一方面,沉积物浅层水合物稳定区的存在起到了缓冲作用,防止/抑制了从地层深部运移至浅层的游离气直接逸散进入海水或者大气中;另一方面,当水合物稳定条件发生变化时,水合物分解释放的游离气可能会进入海水或者大气中,甚至诱发海底滑坡等灾害,造成一定的环境和气候影响。因此,十分有必要对天然气水合物稳定区边界处的碳循环过程进行定量研究,认识天然气水合物稳定底界动态迁移过程中,该处发生的碳循环过程,并定量评估其对环境、气候等的影响。
综上所述,现有技术存在的问题是:
目前的水合物研究中较少涉及天然气水合物系统在全球碳循环系统中的作用,仅有少量研究定性、粗略地描述了在天然气水合物稳定底界垂向迁移过程中发生在水合物层与下伏游离气区(Free Gas Zone:FGZ)之间的水合物循环,属于岩石圈内的碳循环,目前尚未有研究具体针对水合物稳定区边界处,发生在岩石圈-水圈-生物圈之间的碳循环过程。
解决上述技术问题的难度:
发生在天然气水合物稳定区边界处碳循环过程为动态过程,水合物稳定底界随地质过程,如沉积、剥蚀、海平面升降和地下热流异常等而发生垂向迁移,因此较难恢复和重建过去的水合物稳定底界的垂向迁移过程。本发明通过数值模拟的方法,恢复当珊瑚礁底界为古海底时的古水合物稳定底界(古BSR)的位置,进而分析发生在该时期的碳循环过程,并计算珊瑚礁在期间的固碳量、水合物分解释放的碳量等。
解决上述技术问题的意义:
本发明提供了一种研究水合物稳定区边界处碳循环过程的方法,该研究内容属于前沿领域,刻画了地质过程导致水合物稳定底界垂向迁移时,发生在水合物稳定区边界处的碳循环过程。原水合物稳定区边界处冷水珊瑚礁的固碳量被认为是由水合物碳储藏循环至生物圈的碳量,属于海洋生物固碳过程。通过计算水合物分解释放的游离气量和冷水珊瑚礁的固碳量,来估算进入海水中的甲烷量,从而充分认识水合物稳定区边界处,发生在岩石圈-生物圈-水圈-大气圈中的碳循环过程,从而准确评价天然气水合物系统在全球碳循环系统中的作用,帮助正确认识水合物对环境、气候变化等的影响。
发明内容
针对现有技术存在的问题,本发明提供了一种分析天然气水合物稳定区边界处碳循环过程的方法,基于地震解释、数值模拟、水合物分解量计算、珊瑚礁固碳量计算等方法,定量评估水合物稳定区边界处的碳循环过程,从而使人们正确认识天然气水合物系统在全球碳循环系统中的作用,以及水合物对环境、气候变化等的影响。
本发明是这样实现的,一种分析天然气水合物稳定区边界处碳循环过程的方法,包括以下操作步骤:
步骤1:选择天然气水合物分布边界附近发育冷水珊瑚礁条带的地区作为研究区。
步骤2:对BSR进行地震解释。
步骤3:对珊瑚礁顶部层位,即海底进行地震解释,提取海底构造属性平面图来掌握珊瑚礁的平面分布特征。
步骤4:对珊瑚礁底部层位进行地震解释。
步骤5:珊瑚礁顶底所在地层的双程走时为TWT(~50ms),则该地层深度为D(~36.25m)。
步骤6:在VVA软件中,通过刻画珊瑚礁地质体,来计算研究区珊瑚礁的面积A(~160km2)。
步骤7:通过调研,确定珊瑚礁的固碳速率为G,则珊瑚礁固碳量为MCcoral(~3.84*1015g),MCcoral=G*A*T。
步骤8:选择若干条珊瑚礁和BSR均发育的地震剖面,用于模拟当珊瑚礁底界为古海底时的古BSR及古LLGHSZ的位置,读取珊瑚礁顶底层位的深度信息作为古今海底深度。
步骤9:根据天然气水合物气体组分来选择水合物相平衡公式。
步骤10:根据研究区海平面变化曲线,结合现今海底深度,计算当珊瑚礁底界为海底时的原海底深度Dold-sb,从而计算当时古BSR处的静水压力值。
步骤11:根据研究区海水温度变化曲线,结合现今海底温度值Tsb,可恢复当珊瑚礁底界为海底时的原海底温度值Told-sb,从而计算当时古BSR处的温度值Told-bsr:
Told-bsr=Told-sb+GG*Dbsr,式中,GG为当时该地区的地温梯度。
步骤12:在步骤9-11的基础上,通过调整地温梯度参数,模拟当珊瑚礁底界为古海底时的古BSR和古LLGHSZ的位置。
步骤13:当LLGHSZ位置与珊瑚礁不一致时,调整模拟参数条件,重新进行模拟。
步骤14:当LLGHSZ位置与珊瑚礁一致时,说明水合物稳定区垂向运移之前,原水合物分布边界处的甲烷通量为该处的珊瑚礁发育提供了能量。
步骤15:将古今BSR之间的区域作为水合物分解区,计算该区的体积V(~1575*108m3),并估算在水合物稳定底界垂向迁移过程中水合物分解释放的甲烷体积VCH4(1.033*1012m3)和甲烷质量MCH4(7.4*1014g)和碳量MC CH4(5.6*1014g):
步骤16:通过对比步骤7计算的珊瑚礁固碳量和步骤15计算的水合物稳定底界垂向迁移过程中水合物分解释放的碳量,发现珊瑚礁固碳量大于水合物分解释放的碳量。由地震剖面可知,水合物稳定底界下部存在若干气烟囱,因此认为珊瑚礁除以水合物分解的甲烷作为能量外,气烟囱还可将深部的游离气垂向运移至古BSR下部,再沿古BSR侧向运移至古LLGHSZ处,为珊瑚礁提供能量。
步骤17:由以上研究可得出以下认识:LLGHSZ周边的深水珊瑚礁可大量固碳,岩石圈中的碳循环至生物圈中,并避免其进入大气中引起气候变化。
步骤18:针对珊瑚礁的巨大的固碳作用,在现今LLGHSZ条带处,可建议利用人工种植深水珊瑚礁固碳。
进一步,所述步骤5中,沉积速率为SR(~94.5cm/ky),则珊瑚礁的存在时间约为T(~40ky),其中:
D=TWT*Vsw/2000。
式中:Vsw为海水速度,Vsw=1450m/s;T=D/SR。
进一步,所述步骤7珊瑚礁固碳酸钙速率为1000-9000g CaCO3/m2/y,本研究区取5000g CaCO3/m2/y,相当于固碳率600g C/m2/y。
进一步,所述步骤9甲烷水合物相平衡公式如下:
Log10Pe=aTe 2+bTe+c。
式中:Pe和Te为水合物稳定底界处下的压力和温度值,a、b、c为经验常数,分别为a=0.000309℃-2,b=0.040094℃-1,c=0.478626。
进一步,所述步骤10中:
Pold-bsr=ρsw*g*(Dold-sb+Dbsr)。
式中:Dbsr为古海底至古BSR处所对应的深度,ρsw为海水密度,ρsw=1028kg/m3,g是重力加速度,g=9.81m/s2。
进一步,所述步骤15中:
式中:为沉积物平均孔隙度(研究区取~40%),SGH为估算的水合物饱和度(研究区取~10%),E为天然气水合物的容积倍率,取164。MCGH=VCH4*1000/22.4*16.04,这是由于1立方米甲烷为1000/22.4=44.643摩尔,相当于44.643*16.04=716克。
本发明的另一目的在于提供一种实施所述分析天然气水合物稳定区边界处碳循环过程的方法的分析天然气水合物稳定区边界处碳循环过程的控制系统。
本发明的另一目的在于提供一种实现所述分析天然气水合物稳定区边界处碳循环过程的方法的信息数据处理终端。
本发明的另一目的在于提供一种计算机可读存储介质,包括计算机程序,当其在计算机上运行时,使得计算机执行所述的分析天然气水合物稳定区边界处碳循环过程的方法。
综上所述,本发明的优点及积极效果为:
本发明提供了一种分析天然气水合物稳定区边界处碳循环过程的方法,本发明的内容属于前沿领域,刻画了地质过程导致天然气水合物稳定底界垂向迁移时,发生在天然气水合物稳定区边界处的碳循环过程。原天然气水合物稳定区边界处冷水珊瑚礁的固碳量被认为是由天然气水合物碳储藏循环至生物圈的碳量,属于海洋生物固碳过程。通过计算天然气水合物分解释放的游离气量和冷水珊瑚礁的固碳量,来估算进入海水中的甲烷量,从而充分认识天然气水合物稳定区边界处,发生在岩石圈-生物圈-水圈-大气圈中的碳循环过程,从而准确评价天然气水合物系统在全球碳循环系统中的作用,帮助正确认识天然气水合物对环境、气候变化等的影响。以某研究区为例,计算的珊瑚礁固碳量为3.84*1015g,而天然气水合物稳定底界垂向迁移过程中水合物分解释放的碳量5.6*1014g,因此,珊瑚礁固碳量大于水合物分解释放的碳量。分析认为,珊瑚礁除以水合物分解的甲烷作为能量外,水合物稳定底界下部存在的若干气烟囱还可将深部的游离气垂向运移至古BSR下部,再沿古BSR侧向运移至古LLGHSZ处,为珊瑚礁提供能量。
附图说明
图1是本发明实施例提供的一种研究天然气水合物稳定区边界处碳循环过程方法的流程图。
图2是本发明实施例提供的研究区海底倾角属性平面图,显示珊瑚礁条带和现今LLGHSZ的位置。
图3是本发明实施例提供的地震剖面显示珊瑚礁、现今BSR、现今LLGHSZ,数值模拟古BSR和古LLGHSZ、天然气水合物稳定区下伏气烟囱和高振幅异常等地震特征。
图4是本发明实施例提供的海水温度变化曲线和海平面温度变化曲线。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
目前的水合物研究中较少涉及天然气水合物系统在全球碳循环系统中的作用,仅有少量研究定性、粗略地描述了在水合物稳定底界垂向迁移过程中发生在水合物层与下伏游离气区(Free Gas Zone:FGZ)之间的水合物循环,属于岩石圈内的碳循环,目前尚未有研究具体针对水合物稳定区边界处碳循环过程
针对现有技术存在的问题,本发明提供了一种分析天然气水合物稳定区边界处碳循环过程的方法,下面结合附图对本发明作详细的描述。
如图1所示,本发明实施例提供一种分析天然气水合物稳定区边界处碳循环过程的方法,包括以下操作步骤:
S101:选择天然气水合物分布边界附近发育冷水珊瑚礁条带的地区作为研究区。
S102:对BSR进行地震解释。
S103:对珊瑚礁顶部层位,即海底进行地震解释,可提取海底构造属性平面图来掌握珊瑚礁的平面分布特征(如图2所示)。
S104:对珊瑚礁底部层位进行地震解释。
S105:珊瑚礁顶底所在地层的双程走时为TWT(~50ms),则该地层深度为D(~36.25m)。假设沉积速率为SR(~94.5cm/ky),则珊瑚礁的存在时间约为T(~40ky),其中:
D=TWT*Vsw/2000,式中,Vsw为海水速度,Vsw=1450m/s。T=D/SR
S106:在VVA软件中,通过刻画珊瑚礁地质体,来计算研究区珊瑚礁的面积A(~160km2)。
S107:通过调研,确定珊瑚礁的固碳速率为G(珊瑚礁固碳酸钙速率为1000-9000gCaCO3/m2/y,本研究区取5000g CaCO3/m2/y,相当于固碳率600g C/m2/y),则珊瑚礁固碳量为MCcoral(~3.84*1015g),MCcoral=G*A*T。
S108:选择若干条珊瑚礁和BSR均发育的地震剖面(如图3所示),用于模拟当珊瑚礁底界为古海底时的古BSR及古LLGHSZ的位置,读取珊瑚礁顶底层位的深度信息作为古今海底深度。
S109:根据天然气水合物气体组分来选择水合物相平衡公式,如甲烷水合物相平衡公式如下:
Log10Pe=aTe 2+bTe+c,式中Pe和Te为水合物稳定底界处下的压力和温度值,a、b、c为经验常数,分别为a=0.000309℃-2,b=0.040094℃-1,c=0.478626。
S110:根据研究区海平面变化曲线,结合现今海底深度,计算当珊瑚礁底界为海底时的原海底深度Dold-sb,从而计算当时古BSR处的静水压力值:
Pold-bsr=ρsw*g*(Dold-sb+Dbsr),式中,Dbsr为古海底至古BSR处所对应的深度,ρsw为海水密度,ρsw=1028kg/m3,g是重力加速度,g=9.81m/s2。
S111:根据研究区海水温度变化曲线(如图4所示),结合现今海底温度值Tsb,可恢复当珊瑚礁底界为海底时的原海底温度值Told-sb,从而计算当时古BSR处的温度值Told-bsr:
Told-bsr=Told-sb+GG*Dbsr,式中,GG为当时该地区的地温梯度。
S112:在S109-S111的基础上,通过调整地温梯度参数,模拟当珊瑚礁底界为古海底时的古BSR和古LLGHSZ的位置。
S113:当LLGHSZ位置与珊瑚礁不一致时,调整模拟参数条件,重新进行模拟。
S114:当LLGHSZ位置与珊瑚礁一致时,说明水合物稳定区垂向运移之前,原水合物分布边界处的甲烷通量为该处的珊瑚礁发育提供了能量。
S115:将古今BSR之间的区域作为水合物分解区,计算该区的体积V(~1575*108m3),并估算在水合物稳定底界垂向迁移过程中水合物分解释放的甲烷体积VCH4(1.033*1012m3)和甲烷质量MCH4(7.4*1014g)和碳量MC CH4(5.6*1014g):
式中,为沉积物平均孔隙度(研究区取~40%),SGH为估算的水合物饱和度(研究区取~10%),E为天然气水合物的容积倍率,取164。MCGH=VCH4*1000/22.4*16.04,这是由于1立方米甲烷为1000/22.4=44.643摩尔,相当于44.643*16.04=716克。
S116:通过对比S107计算的珊瑚礁固碳量和S115计算的水合物稳定底界垂向迁移过程中水合物分解释放的碳量,发现珊瑚礁固碳量大于水合物分解释放的碳量。由地震剖面可知,水合物稳定底界下部存在若干气烟囱,因此认为珊瑚礁除以水合物分解的甲烷作为能量外,气烟囱还可将深部的游离气垂向运移至古BSR下部,再沿古BSR侧向运移至古LLGHSZ处,为珊瑚礁提供能量。
S117:由以上研究可得出以下认识:LLGHSZ周边的深水珊瑚礁可大量固碳,岩石圈中的碳循环至生物圈中,并避免其进入大气中引起气候变化。
S118:针对珊瑚礁的巨大的固碳作用,可否考虑在现今LLGHSZ条带处,人工种植珊瑚礁,用于固碳。若可行,将会对全球碳循环造成重要影响。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种分析天然气水合物稳定区边界处碳循环过程的方法,其特征在于,所述分析天然气水合物稳定区边界处碳循环过程的方法基于地震解释、数值模拟、水合物分解量计算、珊瑚礁固碳量计算的方法,定量评估水合物稳定区边界处的碳循环过程;
最终获取天然气水合物系统在全球碳循环系统中的作用信息,以及水合物对环境、气候变化影响因素信息;
所述分析天然气水合物稳定区边界处碳循环过程的方法具体包括以下步骤:
步骤1:选择天然气水合物分布边界附近发育冷水珊瑚礁条带的地区作为研究区;
步骤2:对BSR进行地震解释;沉积速率为SR 94.5 cm/ky,珊瑚礁的存在时间为T40 ky,其中:
D =TWT*Vsw /2000;
式中:Vsw为海水速度,Vsw=1450 m/s;T=D/SR;
步骤3:对珊瑚礁顶部层位,即海底进行地震解释,提取海底构造属性平面图来掌握珊瑚礁的平面分布特征;
步骤4:对珊瑚礁底部层位进行地震解释;
步骤5:珊瑚礁顶底所在地层的双程走时为TWT,该地层深度为D;
步骤6:在VVA软件中,通过刻画珊瑚礁地质体,计算研究区珊瑚礁的面积A;
步骤7:通过调研,确定珊瑚礁的固碳速率为G,则珊瑚礁固碳量为MCcoral,MCcoral=G*A*T;T代表珊瑚礁的存在时间,单位是ky;
步骤8:选择若干条珊瑚礁和BSR均发育的地震剖面,用于模拟当珊瑚礁底界为古海底时的古BSR及古LLGHSZ的位置,读取珊瑚礁顶底层位的深度信息作为古今海底深度;
步骤9:根据天然气水合物气体组分来选择水合物相平衡公式;
步骤10:根据研究区海平面变化曲线,结合现今海底深度,计算当珊瑚礁底界为海底时的原海底深度Dold-sb,从而计算当时古BSR处的静水压力值;
步骤11:根据研究区海水温度变化曲线,结合现今海底温度值Tsb,恢复当珊瑚礁底界为海底时的原海底温度值Told-sb,计算当时古BSR处的温度值Told-bsr:
Told-bsr = Told-sb + GG *Dbsr;
式中,GG为当时该地区的地温梯度;Dbsr为古海底至古BSR处所对应的深度;
步骤12:通过步骤9-步骤11,通过调整地温梯度参数,模拟当珊瑚礁底界为古海底时的古BSR和古LLGHSZ的位置;
步骤13:LLGHSZ位置与珊瑚礁不一致时,调整模拟参数条件,重新进行模拟;
步骤14:当LLGHSZ位置与珊瑚礁一致时,说明水合物稳定区垂向运移之前,原水合物分布边界处的甲烷通量为该处的珊瑚礁发育提供了能量;
步骤15:将古今BSR之间的区域作为水合物分解区,计算该区的体积V,并估算在水合物稳定底界垂向迁移过程中水合物分解释放的甲烷体积VCH4和甲烷质量MCH4和碳量MC CH4;
步骤16:通过对比步骤7计算的珊瑚礁固碳量和步骤15计算的水合物稳定底界垂向迁移过程中水合物分解释放的碳量;
步骤17:分析LLGHSZ周边的深水珊瑚礁的固碳作用;
步骤18:在LLGHSZ条带处,利用人工种植深水珊瑚礁进行固碳。
2.如权利要求1中所述的分析天然气水合物稳定区边界处碳循环过程的方法,其特征在于,所述步骤7珊瑚礁的固碳速率为1000-9000 g CaCO3/m2/y。
3.如权利要求1中所述的分析天然气水合物稳定区边界处碳循环过程的方法,其特征在于,所述步骤9水合物相平衡公式如下:
Log10 Pe = aTe 2 + bTe + c;
式中:Pe和Te为水合物稳定底界处下的压力和温度值,a、b、c为经验常数,分别为a =0.000309 °C-2, b = 0.040094 °C-1,c = 0.478626。
4.如权利要求1中所述的分析天然气水合物稳定区边界处碳循环过程的方法,其特征在于,所述步骤10中 :
Pold-bsr= ρsw *g *(Dold-sb+Dbsr);
式中:Dbsr为古海底至古BSR处所对应的深度,ρsw为海水密度,ρsw = 1028 kg/m3,g是重力加速度,g=9.81 m/s2。
5.如权利要求1中所述的分析天然气水合物稳定区边界处碳循环过程的方法,其特征在于,所述步骤15中:
VCH4 = V*ø*SGH*E;
式中:ø为沉积物平均孔隙度,SGH为估算的水合物饱和度,E为天然气水合物的容积倍率,取164;MCCH4= VCH4 *1000/22.4*16.04,这是由于1立方米甲烷为1000/22.4=44.643摩尔,相当于44.643*16.04=716克。
6.一种实施权利要求1~5任意一项所述分析天然气水合物稳定区边界处碳循环过程的方法的分析天然气水合物稳定区边界处碳循环过程的控制系统。
7.一种实现权利要求1~5任意一项所述分析天然气水合物稳定区边界处碳循环过程的方法的信息数据处理终端。
8.一种计算机可读存储介质,包括计算机程序,当其在计算机上运行时,使得计算机执行如权利要求1-5任意一项所述的分析天然气水合物稳定区边界处碳循环过程的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910784394.7A CN110456024B (zh) | 2019-08-23 | 2019-08-23 | 一种分析天然气水合物稳定区边界处碳循环过程的方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910784394.7A CN110456024B (zh) | 2019-08-23 | 2019-08-23 | 一种分析天然气水合物稳定区边界处碳循环过程的方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110456024A CN110456024A (zh) | 2019-11-15 |
CN110456024B true CN110456024B (zh) | 2022-02-15 |
Family
ID=68488821
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910784394.7A Active CN110456024B (zh) | 2019-08-23 | 2019-08-23 | 一种分析天然气水合物稳定区边界处碳循环过程的方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110456024B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113156080B (zh) * | 2021-04-08 | 2022-03-25 | 青岛海洋地质研究所 | 一种模拟底辟活动对水合物成藏影响规律的装置及方法 |
CN113142102B (zh) * | 2021-04-30 | 2022-10-25 | 中国石油大学(华东) | 建立深水珊瑚海底养殖农场的选址及种植方法及装置 |
CN117211739B (zh) * | 2023-09-22 | 2024-04-09 | 青岛海洋地质研究所 | 一种海洋天然气水合物原位种植装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007034663A1 (ja) * | 2005-09-21 | 2007-03-29 | Idemitsu Kosan Co., Ltd. | 試料中の循環炭素含有量の測定方法及び試料中の循環炭素物質の含有率測定法 |
CN106459768A (zh) * | 2014-05-19 | 2017-02-22 | 杰富意钢铁株式会社 | 焦炭的制造方法和焦炭以及配煤的均质性的评价方法 |
CN106613735A (zh) * | 2016-12-23 | 2017-05-10 | 国家海洋局第三海洋研究所 | 一种提高红树林植被固碳能力的混交种植方法 |
AU2018201734A1 (en) * | 2015-08-12 | 2018-05-10 | Commonwealth Scientific And Industrial Research Organisation | Methanogenesis |
CN108370881A (zh) * | 2016-11-23 | 2018-08-07 | 雷学军 | 草本植物应对全球气候变化的方法 |
CN108826831A (zh) * | 2018-08-24 | 2018-11-16 | 杭州中泰深冷技术股份有限公司 | 氮循环制冷的深冷分离一氧化碳气体的装置和工艺 |
Family Cites Families (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5322035A (en) * | 1991-07-23 | 1994-06-21 | Space Biospheres Ventures | Hydrological system for a closed ecological system |
US6578405B2 (en) * | 2001-09-25 | 2003-06-17 | Schlumberger Technology Corporation | Gas seep detection |
US9519072B2 (en) * | 2006-05-11 | 2016-12-13 | Schlumberger Technology Corporation | Method and apparatus for locating gas hydrate |
US8882552B2 (en) * | 2007-12-29 | 2014-11-11 | Kal Karel Lambert | Biophysical geoengineering compositions and methods |
GB201108224D0 (en) * | 2011-05-17 | 2011-06-29 | Univ Heriot Watt | Early warning system for hydrate or clathrate materials |
KR101358013B1 (ko) * | 2012-03-13 | 2014-02-05 | 한국지질자원연구원 | 시추코어 내의 가스를 선상에서 추출, 분석, 저장하기 위한 장치 |
WO2014082018A1 (en) * | 2012-11-23 | 2014-05-30 | Fugro Geoconsulting, Inc. | Method and system for identification of gas hydrates and fee gas in geologic beds |
CN103439356B (zh) * | 2013-06-05 | 2016-08-10 | 中国石油大学(华东) | 用于超临界二氧化碳定压比热测量的实验装置与方法 |
US9645129B2 (en) * | 2013-06-27 | 2017-05-09 | Board Of Regents, The University Of Texas System | Process-based approach for the detection of deep gas invading the surface |
CN105989564B (zh) * | 2015-02-14 | 2022-02-11 | 雷学军 | 大气碳资源及co2当量物质综合开发利用的方法 |
CN105353101A (zh) * | 2015-09-24 | 2016-02-24 | 中国石油天然气股份有限公司 | 一种天然气成藏过程定量恢复的方法 |
CN105548438B (zh) * | 2016-01-26 | 2017-03-22 | 青岛海洋地质研究所 | 天然气水合物释放气体的连续高压氧化实验装置和方法 |
US20180230778A1 (en) * | 2016-06-17 | 2018-08-16 | Prostim Labs, Llc | System and Method for Permanent Storage of Carbon Dioxide in Shale Reservoirs |
CN107251724B (zh) * | 2016-09-23 | 2020-02-07 | 雷学军 | 发展新气候经济实现全球净零碳排放的方法 |
CN106629903B (zh) * | 2016-12-07 | 2019-08-09 | 大连理工大学 | 一种基于lng冷能的烟气水合物法海水淡化系统 |
CN106680872B (zh) * | 2017-02-06 | 2018-11-06 | 中国石油大学(华东) | 一种高效的海洋天然气水合物地震研究方法 |
CN110012785A (zh) * | 2017-07-04 | 2019-07-16 | 雷学军 | 动碳与静碳相互转化的方法 |
CN107462689A (zh) * | 2017-08-24 | 2017-12-12 | 新疆国利衡清洁能源科技有限公司 | 一种煤炭地下气化煤气逸散边界标定系统及其标定方法 |
CN107678055B (zh) * | 2017-08-29 | 2018-08-28 | 广州海洋地质调查局 | 一种海域天然气水合物海底甲烷监测系统及方法 |
CN107632134A (zh) * | 2017-10-23 | 2018-01-26 | 大庆东油睿佳石油科技有限公司 | 一种天然气水合物海水驱光刻玻璃模型微观实验装置 |
CN107861158A (zh) * | 2017-11-07 | 2018-03-30 | 中国石油大学(华东) | 一种确定天然气水合物bsr界面与海底相交位置的研究方法 |
CN107930537A (zh) * | 2017-11-15 | 2018-04-20 | 中国科学院广州能源研究所 | 一种模拟海底甲烷渗漏导致早期成岩作用的反应装置及方法 |
CN108647461B (zh) * | 2018-05-17 | 2022-04-19 | 中国石油大学(华东) | 盐底劈构造带水合物稳定底界的二维数值模拟方法及系统 |
CN109569473B (zh) * | 2018-11-21 | 2020-11-24 | 太原理工大学 | 一种液态金属催化烃生产氢气和炭黑的装置及方法 |
CN109885962B (zh) * | 2019-03-05 | 2023-09-26 | 中国石油大学(华东) | 海域天然气水合物分解诱发海底滑坡的数值模拟预测方法 |
-
2019
- 2019-08-23 CN CN201910784394.7A patent/CN110456024B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007034663A1 (ja) * | 2005-09-21 | 2007-03-29 | Idemitsu Kosan Co., Ltd. | 試料中の循環炭素含有量の測定方法及び試料中の循環炭素物質の含有率測定法 |
CN106459768A (zh) * | 2014-05-19 | 2017-02-22 | 杰富意钢铁株式会社 | 焦炭的制造方法和焦炭以及配煤的均质性的评价方法 |
AU2018201734A1 (en) * | 2015-08-12 | 2018-05-10 | Commonwealth Scientific And Industrial Research Organisation | Methanogenesis |
CN108370881A (zh) * | 2016-11-23 | 2018-08-07 | 雷学军 | 草本植物应对全球气候变化的方法 |
CN106613735A (zh) * | 2016-12-23 | 2017-05-10 | 国家海洋局第三海洋研究所 | 一种提高红树林植被固碳能力的混交种植方法 |
CN108826831A (zh) * | 2018-08-24 | 2018-11-16 | 杭州中泰深冷技术股份有限公司 | 氮循环制冷的深冷分离一氧化碳气体的装置和工艺 |
Non-Patent Citations (5)
Title |
---|
An Areal Assessment of Subseafloor Carbon Cycling in Cold Seeps and Hydrate-Bearing Areas in the Northern South China Sea;Yanping Zhang et al.;《Geofluids》;20190317;第1-14页 * |
Lee, M. W.Amplitude Variation of Bottom Simulating Reflection with Respect to Frequency: Transitional Base or Attenuation * |
Scientific investigations rept.《国家工程技术数字图书馆》.2007,第1-16页. * |
南沙珊瑚礁生态系的碳循环;宋金明 等;《海洋与湖沼》;20031231;第34卷(第6期);第586-592页 * |
珊瑚礁区碳循环研究进展;严宏强 等;《生态学报》;20091231;第29卷(第11期);第6207-6215页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110456024A (zh) | 2019-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110456024B (zh) | 一种分析天然气水合物稳定区边界处碳循环过程的方法及系统 | |
Wallmann et al. | Kinetics of organic matter degradation, microbial methane generation, and gas hydrate formation in anoxic marine sediments | |
Nöthen et al. | Reconstructing changes in seep activity by means of pore water and solid phase Sr/Ca and Mg/Ca ratios in pockmark sediments of the Northern Congo Fan | |
Fantle et al. | Sr isotopes and pore fluid chemistry in carbonate sediment of the Ontong Java Plateau: Calcite recrystallization rates and evidence for a rapid rise in seawater Mg over the last 10 million years | |
Luff et al. | Fluid flow, methane fluxes, carbonate precipitation and biogeochemical turnover in gas hydrate-bearing sediments at Hydrate Ridge, Cascadia Margin: numerical modeling and mass balances | |
Langer et al. | Frozen ponds: production and storage of methane during the Arctic winter in a lowland tundra landscape in northern Siberia, Lena River delta | |
CN110632666B (zh) | 一种碳酸盐岩溶蚀孔洞分布的预测方法及装置 | |
Luff et al. | Simulation of long-term feedbacks from authigenic carbonate crust formation at cold vent sites | |
Paull et al. | Authigenic carbon entombed in methane-soaked sediments from the northeastern transform margin of the Guaymas Basin, Gulf of California | |
Haeckel et al. | Pore water profiles and numerical modelling of biogeochemical processes in Peru Basin deep-sea sediments | |
Malone et al. | Variable methane fluxes in shallow marine systems over geologic time: The composition and origin of pore waters and authigenic carbonates on the New Jersey shelf | |
Belmaker et al. | 10Be dating of Neogene halite | |
Ge et al. | Sr isotopic compositions of cold seep carbonates from the South China Sea and the Panoche Hills (California, USA) and their significance in palaeooceanography | |
US20220187227A1 (en) | Method of estimating a mineral content of a geological structure | |
Paull et al. | Seawater circulation through the flank of the Florida platform: evidence and implications | |
Wang et al. | Expansion of the Southern California oxygen minimum zone during the early-to mid-Holocene due to reduced ventilation of the Northeast Pacific | |
Hay et al. | Evolution of sediment fluxes and ocean salinity | |
Staudigel et al. | An abrupt Middle-Miocene increase in fluid flow into the Leeward Margin Great Bahama Bank, constraints from δ44Ca and Δ47 values | |
Chen et al. | Response of Shallow Gas‐Charged Holocene Deposits in the Yangtze Delta to Meter‐Scale Erosion Induced by Diminished Sediment Supply: Increasing Greenhouse Gas Emissions | |
Mienert et al. | Finding and using the world’s gas hydrates | |
Thomas et al. | Influence of Pleistocene sea-level lows on offshore fresh groundwater reservoirs-a numerical case study on the New Jersey shelf | |
Thomas et al. | Influence of the Last Glacial Maximum on New Jersey shelf offshore fresh groundwater reservoirs-investigating the role of geological heterogeneity | |
Hoffman | Spatial variability of erosion patterns along the eastern margin of the Rio Grande Rift | |
Sim et al. | Are Seismic Velocity Time-Lapse Changes Due to Fluid Substitution or Matrix Dissolution? A CO2 Sequestration Study at Pohokura Field, New Zealand. | |
Pengfei et al. | In situ quantification of oil-gas bubble seep flux from cold seeps at the seabed―a case study of GC600 in the Gulf of Mexico |
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 |