CN108804860B - 气烟囱对天然气水合物生成和分布的影响分析系统及方法 - Google Patents
气烟囱对天然气水合物生成和分布的影响分析系统及方法 Download PDFInfo
- Publication number
- CN108804860B CN108804860B CN201811088537.2A CN201811088537A CN108804860B CN 108804860 B CN108804860 B CN 108804860B CN 201811088537 A CN201811088537 A CN 201811088537A CN 108804860 B CN108804860 B CN 108804860B
- Authority
- CN
- China
- Prior art keywords
- gas
- hydrate
- chimney
- resistivity
- distribution
- 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
- 239000007789 gas Substances 0.000 title claims abstract description 152
- NMJORVOYSJLJGU-UHFFFAOYSA-N methane clathrate Chemical compound C.C.C.C.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O NMJORVOYSJLJGU-UHFFFAOYSA-N 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 29
- 239000013049 sediment Substances 0.000 claims abstract description 67
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 57
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 26
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 13
- 238000013508 migration Methods 0.000 claims abstract description 10
- 230000005012 migration Effects 0.000 claims abstract description 10
- 239000012530 fluid Substances 0.000 claims abstract description 9
- 150000004677 hydrates Chemical class 0.000 claims abstract description 9
- 230000002093 peripheral effect Effects 0.000 claims abstract description 8
- 238000012544 monitoring process Methods 0.000 claims abstract description 7
- 238000012360 testing method Methods 0.000 claims abstract description 7
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Natural products C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 46
- 238000006243 chemical reaction Methods 0.000 claims description 32
- 239000011148 porous material Substances 0.000 claims description 29
- 238000002347 injection Methods 0.000 claims description 23
- 239000007924 injection Substances 0.000 claims description 23
- 238000002474 experimental method Methods 0.000 claims description 17
- 238000005755 formation reaction Methods 0.000 claims description 13
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 claims description 12
- 239000002245 particle Substances 0.000 claims description 10
- 239000003345 natural gas Substances 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 7
- -1 natural gas hydrates Chemical class 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 239000000243 solution Substances 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 6
- 239000011435 rock Substances 0.000 claims description 6
- 229920006395 saturated elastomer Polymers 0.000 claims description 6
- 239000011780 sodium chloride Substances 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 239000012153 distilled water Substances 0.000 claims description 3
- 238000013401 experimental design Methods 0.000 claims description 3
- 230000006641 stabilisation Effects 0.000 claims description 3
- 238000011105 stabilization Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims 1
- 238000011161 development Methods 0.000 abstract description 5
- 238000011160 research Methods 0.000 description 6
- 238000005192 partition Methods 0.000 description 4
- 238000005086 pumping Methods 0.000 description 4
- 230000000035 biogenic effect Effects 0.000 description 3
- 230000008021 deposition Effects 0.000 description 3
- 229910001069 Ti alloy Inorganic materials 0.000 description 2
- 239000000956 alloy Substances 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000011065 in-situ storage Methods 0.000 description 2
- 230000000476 thermogenic effect Effects 0.000 description 2
- 238000004140 cleaning Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 235000021185 dessert Nutrition 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 239000008398 formation water Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Abstract
本发明属于天然气水合物勘探开发技术领域，公开了一种气烟囱对天然气水合物生成和分布的影响的系统及方法，模拟气烟囱作为流体运移通道对天然气水合物稳定区供给气源，通过电阻率监测结果计算水合物饱和度；设置隙水盐度、沉积物粒度、实验温度和压力不同的地质条件参数，再通过设置不同的进气速度；模拟不同地质条件下，不同气源供给强度下气烟囱内部及周边沉积物中的水合物分布。本发明针对不同位置的天然气水合物合成样品，包括气烟囱内部和沉积地层距离气烟囱不同位置处，分别测试其水合物饱和度参数，掌握当气烟囱为天然气水合物体系提供气源时所生成水合物的空间分布特征。
Description
技术领域
本发明属于天然气水合物勘探开发技术领域，尤其涉及一种气烟囱对天然气水合物生成和分布的影响分析系统及方法。
背景技术
目前，业内常用的现有技术是这样的：
充足的气源是形成高饱和度天然气水合物藏的基础。对气源条件的深入分析不仅有助于了解天然气水合物体系的成藏机理，还会影响天然气水合物的资源调查、评价，甚至开发利用。水合物的气源包括原地生成的生物成因气和深部垂向运移而来的热成因气，或者是二者的混合气。一般认为，原地生成的生物成因气气量较少，不足以形成丰富的天然气水合物藏，因此沿运移通道运移来的深部生物成因或热成因气，通常被认为是形成天然气水合物的主要气源。
与天然气水合物相关的流体运移通道包括气烟囱、断层、底辟构造等，前人研究发现天然气水合物体系的发育与流体运移通道关系十分密切。如在墨西哥湾的Green Canyon地区，GC955-Q和GC955-H两口井相距1km。结合测井和地震数据发现，在浅部地层，Q井无水合物发育，而H井在海底下部250-350m的深部发育有裂缝型水合物，这可能是由H井处发育有气烟囱，而Q井处未发育垂向流体运移通道所造成。
但目前关于气烟囱如何影响天然气水合物生成和分布的研究并不深入。有研究者认为，在天然气水合物稳定区GHSZ，气烟囱中会形成高饱和度的水合物，堵塞管道，代表了潜在的天然气水合物甜点区；但也有研究者认为在GHSZ，当游离气在气烟囱垂向上移时会在壁上形成天然气水合物，饱和度向气烟囱外侧降低，形成非渗透性的气烟囱壁。
而在气烟囱中间，由于水合物形成过程伴随的排盐效应导致的地层水盐度增高和孔隙水的缺失等因素，天然气水合物的生成受到抑制，形成中空的通道，促进流体沿气烟囱向上移动。另外，在气烟囱周边的沉积物中，水合物饱和度是否与和气烟囱的距离存在相关性等都是目前尚未进行探索的问题。
综上所述，目前水合物研究中存在的问题是：
(1)气烟囱内部是被高饱和度水合物堵塞，还是为中空通道，不同学者间认识不同，尚未进行实验或取样验证；
(2)在气烟囱周边的沉积物中，水合物饱和度是否与和气烟囱的距离存在相关性等也尚未进行探索。
(3)针对不同位置的天然气水合物合成样品，没有分别测试其水合物饱和度参数，不能掌握当气烟囱为天然气水合物体系提供气源时所生成水合物的空间分布特征。
解决上述技术问题的难度和意义：
为了直观掌握和了解气烟囱内部和周边沉积物中的水合物分布情况，且模拟实际沉积地层粒度粗细相间的特征，需要将反应釜内部用带孔隔板分隔为若干层，并且在不同位置设计若干电阻率传感器来监测水合物合成情况。另外，利用本装置可以设置不同的地质条件参数，如孔隙水盐度、沉积物粒度、实验温度和压力条件等，再通过设置不同的进气速度，可模拟不同地质条件下，不同气源供给强度下气烟囱及周边沉积物中的水合物分布研究。
发明内容
针对现有技术存在的问题，本发明提供了一种气烟囱对天然气水合物生成和分布的影响分析系统及方法。
本发明是这样实现的，一种气烟囱对天然气水合物生成和分布的影响分析方法，所述烟囱对天然气水合物生成和分布的影响方法包括：
模拟气烟囱作为流体运移通道对天然气水合物稳定区供给气源，通过电阻率监测结果计算水合物饱和度；
设置隙水盐度、沉积物粒度、实验温度和压力不同的地质条件参数，再通过设置不同的进气速度；模拟不同地质条件下，分析不同气源供给强度下气烟囱内部及周边沉积物中的水合物分布。
进一步，所述烟囱对天然气水合物生成和分布的影响方法具体包括：
步骤一，通过测量海底沉积物样品，确定实验所用沉积物粒径、孔隙度、孔隙水盐度、沉积物体积含水量、甲烷气体流速，以及温度、围压值和轴压值实验参数；
步骤二:准备实验材料：纯度达到99.99％的甲烷气体、盐度为35psu的氯化钠溶液，粒度分别为0.125-0.18mm、0.045-0.1mm和0.25-0.5mm的沉积物；
步骤三:采用蒸馏水清洗高压反应釜，清洗3-5遍；
步骤四:注入甲烷气体并保持一定压力静置，实验装置不漏气，再关闭注气口阀门；
步骤五:恒温箱温度设置为1℃，整个实验过程中保持不变；
步骤六:将0.125-0.18mm和0.045-0.1mm不同粒度沉积物分层填入样品腔中，用带孔隔板将不同粒度沉积物层分隔开，气烟囱内放置0.25-0.5mm的沉积物，并按照实验设计在不同位置设置电阻率和温度传感器；
步骤七:加入盐度为35psu的NaCl溶液，沉积物孔隙空间被100％孔隙水饱和，上下左右摇晃在沉积物中均匀分布；
步骤八:通过轴压和围压控制系统分别控制样品腔的轴压和围压；
步骤九:静置1小时，并测试电阻率参数在期间的变化；
步骤十:设置注气气体流速为0-2000mL/min，精度为0.25F.S，气体增压上限为25MPa；向气烟囱注入甲烷气，开始水合物合成，并通过反应釜顶部排气系统控制反应釜内部压力值；
步骤十一：持续稳定注气，直至各电阻率传感器所测电阻率值保持不变或者气烟囱堵塞无法注气，实验结束；
步骤十二：收集并处理实验数据，利用Archie公式计算水合物饱和度。
进一步，利用Archie公式计算水合物饱和度，具体包括：
岩石孔隙100％被孔隙水饱和的沉积物中，电阻率R0与孔隙水电阻率Rw成正比，比例系数为地层因子F：
式中，a为初始沉积物岩性系数；是储层孔隙度；m为胶结系数；
含水合物沉积物中，电阻率Rt与100％含水的沉积物电阻率R0之比为电阻率增大指数I：
式中，b为含水合物沉积物岩性系数；Sw是储层含水饱和度；n是饱和度指数；
通过不同途径计算水合物饱和度Sh：
式中，a、b、m、n通过岩电实验确定，a和m与岩性和孔隙结构相关，b和n与岩石的润湿性和含油气性相关；b取值1，n取值1.94。
进一步，步骤六中，电阻率和温度传感器分别放置在气烟囱的上、中、下位置，以及粗细沉积地层中距离气烟囱的不同位置处，用来监测水合物的分布情况。
本发明的另一目的在于提供一种计算机程序，所述计算机程序运行实现所述的气烟囱对天然气水合物生成和分布的影响分析方法。
本发明的另一目的在于提供一种终端，所述终端至少搭载实现所述气烟囱对天然气水合物生成和分布的影响分析方法的控制器。
本发明的另一目的在于提供一种计算机可读存储介质，包括指令，当其在计算机上运行时，使得计算机执行所述的气烟囱对天然气水合物生成和分布的影响分析方法。
本发明的另一目的在于提供一种实现所述气烟囱对天然气水合物生成和分布的影响分析方法的气烟囱对天然气水合物生成和分布的影响分析系统，包括：
高压反应釜；
高压反应釜放置于恒温箱内，高压反应釜中下部安装有气烟囱；
高压反应釜焊接的多块带孔隔板将高压反应釜内部空间分隔为若干层；
高压反应釜和气烟囱内部安装有若干电阻率传感器和温度传感器；电阻率传感器和温度传感器连接至恒温箱外部的数据采集器上；
数据采集器连接电脑，气烟囱下部通过注气口连接阀门，阀门连接增压泵，增压泵连接甲烷气瓶；
高压反应釜顶部通过导管连接阀门，阀门连接压力计，压力计连接排气瓶。
进一步，增压阀用于在实验过程中调节甲烷气体的供给速度，气体流速为0-2000mL/min，精度为0.25F.S，气体增压上限为25MPa；
恒温箱的温度控制系统为恒温空气浴，温度为-80℃～室温，恒温波动度±0.1℃；
高压反应釜为圆柱体外形，采用TC4钛合金材料，内径为40cm，高为40cm，容积为50L，最高工作压力为30MPa；
气烟囱为壁上布满小孔的塑料管，内部填充粗粒沉积物代表高渗透性气体运移通道；
带孔隔板放入0.1-0.25mm和0.05-0.1mm不同粒度沉积物代表粗细沉积物互层；
带孔隔板上的小孔孔径小于沉积物最小粒径；电脑通过相关软件实时监控和保存测试数据；
进一步，所述的气烟囱对天然气水合物生成和分布的影响分析系统进一步包括：
围压控制系统，通过手摇注水泵向围压腔内注水或者从围压腔内抽水，控制高压反应釜的围压值，当围压压力传感器测得的压力值大于围压设定值时，通过手摇注水泵从围压腔抽水，当围压压力传感器测得的压力值小于围压设定值时，通过手摇注水泵向围压腔注水；
轴压控制系统，通过手摇注水泵向轴压腔内注水或者从轴压腔内抽水，控制高压反应釜的轴压值，当轴压压力传感器测得的压力值大于轴压设定值时，通过手摇注水泵从轴压腔抽水，当轴压压力传感器测得的压力值小于围压设定值时，通过手摇注水泵向轴压腔注水。
综上所述，本发明的优点及积极效果为：
本发明提供一种研究气烟囱对天然气水合物生成和分布的影响的装置及方法，确定当气烟囱对天然气水合物体系进行气源供给时，气烟囱内部及周边沉积物中的天然气水合物分布特征。针对不同位置的天然气水合物合成样品，包括气烟囱内部和沉积地层距离气烟囱不同位置处，分别测试其水合物饱和度参数，掌握当气烟囱为天然气水合物体系提供气源时所生成水合物的空间分布特征。
本发明的优点还有：
目前，前人对气烟囱和周边沉积物的水合物分布特征研究尚不深入，不同学者间存在争议且缺少直接证据，而本发明首次模拟气烟囱作为流体运移通道对天然气水合物稳定区供给气源，通过电阻率监测结果来计算水合物饱和度，从而可以直观认识和了解气烟囱内部和周边沉积物中的水合物分布情况；
本发明可以设置不同的地质条件参数，如孔隙水盐度、沉积物粒度、实验温度和压力条件等，再通过设置不同的进气速度，可以模拟不同地质条件下，不同气源供给强度下气烟囱及周边沉积物中的水合物分布研究。
附图说明
图1是本发明实施例提供的气烟囱对天然气水合物生成和分布的影响分析系统示意图。
图2是本发明实施例提供的高压反应釜结构示意图。
图中：1、甲烷气瓶；2、增压泵；3、阀门；4、恒温箱；5、高压反应釜；6、气烟囱；7、带孔隔板；8、电阻率传感器；9、温度传感器；10、数据采集器；11、电脑；12、阀门；13、压力计；14、排气瓶；15、注气口；16、围压腔；17、手摇注水泵；18、压力计；19、轴压腔；20、压力计；21、手摇注水泵。
图3是本发明实施例提供的气烟囱对天然气水合物生成和分布的影响分析方法流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
针对不同位置的天然气水合物合成样品，没有分别测试其水合物饱和度参数，不能掌握当气烟囱为天然气水合物体系提供气源时所生成水合物的空间分布特征。
如图1-图2所示，本发明实施例提供的气烟囱对天然气水合物生成和分布的影响分析系统，由恒温箱、高压反应釜和数据采集系统构成。高压反应釜5放置于恒温箱4内，高压反应釜5中下部设计有气烟囱6，高压反应釜5设计有带孔隔板7将内部空间分隔为若干层，高压反应釜5和气烟囱6内部设置若干电阻率传感器8和温度传感器9，电阻率传感器8和温度传感器9连接至恒温箱4外部的数据采集器10上，数据采集器10连接电脑11，气烟囱6下部通过注气口15连接阀门3，阀门3连接增压泵2，增压泵2连接甲烷气瓶1，高压反应釜5顶部通过导管连接阀门12，阀门12连接压力计13，压力计13连接排气瓶14。
增压阀2可用来在实验过程中调节甲烷气体的供给速度，代表向上运移的不同强度的深部流体，气体流速范围为0-2000mL/min，精度为0.25F.S，气体增压上限为25MPa；
恒温箱4的温度控制系统为恒温空气浴，温度控制范围为-80℃至室温，最佳恒温波动度可达±0.1℃；
高压反应釜5为圆柱体外形，采用TC4钛合金材料制作，内径为40cm，高为40cm，容积约为50L，最高工作压力为30MPa；
气烟囱6为壁上布满小孔的塑料管，实验时内部填充粗粒沉积物代表高渗透性气体运移通道；
带孔隔板7将高压反应釜5的空间分隔为若干层，实验时放入不同粒度沉积物代表粗细沉积物互层，如0.1-0.25mm和0.05-0.1mm，也可选择其他粒径的沉积物；
带孔隔板7上的小孔大量均匀分布，且孔径小于沉积物最小粒径；
若干个电阻率传感器8和温度传感器9分别放置于气烟囱内部和各沉积地层距离气烟囱的不同位置处，用来监测水合物的生成和分布情况；
电脑11可通过相关软件实时监控和保存测试数据；
围压控制系统通过压力计18监测围压，利用手摇注水泵17对围压腔16进行注水或抽水，从而控制高压反应釜5的围压值，该方法对围压腔的体积改变较小；
轴压控制系统通过压力计20监测轴压，利用手摇注水泵21对轴压腔19进行注水或抽水，控制高压反应釜5的轴压值，该方法对轴压腔的体积改变较小。
如图3所示，本发明实施例提供的气烟囱对天然气水合物生成和分布的影响分析方法，包括以下步骤：
S101:通过实际测量海底沉积物样品和查阅文献，确定实验所用沉积物粒径、孔隙度、孔隙水盐度、沉积物体积含水量、甲烷气体流速，以及温度、围压值和轴压值等实验参数；
S102:准备实验材料：纯度达到99.99％的甲烷气体、盐度为35psu的氯化钠溶液，粒度分别为0.125-0.18mm、0.045-0.1mm和0.25-0.5mm的沉积物；
S103:采用蒸馏水清洗高压反应釜，清洗3-5遍；
S104:注入甲烷气体并保持一定压力静置，确保实验装置不漏气，再关闭注气口阀门；
S105:恒温箱温度设置为1℃，在整个实验过程中保持不变；
S106:将不同粒度沉积物(0.125-0.18mm和0.045-0.1mm)分层填入样品腔中，用带孔隔板将不同粒度沉积物层分隔开，气烟囱内放置0.25-0.5mm的沉积物，并按照实验设计在不同位置设置电阻率和温度传感器；
S107:加入盐度为35psu的NaCl溶液，沉积物孔隙空间被100％孔隙水饱和，轻轻上下左右摇晃使其开始时在沉积物中均匀分布；
S108:参考海底沉积物的实际地质条件，通过轴压和围压控制系统分别控制样品腔的轴压和围压；
S109:静置1小时，并测试电阻率参数在期间的变化；
S110:设置一定的注气速度，开始持续不断向气烟囱注入甲烷气，开始水合物合成，并通过反应釜顶部排气系统控制反应釜内部压力值；
S111：持续稳定注气，直至各电阻率传感器所测电阻率值保持不变或者气烟囱堵塞无法注气，实验结束；
S112：收集并处理实验数据，利用Archie公式计算水合物饱和度。
计算过程如下：
对岩石孔隙100％被孔隙水饱和的沉积物而言，其电阻率R0与孔隙水电阻率Rw成正比，其比例系数即为地层因子F：
式中，a为初始沉积物岩性系数；是储层孔隙度；m为胶结系数。
对含水合物沉积物而言，其电阻率Rt与100％含水的沉积物电阻率R0之比可定义为电阻率增大指数I：
式中，b为含水合物沉积物岩性系数；Sw是储层含水饱和度；n是饱和度指数。
由公式(1)和(2)，可通过不同途径计算水合物饱和度Sh：
式中，a、b、m、n通常可通过岩电实验确定，a和m与岩性和孔隙结构相关，b和n与岩石的润湿性和含油气性相关。
本发明拟采用公式(4)计算水合物饱和度，R0和Rt可通过实验测量得到，b和n可通过调研工区岩电实验等结果确定，本发明拟将b取值1，n取值1.94。另外，在具体工区应用时可将该计算结果与岩心分析结果进行对比，修正饱和度计算中所用到的阿尔奇参数。
步骤S102可根据实验目的准备不同盐度的孔隙水和不同粒径的沉积物，通过水合物合成实验来分析不同盐度下或不同粒度沉积物地层中的水合物分布情况；
步骤S105设置不同的实验温度参数，来分析不同温度下气烟囱周边的水合物合成和分布情况；
步骤S106中分层放入0.125-0.18mm和0.045-0.1mm沉积物用来代表粗细相间的沉积地层，气烟囱中放入0.25-0.5mm的粗粒沉积物用来指示气烟囱的高渗透性；
步骤S106中沉积物并非填满整个样品腔，而是在顶部预留一部分空间，用来排出气烟囱持续供给的甲烷气；
步骤S106中将电阻率和温度传感器分别放置在气烟囱的上、中、下位置，以及粗细沉积地层中距离气烟囱的不同位置处，用来监测水合物的分布情况；
步骤S106中由于设置了很多个电阻率和温度传感器，因此认为每个数据点处的沉积物孔隙度和水合物饱和度是均匀分布的，符合阿尔奇公式的适用范围；
步骤S109在水合物开始合成前静置沉积物，并监测电阻率变化，目的是分析孔隙水在垂向上是否为不均匀分布；
步骤S110中通过设置不同的注气速度，分析气源供给强度对水合物合成和分布的影响；
步骤S111中的实验结束条件不确定，是因为目前尚不确定气烟囱为水合物稳定区供给气源时水合物在气烟囱内部及周边的分布情况，因此需要考虑各种可能发生的情况。
本发明拟采用公式(4)计算水合物饱和度，R0和Rt可通过实验测量得到，b和n可通过调研工区岩电实验等结果确定，本发明拟将b取值1，n取值1.94。另外，在具体工区应用时可将该计算结果与岩心分析结果进行对比，修正饱和度计算中所用到的阿尔奇参数。
所述步骤S102、S105、S106、S107、S108、S110中可以设置不同的地质条件参数，如孔隙水盐度、沉积物粒度、实验温度和压力条件等，还可以设置不同的注甲烷气速度，因此该实验方法可以模拟不同地质条件下、气烟囱的不同气源供给强度下，水合物在气烟囱内部及周边沉积物中的分布特征。
步骤S112中利用Archie公式计算水合物饱和度，该计算可利用不同公式进行，公式的选择取决于实验中测试的电阻率参数，以及模拟目标区的岩电实验结果中可用的阿尔奇参数。
本发明实施例提供的气烟囱对天然气水合物生成和分布的影响分析系统，进一步包括：
围压控制系统，通过手摇注水泵向围压腔内注水或者从围压腔内抽水，控制高压反应釜的围压值，当围压压力传感器测得的压力值大于围压设定值时，通过手摇注水泵从围压腔抽水，当围压压力传感器测得的压力值小于围压设定值时，通过手摇注水泵向围压腔注水；
轴压控制系统，通过手摇注水泵向轴压腔内注水或者从轴压腔内抽水，控制高压反应釜的轴压值，当轴压压力传感器测得的压力值大于轴压设定值时，通过手摇注水泵从轴压腔抽水，当轴压压力传感器测得的压力值小于围压设定值时，通过手摇注水泵向轴压腔注水。
在上述实施例中，可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用全部或部分地以计算机程序产品的形式实现，所述计算机程序产品包括一个或多个计算机指令。在计算机上加载或执行所述计算机程序指令时，全部或部分地产生按照本发明实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中，或者从一个计算机可读存储介质向另一个计算机可读存储介质传输，例如，所述计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线(DSL)或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输)。所述计算机可读取存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质，(例如，软盘、硬盘、磁带)、光介质(例如，DVD)、或者半导体介质(例如固态硬盘SolidState Disk(SSD))等。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (4)
1.一种气烟囱对天然气水合物生成和分布的影响分析方法，其特征在于，所述气烟囱对天然气水合物生成和分布的影响分析方法包括：
模拟气烟囱作为流体运移通道对天然气水合物稳定区供给气源，通过电阻率监测结果计算水合物饱和度；
设置孔隙水盐度、沉积物粒度、实验温度和压力不同的地质条件参数，再通过设置不同的进气速度；模拟不同地质条件下，分析不同气源供给强度下气烟囱内部及周边沉积物中的水合物分布；
所述气烟囱对天然气水合物生成和分布的影响分析方法具体包括：
步骤一，通过测量海底沉积物样品，确定实验所用沉积物粒径、孔隙度、孔隙水盐度、沉积物体积含水量、甲烷气体流速，以及温度、围压值和轴压值实验参数；
步骤二:准备实验材料：纯度达到99.99％的甲烷气体、盐度为35psu的氯化钠溶液，粒度分别为0.125-0.18mm、0.045-0.1mm和0.25-0.5mm的沉积物；
步骤三:采用蒸馏水清洗高压反应釜，清洗3-5遍；
步骤四:注入甲烷气体并保持一定压力静置，实验装置不漏气，再关闭注气口阀门；
步骤五:恒温箱温度设置为1℃，整个实验过程中保持不变；
步骤六:将0.125-0.18mm和0.045-0.1mm不同粒度沉积物分层填入样品腔中，用带孔隔板将不同粒度沉积物层分隔开，气烟囱内放置0.25-0.5mm的沉积物，并按照实验设计在不同位置设置电阻率和温度传感器；
步骤七:加入盐度为35psu的NaCl溶液，沉积物孔隙空间被100％孔隙水饱和，上下左右摇晃在沉积物中均匀分布；
步骤八:通过轴压和围压控制系统分别控制样品腔的轴压和围压；
步骤九:静置1小时，并测试电阻率参数在期间的变化；
步骤十:设置注气气体流速为0-2000mL/min，精度为0.25F.S，气体增压上限为25MPa；向气烟囱注入甲烷气，开始水合物合成，并通过反应釜顶部排气系统控制反应釜内部压力值；
步骤十一：持续稳定注气，直至各电阻率传感器所测电阻率值保持不变或者气烟囱堵塞无法注气，实验结束；
步骤十二：收集并处理实验数据，利用Archie公式计算水合物饱和度；
利用Archie公式计算水合物饱和度，具体包括：
岩石孔隙100％被孔隙水饱和的沉积物中，电阻率R0与孔隙水电阻率Rw成正比，比例系数为地层因子F：
式中，a为初始沉积物岩性系数；是储层孔隙度；m为胶结系数；
含水合物沉积物中，电阻率Rt与100％含水的沉积物电阻率R0之比为电阻率增大指数I：
式中，b为含水合物沉积物岩性系数；Sw是储层含水饱和度；n是饱和度指数；
通过不同途径计算水合物饱和度Sh：
式中，a、b、m、n通过岩电实验确定，a和m与岩性和孔隙结构相关，b和n与岩石的润湿性和含油气性相关；b取值1，n取值1.94。
2.如权利要求1所述的气烟囱对天然气水合物生成和分布的影响分析方法，其特征在于，步骤六中，电阻率和温度传感器分别放置在气烟囱的上、中、下位置，以及粗细沉积地层中距离气烟囱的不同位置处，用来监测水合物的分布情况。
3.一种终端，其特征在于，所述终端至少搭载实现权利要求1～2任意一项所述气烟囱对天然气水合物生成和分布的影响分析方法的控制器。
4.一种计算机可读存储介质，包括指令，当其在计算机上运行时，使得计算机执行如权利要求1-2任意一项所述的气烟囱对天然气水合物生成和分布的影响分析方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811088537.2A CN108804860B (zh) | 2018-09-18 | 2018-09-18 | 气烟囱对天然气水合物生成和分布的影响分析系统及方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811088537.2A CN108804860B (zh) | 2018-09-18 | 2018-09-18 | 气烟囱对天然气水合物生成和分布的影响分析系统及方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108804860A CN108804860A (zh) | 2018-11-13 |
| CN108804860B true CN108804860B (zh) | 2023-11-24 |
Family
ID=64082079
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811088537.2A Active CN108804860B (zh) | 2018-09-18 | 2018-09-18 | 气烟囱对天然气水合物生成和分布的影响分析系统及方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108804860B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110988307A (zh) * | 2019-12-26 | 2020-04-10 | 中国科学院广州能源研究所 | 一种评价海洋天然气水合物储层优劣的方法 |
| CN113960672B (zh) * | 2021-10-22 | 2023-09-22 | 中国石油大学(华东) | 评价多边形断层在流体运移过程中作用的方法、系统、介质 |
| CN114733472A (zh) * | 2022-04-24 | 2022-07-12 | 中国科学院长春应用化学研究所 | 一种适用于多级反应的反应釜 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105403672A (zh) * | 2015-11-25 | 2016-03-16 | 中国科学院广州能源研究所 | 模拟天然气水合物开采过程地层形变的实验装置和方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9519072B2 (en) * | 2006-05-11 | 2016-12-13 | Schlumberger Technology Corporation | Method and apparatus for locating gas hydrate |
| EP3221022B8 (en) * | 2014-11-17 | 2020-03-04 | ExxonMobil Upstream Research Company | Distillation tower with a liquid collection system |
-
2018
- 2018-09-18 CN CN201811088537.2A patent/CN108804860B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105403672A (zh) * | 2015-11-25 | 2016-03-16 | 中国科学院广州能源研究所 | 模拟天然气水合物开采过程地层形变的实验装置和方法 |
Non-Patent Citations (2)
| Title |
|---|
| Probable patterns of gas flow and hydrate accretion at the base of the hydrate stability zone;Davies, Richard J.;GEOLOGY;第42卷(第12期);第1055-1058页 * |
| 天然气水合物成藏体系研究进展;卜庆涛等;新能源进展;第3卷(第6期);第435-443页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108804860A (zh) | 2018-11-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108804860B (zh) | 气烟囱对天然气水合物生成和分布的影响分析系统及方法 | |
| CN109307755B (zh) | 气藏水侵与排水采气物理模拟实验装置及方法 | |
| KR101350340B1 (ko) | 가스 하이드레이트 생산모사시스템 및 이를 이용한 생산모사방법 | |
| US9970267B2 (en) | Experimental device for simulating exploitation of natural gas hydrate in permeable boundary layer | |
| Delli et al. | Prediction performance of permeability models in gas-hydrate-bearing sands | |
| Shen et al. | Coupling effect of porosity and hydrate saturation on the permeability of methane hydrate-bearing sediments | |
| CN103471976B (zh) | 一种测量含水合物多孔沉积物渗透率的装置 | |
| WO2017088226A1 (zh) | 模拟天然气水合物开采过程地层形变的实验装置和方法 | |
| CN110847873B (zh) | 一种低渗天然气水合物藏原位水力射流开采装置及方法 | |
| AU2017243937B2 (en) | Device for determining petrophysical parameters of an underground formation | |
| WO2014176794A1 (zh) | 天然气水合物开采地层稳定性三维模拟装置 | |
| CN106896212B (zh) | 监测深水钻井液侵入过程水合物储层物性变化的装置 | |
| CN101376854A (zh) | 模拟三维条件下气体水合物成藏过程的方法及装置 | |
| CN107842341A (zh) | 一种天然气水合物开采监测模拟装置及方法 | |
| NO342878B1 (no) | Fremgangsmåte, program og datasystem for å harmonisere modelldata for hydrokarbonreservoarer | |
| CN111855715A (zh) | 一种盐下白云岩水岩反应的模拟与评价方法及系统 | |
| Li et al. | Permeability experiments on the methane hydrate in quartz sands and its model verification | |
| Choi et al. | An experimental system and procedure of unsteady-state relative permeability test for gas hydrate-bearing sediments | |
| CN110309611B (zh) | 基于气水厚度分布的气水两相渗流规律预测方法及系统 | |
| CN110905496B (zh) | 一种气驱超覆模拟装置及其使用方法 | |
| Li et al. | In situ estimation of relative permeability from resistivity measurements | |
| CN115468824A (zh) | 一种气体水合物制备及其物性参数原位检测的装置和方法 | |
| Goda et al. | Wettability quantification–prediction of wettability for australian formations | |
| CN107288633A (zh) | 具有重力分异现象的凝析气藏地质储量计算方法 | |
| CN208834300U (zh) | 气烟囱对天然气水合物生成和分布的影响分析装置 |
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 |