CN106093299B - 一种致密气储层钻井液伤害评价实验方法 - Google Patents
一种致密气储层钻井液伤害评价实验方法 Download PDFInfo
- Publication number
- CN106093299B CN106093299B CN201610390509.0A CN201610390509A CN106093299B CN 106093299 B CN106093299 B CN 106093299B CN 201610390509 A CN201610390509 A CN 201610390509A CN 106093299 B CN106093299 B CN 106093299B
- Authority
- CN
- China
- Prior art keywords
- rock core
- drilling fluid
- core
- rock
- sample
- 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.)
- Expired - Fee Related
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 148
- 238000005553 drilling Methods 0.000 title claims abstract description 145
- 230000006378 damage Effects 0.000 title claims abstract description 129
- 238000011156 evaluation Methods 0.000 title claims abstract description 53
- 238000002474 experimental method Methods 0.000 title claims abstract description 34
- 239000011435 rock Substances 0.000 claims abstract description 246
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 105
- 238000012360 testing method Methods 0.000 claims abstract description 94
- 239000000706 filtrate Substances 0.000 claims abstract description 82
- 230000035699 permeability Effects 0.000 claims abstract description 80
- 208000027418 Wounds and injury Diseases 0.000 claims abstract description 68
- 208000014674 injury Diseases 0.000 claims abstract description 67
- 230000009545 invasion Effects 0.000 claims abstract description 65
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 55
- 238000003384 imaging method Methods 0.000 claims abstract description 39
- 239000008398 formation water Substances 0.000 claims abstract description 25
- 239000007788 liquid Substances 0.000 claims abstract description 18
- 238000009826 distribution Methods 0.000 claims abstract description 17
- 230000011218 segmentation Effects 0.000 claims abstract description 15
- 238000001228 spectrum Methods 0.000 claims abstract description 5
- 238000013215 result calculation Methods 0.000 claims abstract description 3
- 238000000034 method Methods 0.000 claims description 60
- 239000007789 gas Substances 0.000 claims description 51
- 238000001035 drying Methods 0.000 claims description 34
- 239000007790 solid phase Substances 0.000 claims description 14
- 230000005311 nuclear magnetism Effects 0.000 claims description 13
- 238000004458 analytical method Methods 0.000 claims description 12
- 238000006073 displacement reaction Methods 0.000 claims description 12
- 238000005516 engineering process Methods 0.000 claims description 12
- 238000005119 centrifugation Methods 0.000 claims description 11
- 238000004088 simulation Methods 0.000 claims description 11
- 239000007791 liquid phase Substances 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 10
- OKKJLVBELUTLKV-UHFFFAOYSA-N Methanol Chemical compound OC OKKJLVBELUTLKV-UHFFFAOYSA-N 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 9
- 229920006395 saturated elastomer Polymers 0.000 claims description 9
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 239000004576 sand Substances 0.000 claims description 7
- 239000011555 saturated liquid Substances 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 5
- 238000001225 nuclear magnetic resonance method Methods 0.000 claims description 5
- 238000005303 weighing Methods 0.000 claims description 5
- 239000003153 chemical reaction reagent Substances 0.000 claims description 4
- 238000000703 high-speed centrifugation Methods 0.000 claims description 4
- 229910052757 nitrogen Inorganic materials 0.000 claims description 4
- 239000012071 phase Substances 0.000 claims description 4
- 239000012047 saturated solution Substances 0.000 claims description 4
- 238000003860 storage Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000003825 pressing Methods 0.000 claims description 3
- 150000003839 salts Chemical class 0.000 claims description 3
- 238000005520 cutting process Methods 0.000 description 12
- 238000011160 research Methods 0.000 description 7
- 239000002245 particle Substances 0.000 description 6
- 230000000903 blocking effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 239000003795 chemical substances by application Substances 0.000 description 4
- 239000003208 petroleum Substances 0.000 description 3
- 241001074085 Scophthalmus aquosus Species 0.000 description 2
- CDBYLPFSWZWCQE-UHFFFAOYSA-L Sodium Carbonate Chemical compound [Na+].[Na+].[O-]C([O-])=O CDBYLPFSWZWCQE-UHFFFAOYSA-L 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000008595 infiltration Effects 0.000 description 2
- 238000001764 infiltration Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- HEMHJVSKTPXQMS-UHFFFAOYSA-M Sodium hydroxide Chemical class [OH-].[Na+] HEMHJVSKTPXQMS-UHFFFAOYSA-M 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000007605 air drying Methods 0.000 description 1
- 239000003513 alkali Substances 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 238000009388 chemical precipitation Methods 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 239000002734 clay mineral Substances 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000013100 final test Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000036571 hydration Effects 0.000 description 1
- 238000006703 hydration reaction Methods 0.000 description 1
- 239000010977 jade Substances 0.000 description 1
- 229920005610 lignin Polymers 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000009527 percussion Methods 0.000 description 1
- 229920002401 polyacrylamide Polymers 0.000 description 1
- 238000004886 process control Methods 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 235000017550 sodium carbonate Nutrition 0.000 description 1
- 229910000029 sodium carbonate Inorganic materials 0.000 description 1
- 230000008961 swelling Effects 0.000 description 1
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
- G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
-
- 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/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N2015/0813—Measuring intrusion, e.g. of mercury
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Pathology (AREA)
- Biochemistry (AREA)
- Engineering & Computer Science (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Geochemistry & Mineralogy (AREA)
- Remote Sensing (AREA)
- Dispersion Chemistry (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
本发明公开了一种致密气储层钻井液伤害评价实验方法，其步骤为：岩心样品烘干，做标记，进行均匀分段切片处理，对岩心进行干样重量等测试；将各小段岩心放入流体中浸润饱和；对岩心进行换向离心操作，建立起各小段岩心样品原始含水饱和度Swi；对模拟地层水和钻井液滤液的核磁共振T2谱进行测试并判断差异；对岩心的内部流体分布状态做核磁成像测试；将岩心放入钻井液循环仪的岩心夹持器进行钻井液循环后测试并计算各小段岩心的含水饱和度Swi和气体渗透率K2，用测试结果计算对比钻井液侵入后造成的渗透率伤害率，评价钻井液侵入导致的不同深度段的伤害程度。本发明是对原有的切片法进行改进，建立了新的适用于致密、超致密储层钻井液伤害评价指标。
Description
技术领域
本发明涉及油气田开发领域，尤其涉及一种致密砂岩气藏钻井过程中钻井液伤害评价实验方法。
背景技术
在正压差钻井过程中，钻井液中的固相颗粒和液相不可避免地会在压差作用下侵入储层，改变井筒附近储层流体分布的原始状态，并可能导致储层伤害。潜在的伤害主要包括四种类型：1)滤液对储层的伤害：在井内液柱压力与储层压力差的作用下，滤液进入储层。滤液能使储层内的粘土发生水化膨胀，或使储层中粘土矿物等微粒发生分散运移堵塞储层孔道而产生伤害；滤液进入储层内产生水锁效应堵塞储层而产生伤害；滤液的组成与地层水不配伍时，产生各种化学沉淀堵塞储层。在井壁上未形成泥饼之前，滤液对油、储层的伤害更为严重；2)钻井液的固相颗粒对储层的伤害：固相颗粒在井壁上未形成泥饼之前，侵入储层孔道堵塞使储层的渗透率大大降低；3)压差对储层的伤害：通常钻井过程中钻井液柱压力大于地层压力，因此，压差是客观存在的，对储层的伤害也是必然的。压差越大，滤液固相颗粒侵入油、储层越深，污染越严重，伤害程度越大；4)钻井液处理剂对储层的伤害：常规钻井液处理剂中的烧碱、纯碱、碱溶性有机处理剂、木质素衍生物、维纤素衍生物、腐植酸及其衍生物、聚丙烯酰胺及其衍生物等，都能改善钻井液的性能，提高钻井速度，稳定井眼，保证井下安全。但各种处理剂都会不同程度的伤害储层。
钻井液侵入伤害一旦形成将影响油气藏的开发。如何准确定量的评价伤害程度，为钻井工程提供合理的施工建议，是油气勘探开发过程中面临的一个难题。一般通过研究钻井液对油气层滤失量、侵入深度、渗透率伤害程度的大小，来客观反映钻井液中对油气层的伤害程度。
由于钻井液伤害涉及内容多且复杂，对于钻井液造成储层伤害的室内研究，早期不同研究者的研究结果都是根据各自的设备和实验方法得到的，相互之间可比性很差，且同一实验的再现性差。虽然在2002年国内已经形成了一套钻井液完井液伤害油层室内评价方法行业标准，但该标准与其他研究者的评价过程中仍具有一定局限性。
目前钻井液伤害程度室内试验评价的方法主要有两种：
第一种是世界范围内普遍采用的岩心伤害率(DR)作为储层伤害的室内试验评价指标。实验用整根岩心污染前后渗透率的比值来描述钻井液对储层的伤害，而伤害深度的评价则可以通过分段截割法和电阻率测试法实现。
分段截割法是在钻井液伤害试验后，对测试长岩心进行连续截割，测量剩余段的渗透率，直到渗透率为一常数。由岩心的总长度减去渗透率为一常数的剩余长度，得到固相侵入岩心的深度。1989年Vol13，石油大学学报(自然科学版)，夏俭英在《研究泥浆固相侵入岩心深度的分段截割法》中分别用钻井液和滤液污染岩心。在污染前后分别测量其渗透率，然后对比分析，得出固相和液相引起岩心伤害的程度。采用分段截割法连续截割岩心，评价钻井液固相颗粒侵入岩心的深度。
电阻率测试法，依据饱和流体岩石电阻率为定值，随之而流体侵入岩石电阻率发生变化。通过实验装置在实验岩心轴向方向上布设电极。通过计算机自动采集数据实现连续监测整个侵入过程中的电阻率变化，确定侵入深度。通过测试污染前后渗透率确定污染程度。1994年Vol16，西南石油学院学报，张开洪在《泥浆滤液侵入对岩石物性及电性影响的实验研究》，1996年Vol18，西南石油学院学报，陈可贵等在《长岩心实验确定地层侵入深度和污染深度的研究》，1999年Vol23，测井技术，胡国恒等在《钻井液侵入对储层电性物性影响实验研究》中都依据此方法通过连续监测侵入过程中的电阻率变化,并确定侵入深度，利用渗透率变化评价了伤害程度。
第二种是Marx(1987)提出用各段岩心伤害率(SDR)确定伤害深度的方法。包括切片法、分段截割法和梯度法。
切片法是把被污染岩心切为几段，分别测定各段岩心污染前后的渗透率，由渗透率恢复率确定伤害程度。2012年，陈光在硕士学位论文《钻井液伤害储层室内评价实验》中采用改进的切片法，测试污染前后各段岩心渗透率测，研究钻井液在一定压力和温度下对地层的污染规律。
梯度法是通过测定沿岩心轴向的压力降来确定各段岩心在钻井液侵入过程中的渗透率，明确伤害深度和程度。1998年Vol21，钻采工艺，王松等在《高温高压钻井液完井液对储层伤害的室内研究》中利用高温高压动失水仪和多点液体渗透率梯度测定仪，评价了常规处理剂对人造岩心的伤害程度。1999年Vol24，地球科学—中国地质大学学报，张琰等在《储层伤害室内评价方法的研究》中修正了油气层伤害室内评价指标岩心伤害率DR的定义，用污染后伤害段岩心的渗透率取代污染后的岩心渗透率，作为岩心伤害程度的评价参数。分别利用切片法和改进的梯度测量法确定岩心伤害率。
此外还有学者将核磁共振技术引入到钻井液伤害的评价中，2003年Vol20，钻井液与完井液，范振中等《核磁共振技术用于评价钻井液对岩心的伤害》利用核磁共振技术测量污染前后的孔道直径变化，分析三维空间内岩心的污染程度。但仅通过信号幅度识别孔径大小变化趋势，并无伤害程度的定量分析。
2002年国内形成了首个关于钻井液伤害评价的行业标准SY/T6540-2002《钻井液完井液伤害油层室内评价方法》，2009年Vol26，钻井液与完井液，孟小海在《钻井液与水泥浆伤害油层室内综合动态评价方法》中依据标准，对滤失量、侵入深度、渗透率伤害和恢复程度等参数进行了分析。指出根据岩心的性质，采用不同的方法确定岩心污染状况。
上述的室内评价实验方法与标准，虽然在一定程度上定量评价了钻井液对储层伤害程度与伤害深度，但随着实践的不断展开与分析精度要求的提高，各种实验方法的局限性也显现出来：
采用岩心伤害率(DR)作为伤害指标的方法，实验用整根岩心污染前后渗透率的比值反映钻井液对储层的伤害程度，岩心长度无统一标准。这种方法认为整根岩心都受到了伤害，因此评价不能真实地反映储层伤害程度，此外岩心长度的差异会对评价结果有较大影响；分段截割确定固相侵入深度时会造成岩心厚度损失和破坏，实验重复性和再现性差，且岩心非均一性强，截割剩余岩样渗透率测量趋于常数过程控制难度较大；电阻率法确定侵入深度，测量得到的电阻率和渗透率是两测点间岩心段的平均值，对于非均匀侵入的岩心，由于截面上的值一般与测量平均值不等。实验得到的侵入深度和伤害污染程度难以反映客观事实。
用各段岩心伤害率(SDR)确定伤害深度的方法中切片法需要破坏岩心，且污染前各段岩心渗透率只能看作均一的，各段的伤害程度评价不真实；梯度法虽然不需要破坏岩心，还可以对污染过程进行实时监测，但实验对压力测量精度要求很高。
此外上述实验中还具有一些共性问题：钻井液伤害渗透率变化初始值选择的是干岩心或者岩心饱和束缚水状态时的渗透率，而非地层原始含水饱和度Swi状态渗透率，评价结果不能反映地层条件下的伤害程度；多数实验采用的低压岩心饱和地层水对于致密砂岩并不适用；建立岩心含水饱和度的风干、烘干、气驱法，难以保障水在岩心中均匀分布，对钻井液侵入伤害结果的评价具有影响；对于致密、超致密砂岩钻井液侵入损害评价，往往因为损害后局部渗透率降低较大而引起渗透率测试困难，形成岩心整体损害程度强的假象，且测试时间过长，岩心样品的含水饱和度会发生较大的变化，影响测试结果；实验评价指标主要为渗透率变化率，无法明确钻井液液相、固相的侵入深度。
发明内容
针对上述问题，为克服现有技术的不足，本发明的目的是提供一种致密气储层钻井液伤害评价的实验方法，该方法考虑了长岩心整体评价的不客观，分段截割的破坏性，电阻率测量的评价方法的误差性，对原有的切片法进行了改进，同时考虑到致密气储层砂岩低孔渗特性，在岩样处理、岩心含水饱和度建立、钻井液侵入、钻井液滤液侵入深度评价等方面做出了创新性发明，建立了新的适用于致密、超致密储层钻井液伤害评价指标。
一种致密气储层钻井液伤害评价实验方法，包含以下步骤：
S1、岩心样品烘干，做标记，进行均匀分段切片处理；对分段切片的各段岩心进行干样重量G0、长度L、直径D、气体渗透率K0、气体孔隙度测试；
S2、在高压真空环境下，将各小段岩心放入密度为ρ1的流体中浸润饱和，测试各段岩心饱和水后重量G1，计算各段岩心表观体积Vt，有效孔隙体积Vp和岩样有效孔隙度
S3、利用高速离心和核磁共振技术，对岩心进行换向离心操作，使其含水率达到设定值，建立起各小段岩心样品原始含水饱和度Swi；
S4、对模拟地层水和钻井液滤液的核磁共振T2时间进行测试并判断两者之间的差异；利用核磁成像技术对各小段岩心在原始含水饱和度Swi状态的内部流体分布状态做核磁成像测试，并对测试图像进行伪彩处理；测试各小段岩心的原始含水饱和度Swi状态下的气体渗透率K1和样品重量G2，利用称重法校正含水饱和度；
S5、将各小段岩心按照切开前的顺序放入钻井液循环仪的岩心夹持器，利用密度为的ρ2钻井液进行钻井液循环；循环后，计算各小段岩心的含水饱和度Sqwi，测试气体渗透率K2和各小段岩心重量G3，利用称重法校正核磁测试含水饱和度，依据渗透率变化率公式Ir＝(K1-K2)/K1×100％，用测试结果计算对比各小段岩心钻井液侵入后造成的渗透率伤害率Ir，评价钻井液侵入导致的不同深度段的伤害程度；
S6、对钻井液侵入模拟的各段岩心进行T2加权核磁成像测试，测试参数与步骤S4中相同；将钻井液侵入后的各段岩心核磁成像的图像与原始含水饱和度Swi状态核磁成像测试图像对比；分析钻井液侵入前后各段岩心内部流体分布变化，明确钻井液液相侵入深度Ll；
S7、岩心烘干处理，对各小段岩心进行扫描电镜分析，以此确定钻井液固相侵入深度Ls，分析岩心内部微观特征差异。
进一步的，所述步骤S1具体包括如下内容：
在井筒取芯的钻井液侵入方向上钻取柱塞样品，保持岩样端面和柱面平整，端面垂直柱面，柱塞结构完整；柱塞样品直径为2.5cm，柱塞样品长度不小于直径的1.5倍；
对岩样采用甲醇试剂进行除盐处理，处理后进行岩样烘干，烘干温度控制在60～65℃，相对湿度控制在40％～50％；岩样烘干至恒重，烘干时间为48h，48h后每2h称量一次，两次称量的差值小于10mg；
将烘干的岩心沿岩心轴向在侧面画一个由左向右的箭头指示方向，将岩心均匀切成n等份，每段长度1.4～1.6cm，由左至右分别标号，标号的数值为1、2、…、n；
分别称取各小段岩心的干样重G0n，测量各小段岩心的长度Ln直径Rn，n的数值为1、2、…、n，各小段岩心干样测试气体渗透率K0n和气体孔隙度n的数值为1、2、…、n，测试围压3.5MPa，驱替压差1MPa，根据实验岩心的渗透率特征，如果渗透率过低，测试压差可相应增加，并保证围压是驱替压力的2倍以上。
进一步的，所述步骤S2具体包括如下内容：
将烘干后的岩样抽真空20h以上，将密度为ρ1的模拟地层水放入进行饱和；
采取25MPa高压作为饱和压力，加压时间不低于20h，保证致密砂岩岩样充分饱和模拟地层水；
岩样在饱和液中浸泡至少24h以上，测定饱和液体后各小段岩心质量G1n，n的数值为1、2、…、n；
计算岩样的有效孔隙体积与有效孔隙度：
Vp＝(G1–G0)/ρ1；
Vt＝L×A；
A＝π×(D/2)2；
G0-干岩样质量，g；
G1-岩心饱和水后重量，g；
ρ1-测试温度下的饱和液体的密度，g/cm3；
Vρ-岩样有效孔隙体积，cm3；
Vt-岩心表观体积，cm3；
L-样品长度，cm；
D-样品直径，cm；
A-样品端面积，cm2；
-岩样有效孔隙度，％；
进一步的，所述步骤S3具体包括如下内容：
利用高速离心机对饱和模拟地层水的岩心进行离心操作，依次在100psi、200psi、300psi、400psi、500psi、600psi离心力下进行恒速离心，离心过程中每8min调换岩心方向，离心时间40min；在每个离心速度下完成离心操作后，通过核磁T2谱测试计算含水饱和度，当测试值接近岩心原始含水饱和度Swi后，不再进行更高一级的离心速度离心；
对核磁建立的原始含水饱和度Swi状态的各段岩心称重G2n，n的数值为1、2、…、n，应用称重法，根据Swi＝{1-[(G2-G0)/ρ1]/Vp}×100％计算含水饱和度，校正核磁共振法建立的原始含水饱和度Swi；
如果储层原始含水饱和度Swi小于30％，还需联合N2气驱替、冰箱冷藏方式建立含水饱和度。
进一步的，所述步骤S4具体包括如下内容：
分别对10mL模拟地层水和10mL钻井液滤液进行核磁共振T2谱测试，获得模拟地层水的T2时间：T2,D-H2O，和钻井液滤液T2时间：T2,Z-H2O，明确两种流体的T2时间差异；其测试参数包括等待时间TW—2500ms和回波时间TE—100ms；
在明确两种流体的核磁共振T2时间差异后，利用核磁成像测试各小段岩心原始含水饱和度Swi状态的内部流体分布状态，为便于对比钻井液侵入后的流体分布，需要选取钻井液滤液信号较强的测试参数进行测量，选取等待时间TW为500ms、回波时间TE为10ms；对成像获取的灰度图像进行伪彩处理，通过图片中亮度分析流体分布；
测试各段岩心的原始含水饱和度Swi状态下的气体渗透率K1n，n的数值为1、2、…、n，测试介质为N2气，测试围压3.5MPa，驱替压差1MPa，对于离心法建立的岩心含水饱和度，如果与地层原始含水饱和度相差±2％，可以认为所测渗透率即为地层状态下岩心渗透率，若相差较大，则通过测量的渗透率通过克里金插值计算岩心原始含水饱和度对应的渗透率值K1n，n的数值为1、2、…、n。如果测试时间超过30min，则须用称重法校正含水饱和度。
进一步的，所述步骤S5具体包括如下内容：
将各小段岩心按照切开前的顺序排列，依次放入钻井液循环仪的岩心夹持器，岩心标记的箭头指向钻井液侵入的方向，设置围压10MPa；
将井场取回的密度为ρ2钻井液盛入钻井液循环仪釜体，在循环压差△P＝3.5MPa条件下循环，也可以根据岩心样品的差异，根据地层压力系数与钻井液密度差异计算相应的循环压差。循环时间根据现场钻井时长选择，为评价短时间和长时间钻井液对气层的伤害，选择3小时和96小时两个实验时间；实验过程中记录不同时段对应的钻井液滤失量；
钻井液循环后，利用核磁共振法测试每小段岩心的含水饱和度Sqwin，n的数值为1、2、…、n，并称取各小段岩心钻井液伤害后的重量G3n，n的数值为1、2、…、n，用称重法，根据Sqwi＝{[(G3-G2)/ρ2]/Vp}×100％+Swi校正核磁测试含水饱和度；
对各段岩心进行气体渗透率测试，K2n，n的数值为1、2、…、n；测试介质为氮气，测试围压3.5MPa，驱替压差1MPa；依据渗透率变化率公式Ir＝(K1-K2)/K1×100％，计算对比各小段岩心钻井液侵入后造成的渗透率伤害率Irn，n的数值为1、2、…、n；并依据渗透率伤害率Irn评价钻井液侵入导致的不同深度段的伤害程度，其评价标准如下：钻井液伤害率为Irn≤5％时伤害程度为无伤害，钻井液伤害率为5％<Irn≤30％时伤害程度为弱，钻井液伤害率为30％<Irn≤50％时伤害程度为中等偏弱，钻井液伤害率为50％<Irn≤70％时伤害程度为中等偏强，钻井液伤害率为70％<Irn≤90％时伤害程度为强，钻井液伤害率为90％<Irn时伤害程度为极强。
进一步的，所述步骤S6具体包括如下内容：
对钻井液侵入模拟的各段岩心进行T2加权核磁成像测试，即测试参数具有以下特点：等待时间小于模拟地层水T2时间，Tw＜T2,D-H2O，回波时间TE尽可能小，这样有利于压制岩心中模拟地层水的信号，使模拟地层水和钻井液滤液信号有效区分；等待时间TW—500ms、回波时间TE—10ms；
对核磁成像的灰度图像进行伪彩处理，彩色区域为水信号，颜色较暗的为模拟地层水，颜色亮的信号为钻井液滤液；
将侵入后的各段岩心核磁成像的图像与原始含水饱和度Swi状态所得到的核磁成像测试图像对比；分析钻井液侵入前后各段岩心内部流体分布变化，明确钻井液液相侵入深度Ll。
进一步的，所述步骤S7具体包括如下内容：
对实验岩心进行烘干，烘干温度控制在60～65℃，相对湿度控制在40％～50％；岩样烘干至恒重，烘干时间48h，48h后每2h称量一次，两次称量的差值小于10mg；
对按岩心切片前整体箭头标记方向排列，由左至右，左侧为钻井液侵入方向，将各小段岩心沿轴向中线剖开，对各段岩心剖开面进行连续扫描电镜分析，确定钻井液固相侵入深度Ls，分析岩心内部微观特征差异。
本发明的优点是：
(1)该方法考虑了长岩心整体评价的不客观、分段截割的破坏性、电阻率测量的评价方法的误差性等问题，对原有的切片法进行了改进；
(2)考虑到致密气储层砂岩低孔渗特性，在岩样处理、岩心含水饱和度建立、钻井液侵入、钻井液滤液侵入深度评价等方面做出了创造性的改进，建立了新的适用于致密、超致密储层钻井液伤害评价指标，比现有技术方案更加高效和准确。
附图说明
图1是本发明实施例提供的15号样品三段切片示意图；
图2是地层水、钻井液、钻井液滤液核磁共振T2谱图；
图3是15号样品初始地层含水与钻井液侵入后的核磁共振T2谱测试图；
图4是15号样品原始含水饱和度Swi时砂岩水分布样品核磁成像测试图；
图5是15号样品钻井液侵入后Sqwi时砂岩内部水分布核磁成像测试图；
图6是15号样品钻井液侵入后Sqwi时各小段岩心连续扫描电镜分析图。
具体实施方式
下面结合附图中的实施例对本发明作进一步的详细说明，但并不构成对本发明的任何限制。
下述实施例中所使用的实验方法如无特殊说明，均为常规方法。
下述实施例中所用的材料、试剂等，如无特殊说明，均可从商业途径得到。
实施例：
以克拉苏气田致密砂岩岩心为例，根据试井资料，储层原始含水饱和度Swi为25～35％。选取6块致密砂岩岩心用该发明方法评价钻井液侵入伤害程度。此处采用15号岩心样本作为实验对象，采用本发明提供的致密气储层钻井液伤害评价的实验方法，包括以下步骤：
如图1所示，步骤1、将岩心样品烘干，做标记，进行均匀分段切片处理。对分段切片的各段岩心的进行干样重量G0、长度L、直径D、气体渗透率K0、气体孔隙度测试。
步骤2、各小段岩心高压真空饱和流体。
步骤3、利用高速离心、核磁共振技术，通过岩心换向离心建立各小段岩心样品原始含水饱和度Swi。
步骤4、判断模拟地层水和钻井液滤液的T2时间差异。选取钻井液滤液信号较强的测试参数对各小段岩心原始含水饱和度Swi状态的内部流体分布状态核磁成像测试。测试各小段岩心的原始含水饱和度Swi状态下的气体渗透率K1。
步骤5、将各小段岩心按照切开前的顺序放入钻井液循环仪的岩心夹持器，进行钻井液循环。循环后后，计算每小段岩心的含水饱和度Sqwi，测试气体渗透率K2。依据渗透率变化率公式Ir＝(K1-K2)/K1×100％，计算对比各小段岩心钻井液侵入后造成的渗透率伤害率Ir，评价钻井液侵入导致的不同深度段的伤害程度。
步骤6、对钻井液侵入模拟的各段岩心进行T2加权核磁成像测试，测试参数与步骤4中相同。对核磁成像的灰度图像做伪彩处理，将侵入后的各段岩心核磁成像的图像与原始含水饱和度Swi状态核磁成像测试图像对比。分析钻井液侵入前后各段岩心内部流体分布变化，明确钻井液液相侵入深度Ll。
步骤7、岩心烘干处理、对各小段岩心进行扫描电镜分析，确定钻井液固相侵入深度Ls，分析岩心内部微观特征差异。
如图1所示，步骤1岩心样品烘干，做标记，进行均匀分段切片处理包括以下步骤：
步骤11、在井筒取芯的钻井液侵入方向上钻取柱塞样品，岩样端面、柱面保持平整，端面垂直柱面，柱塞应具有结构完整性；直径一般为2.5cm，长度不小于直径的1.5倍。
步骤12、对岩样采用甲醇等试剂进行除盐处理，处理后进行岩样烘干，烘干温度控制在60～65℃，相对湿度控制在40％～50％；岩样烘干至恒重，烘干时间48h，48h后每2h称量一次，两次称量的差值小于10mg。
步骤13、将烘干的岩心沿岩心轴向在侧面画一个由左向右的箭头指示方向(如图1中箭头方向)，将岩心均匀切成n等份，每段长度1.4～1.6cm，由左至右分别标号n(n＝1、2、……)；
步骤14、分别称取各小段岩心的干样重G0n，测量各小段岩心的长度Ln直径Rn(n＝1、2、……)、各小段岩心干样测试气体渗透率K0n和气体孔隙度测试围压3.5MPa，驱替压差1MPa，根据实验岩心的渗透率特征，如果渗透率过低，测试压差可相应增加，并保证围压是驱替压力的2倍以上。
如图3所示，步骤2各小段岩心高压真空饱和流体包括以下步骤。
步骤21、将烘干后的岩样抽真空20h以上，将密度为ρ1的模拟地层水放入进行饱和；
步骤22、采取25MPa高压作为饱和压力，加压时间不低于20h，保证致密砂岩岩样充分饱和模拟地层水。
步骤23、岩样在饱和液中浸泡至少24h以上，测定饱和液体后各小段岩心质量G1n(n＝1、2、……)。
步骤24、计算岩样的有效孔隙体积和有效孔隙度：
Vp＝(G1–G0)/ρ1
Vt＝L×A
A＝π×(D/2)2
G0-干岩样质量，g；
G1-岩心饱和水后重量，g；
ρ1-测定温度下饱和液体的密度，g/cm3；
Vρ-岩样有效孔隙体积，cm3；
Vt-岩心表观体积，cm3；
L-样品长度，cm；
D-样品直径，cm；
A-样品端面积，cm2；
-岩样有效孔隙度，％；
如图2所示，步骤3利用高速离心、核磁共振技术，通过岩心换向离心建立各小段岩心样品原始含水饱和度Swi包括以下步骤。
步骤31、利用高速离心机对饱水岩心离心，依次在100psi、200psi、300psi、400psi、500psi、600psi的离心力下，进行恒速离心，离心过程中每8min调换岩心方向，离心时间40min。每个离心速度下离心后，通过核磁T2谱测试计算含水饱和度，直到接近岩心原始含水饱和度Swi后不在进行高一级的离心速度离心。
步骤32、对核磁建立的原始含水饱和度Swi状态的各段岩心称重G2n(n＝1、2、……)，应用称重法，根据Swi＝{1-[(G2-G0)/ρ1]/Vp}×100％计算含水饱和度，校正核磁共振法建立的原始含水饱和度Swi。
G0-干岩样质量，g；
G2-岩样原始含水饱和度Swi时的质量，g；
ρ1-在测定饱和温度下饱和液体的密度，g/cm3；
Vρ-岩样有效孔隙体积，cm3；
步骤33、如果储层原始含水饱和度Swi小于30％，还需结合N2气驱替、冰箱冷藏等联合作用建立含水饱和度。
步骤4判断模拟地层水和钻井液滤液的T2时间差异。核磁成像测试各小段岩心原始含水饱和度Swi状态的内部流体分布状态。核磁共振测试各小段岩心的原始含水饱和度Swi状态下的气体渗透率K1包括以下步骤。
步骤41、分别对10mL模拟地层水和10mL钻井液滤液进行核磁共振T2谱测试，获得模拟地层水T2时间T2,D-H2O和钻井液滤液T2时间T2,Z-H2O，明确两种流体的T2时间差异。测试参数：等待时间TW—2500ms、回波时间TE—100ms。
步骤42、在明确两种流体的核磁共振T2时间差异后，利用核磁成像测试各小段岩心原始含水饱和度Swi状态的内部流体分布状态，为便于对比钻井液侵入后的流体分布，选取钻井液滤液信号较强的测试参数：等待时间TW—500ms、回波时间TE—10ms。对成像获取的灰度图像进行伪彩处理，通过图片中亮度分析流体分布。
步骤43、测试各段岩心的原始含水饱和度Swi状态下的气体渗透率K1n(n＝1、2、……)，测试介质氮气，测试围压3.5MPa，驱替压差1MPa。对于离心法建立的岩心含水饱和度，如果与地层原始含水饱和度相差±2％，可以认为所测渗透率即为地层状态下岩心渗透率，如果相差较大，则通过测量的渗透率通过克里金插值计算岩心原始含水饱和度对应的渗透率值。如果测试时间超过30min，则须用称重法校正含水饱和度。
如图3所示，步骤5将各小段岩心按照切开前的顺序放入钻井液循环仪的岩心夹持器，进行钻井液循环。循环后后，测试、计算每小段岩心的含水饱和度Sqwi、气体渗透率K2。依据渗透率变化率公式Ir＝(K1-K2)/K1×100％，计算对比各小段岩心钻井液侵入后造成的渗透率伤害率Ir，评价钻井液侵入导致的不同深度段的伤害程度包括以下步骤。
步骤51、将各小段岩心按照切开前的顺序排列，依次放入钻井液循环仪的岩心夹持器，岩心标记的箭头指向钻井液侵入的方向，设置围压10MPa。
步骤52、将井场取回的密度ρ2的钻井液盛入钻井液循环仪釜体，在循环压差△P＝3.5MPa条件下循环，也可以根据岩心样品的差异和地层压力系数与钻井液密度差异计算相应的循环压差。循环时间根据现场钻井时长选择，为评价短时间和长时间钻井液对气层的伤害，选择3小时和96小时两个实验时间。实验过程中记录不同时段对应的钻井液滤失量。
步骤53、钻井液循环后，利用核磁共振法测试每小段岩心的含水饱和度Sqwin(n＝1、2、……)，并称取各小段岩心钻井液伤害后的重量G3n(n＝1、2、……)，用称重法，根据Sqwi＝{[(G3-G2)/ρ2]/Vp}×100％+Swi校正核磁测试含水饱和度。
Sqwi-钻井液侵入后含水饱和度，％；
G2-岩样原始含水饱和度Swi时的质量，g；
G3-钻井液侵入后岩样的质量，g；
ρ2-测定温度下钻井液滤液的密度，g/cm3；
Vρ-岩样有效孔隙体积，cm3；
步骤54、对各段岩心进行气体渗透率测试，K2n(n＝1、2、……)。测试介质为氮气，测试围压3.5MPa，驱替压差1MPa。依据渗透率变化率公式Ir＝(K1-K2)/K1×100％，计算对比各小段岩心钻井液侵入后造成的渗透率伤害率Irn，n的数值为1、2、…、n。并依据渗透率伤害率Irn评价钻井液侵入导致的不同深度段的伤害程度。
Ir-渗透率变化率，％；
K1-岩心原始含水饱和度Swi状态下气体渗透率，mD；
K2-钻井液侵入后岩心气体渗透率，mD；
表1钻井液伤害率与伤害程度评价关系表
| 钻井液伤害率％ | 伤害程度评价 |
| Ir≤5 | 无 |
| 5<Ir≤30 | 弱 |
| 30<Ir≤50 | 中等偏弱 |
| 50<Ir≤70 | 中等偏强 |
| 70<Ir≤90 | 强 |
| 90<Ir | 极强 |
如图4和图5所示，步骤6对钻井液侵入模拟的各段岩心进行T2加权核磁成像测试，对核磁成像的灰度图像进行伪彩处理，将侵入后的各段岩心核磁成像的图像与原始含水饱和度Swi状态核磁成像测试图像对比。分析钻井液侵入前后各段岩心内部流体分布变化，明确钻井液液相侵入深度包括以下步骤。
步骤61、对钻井液侵入模拟的各段岩心进行T2加权核磁成像测试，即测试参数具有以下特点：等待时间小于模拟地层水T2时间，Tw＜T2,D-H2O，回波时间TE尽可能小，这样有利于压制岩心中模拟地层水的信号，使模拟地层水和钻井液滤液信号有效区分。等待时间TW—500ms、回波时间TE—10ms。
步骤62、对核磁成像的灰度图像进行伪彩处理，彩色区域为水信号，颜色较暗的为模拟地层水，颜色亮的信号为钻井液滤液。
步骤63、将侵入后的各段岩心核磁成像的图像与原始含水饱和度Swi状态核磁成像测试图像对比。分析钻井液侵入前后各段岩心内部流体分布变化，明确钻井液液相侵入深度Ll。
如图6所示，步骤7岩心烘干处理、对各小段岩心进行扫描电镜分析，对钻井液侵入后的微观孔隙结构进行观察包括以下步骤。
步骤71、对实验岩心进行烘干，烘干温度控制在62℃，相对湿度控制在40％～50％；岩样烘干至恒重，烘干时间48h，48h后每2h称量一次，两次称量的差值小于10mg。
步骤72、对岩心切片前，整体按箭头标记方向排列，由左至右，左侧为钻井液侵入方向,将各小段岩心沿轴向中线剖开，对各段岩心剖开面进行连续扫描电镜分析，确定钻井液固相侵入深度Ls，分析岩心内部微观特征差异。
通过上述步骤，最终得到的15号岩心样本钻井液伤害评价结果，其他岩心样本也如同样的方法进行测试，最终测试结果如下表所示：
表2克深地区致密砂岩钻井液伤害评价
以上所述仅是本发明的优选实施方式，应当指出，本发明并不局限于上述方式，在不脱离本发明原理的前提下，还能进一步改进，这些改进也应视为本发明的保护范围。
Claims (8)
1.一种致密气储层钻井液伤害评价实验方法，其特征在于，包含以下步骤：
S1、岩心样品烘干，做标记，进行均匀分段切片处理；对分段切片的各段岩心进行干样重量G0、长度L、直径D、气体渗透率K0、气体孔隙度进行测试；
S2、在高压真空环境下，将各小段岩心放入密度为ρ1的流体中浸润饱和，测试各段岩样饱和液后的质量G1，计算各段岩心表观体积Vt，有效孔隙体积Vp和岩样有效孔隙度
S3、利用高速离心和核磁共振技术，对岩心进行换向离心操作，使其含水率达到设定值，建立起各小段岩心样品原始含水饱和度Swi，测试各段岩心原始含水饱和度Swi时的重量G2，利用称重法校正含水饱和度；
S4、对模拟地层水和钻井液滤液的核磁共振T2谱进行测试并判断两者之间的差异；利用核磁成像技术对各小段岩心在原始含水饱和度下的内部流体分布状态做核磁成像测试，并对测试图像进行伪彩处理；测试各小段岩心的原始含水饱和度状态下的气体渗透率K1；
S5、将各小段岩心按照切开前的顺序放入钻井液循环仪的岩心夹持器，进行钻井液循环；循环后，计算各小段岩心的含水饱和度Sqwi，测试含水饱和度Sqwi状态下的各小段岩心重量G3和气体渗透率K2，用测试结果计算对比各小段岩心钻井液侵入后造成的渗透率伤害率Ir，评价钻井液侵入导致的不同深度段的伤害程度；
S6、对钻井液侵入模拟的各段岩心进行T2加权核磁成像测试，测试参数与步骤S4中相同；将钻井液侵入后的各段岩心核磁成像的图像与原始含水饱和度下核磁成像测试图像对比；分析钻井液侵入前后各段岩心内部流体分布变化，明确钻井液液相侵入深度Ll；
S7、岩心烘干处理，对各小段岩心进行扫描电镜分析，以此确定钻井液固相侵入深度Ls，分析岩心内部微观特征差异。
2.根据权利要求1所述的一种致密气储层钻井液伤害评价实验方法，其特征在于，所述步骤S1具体包括如下内容：
在井筒取芯的钻井液侵入方向上钻取柱塞样品，保持岩样端面和柱面平整，端面垂直柱面，柱塞结构完整；柱塞样品直径为2.5cm，柱塞样品长度不小于直径的1.5倍；
对岩样采用甲醇试剂进行除盐处理，处理后进行岩样烘干，烘干温度控制在60～65℃，相对湿度控制在40％～50％；岩样烘干至恒重，烘干时间48h，48h后每2h称量一次，两次称量的差值小于10mg；
将烘干的岩心沿岩心轴向在侧面画一由左向右的箭头指示方向，将岩心均匀切成n等份，每段长度1.4～1.6cm，由左至右分别标号，标号的数值为1、2、…、n；
分别称取各小段岩心的干样重量G0n，测量各小段岩心的长度Ln直径Rn，n的数值为1、2、…、n，各小段岩心干样测试气体渗透率K0n和气体孔隙度n的数值为1、2、…、n，测试围压3.5MPa，驱替压差1MPa，根据实验岩心的渗透率特征，如果渗透率过低，测试压差可相应增加，并保证围压是驱替压力的2倍以上。
3.根据权利要求1所述的一种致密气储层钻井液伤害评价实验方法，其特征在于，所述步骤S2具体包括如下内容：
将烘干后的岩样抽真空20h以上，将模拟地层水放入进行饱和；
采取25MPa高压作为饱和压力，加压时间不低于20h，保证致密砂岩岩样充分饱和模拟地层水；
岩样在饱和液中浸泡至少24h以上，测定饱和液体后各小段岩心质量G1n，n的数值为1、2、…、n；
计算岩样的有效孔隙体积Vp与有效孔隙度
Vp＝(G1–G0)/ρ1；
Vt＝L×A；
A＝π×(D/2)2；
G0-干岩样质量，g；
G1-岩心饱和水后重量，g；
ρ1-测试温度下的饱和液体的密度，g/cm3；
Vρ-岩样有效孔隙体积，cm3；
Vt-岩心表观体积，cm3；
L-样品长度，cm；
D-样品直径，cm；
A-样品端面积，cm2；
-岩样有效孔隙度，％。
4.根据权利要求1所述的一种致密气储层钻井液伤害评价实验方法，其特征在于，所述步骤S3具体包括如下内容：
利用高速离心机对饱和模拟地层水的岩心进行离心操作，依次在100psi、200psi、300psi、400psi、500psi、600psi离心力下进行恒速离心，离心过程中每8min调换岩心方向，离心时间40min；在每个离心速度下完成离心操作后，通过核磁T2谱测试计算含水饱和度，由于含水饱和度建立难以保证与初始含水饱和度完全相同，当测试值接近岩心原始含水饱和度Swi后，不再进行更高一级的离心速度离心；
对核磁建立的原始含水饱和度状态的各段岩心称重G2n，n的数值为1、2、…、n，应用称重法，根据Swi＝{1-[(G2-G0)/ρ1]/Vp}×100％计算含水饱和度，校正核磁共振法建立的原始含水饱和度；
如果储层原始含水饱和度Swi小于30％，还需联合N2气驱替、冰箱冷藏方式建立含水饱和度。
5.根据权利要求1所述的一种致密气储层钻井液伤害评价实验方法，其特征在于，所述步骤S4具体包括如下内容：
分别对10mL模拟地层水和10mL钻井液滤液进行核磁共振T2谱测试，获得模拟地层水的T2时间：T2,D-H2O，和钻井液滤液T2时间：T2,Z-H2O，明确两种流体的T2时间差异；其测试参数包括等待时间TW—2500ms和回波时间TE—100ms；
在明确两种流体的核磁共振T2时间差异后，利用核磁成像测试各小段岩心原始含水饱和度Swi状态的内部流体分布状态，为便于对比钻井液侵入后的流体分布，需要选取钻井液滤液信号较强的测试参数进行测量，选取等待时间TW为500ms、回波时间TE为10ms；对成像获取的灰度图像进行伪彩处理，通过图片中亮度分析流体分布；
测试各段岩心的原始含水饱和度Swi状态下的气体渗透率K1n，n的数值为1、2、…、n，测试介质为N2气，测试围压3.5MPa，驱替压差1MPa；对于离心法建立的岩心含水饱和度，如果与地层原始含水饱和度相差±2％，则判断所测渗透率即为地层状态下岩心渗透率，若相差较大，则通过测量的渗透率通过克里金插值计算岩心原始含水饱和度对应的渗透率值K1n，n的数值为1、2、…、n，如果测试时间超过30min，则须用称重法校正含水饱和度。
6.根据权利要求1所述的一种致密气储层钻井液伤害评价实验方法，其特征在于，所述步骤S5具体包括如下内容：
将各小段岩心按照切开前的顺序排列，依次放入钻井液循环仪的岩心夹持器，岩心标记的箭头指向钻井液侵入的方向，设置围压10MPa；
将井场取回的钻井液盛入釜体，循环压差在压差△P＝3.5MPa条件下循环，也可以根据岩心样品的差异，根据地层压力系数与钻井液密度差异计算相应的循环压差；循环时间根据现场钻井时长选择，为评价短时间和长时间钻井液对气层的伤害，选择需要记录的实验时间；实验过程中记录不同时段对应的钻井液滤失量；
钻井液循环后，利用核磁共振法测试每小段岩心的含水饱和度Sqwin，n的数值为1、2、…、n，并称取各小段岩心钻井液伤害后的重量G3n，n的数值为1、2、…、n，用称重法，根据Sqwi＝{[(G3-G2)/ρ2]/Vp}×100％+Swi校正核磁测试含水饱和度；
对各段岩心进行气体渗透率测试，K2n，n的数值为1、2、…、n；测试介质为氮气，测试围压3.5MPa，驱替压差1MPa；依据渗透率变化率公式Ir＝(K1-K2)/K1×100％，计算对比各小段岩心钻井液侵入后造成的渗透率伤害率Irn，n的数值为1、2、…、n；并依据渗透率伤害率Irn评价钻井液侵入导致的不同深度段的伤害程度，其评价标准如下：钻井液伤害率为Irn≤5％时伤害程度为无伤害，钻井液伤害率为5％<Irn≤30％时伤害程度为弱，钻井液伤害率为30％<Irn≤50％时伤害程度为中等偏弱，钻井液伤害率为50％<Irn≤70％时伤害程度为中等偏强，钻井液伤害率为70％<Irn≤90％时伤害程度为强，钻井液伤害率为90％<Irn时伤害程度为极强。
7.根据权利要求1所述的一种致密气储层钻井液伤害评价实验方法，其特征在于，所述步骤S6具体包括如下内容：
对钻井液侵入模拟的各段岩心进行T2加权核磁成像测试，即测试参数具有以下特点：等待时间小于模拟地层水T2时间，Tw＜T2,D-H2O，回波时间TE尽可能小，这样有利于压制岩心中模拟地层水的信号，使模拟地层水和钻井液滤液信号有效区分；等待时间TW—500ms、回波时间TE—10ms；
将核磁成像的灰度图像进行伪彩处理，彩色区域为水信号，颜色较暗的为模拟地层水，颜色亮的信号为钻井液滤液；
将侵入后的各段岩心核磁成像的图像与原始含水饱和度Swi状态所得到的核磁成像测试图像对比；分析钻井液侵入前后各段岩心内部流体分布变化，明确钻井液液相侵入深度Ll。
8.根据权利要求1所述的一种致密气储层钻井液伤害评价实验方法，其特征在于，所述步骤S7具体包括如下内容：
对实验岩心进行烘干，烘干温度控制在60～65℃，相对湿度控制在40％～50％；岩样烘干至恒重，烘干时间48h，48h后每2h称量一次，两次称量的差值小于10mg；
对岩心切片前，整体按箭头标记方向排列，由左至右，左侧为钻井液侵入方向，将各小段岩心沿轴向中线剖开，对各段岩心剖开面进行连续扫描电镜分析，确定钻井液固相侵入深度Ls，分析岩心内部微观特征差异。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610390509.0A CN106093299B (zh) | 2016-06-02 | 2016-06-02 | 一种致密气储层钻井液伤害评价实验方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610390509.0A CN106093299B (zh) | 2016-06-02 | 2016-06-02 | 一种致密气储层钻井液伤害评价实验方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106093299A CN106093299A (zh) | 2016-11-09 |
| CN106093299B true CN106093299B (zh) | 2019-06-11 |
Family
ID=57447363
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610390509.0A Expired - Fee Related CN106093299B (zh) | 2016-06-02 | 2016-06-02 | 一种致密气储层钻井液伤害评价实验方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106093299B (zh) |
Families Citing this family (40)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10422916B2 (en) * | 2017-08-10 | 2019-09-24 | Saudi Arabian Oil Company | Methods and systems for determining bulk density, porosity, and pore size distribution of subsurface formations |
| CN107631974B (zh) * | 2017-09-07 | 2019-11-12 | 中国海洋石油集团有限公司 | 一种致密砂岩蜡封岩心含水饱和度的获取方法 |
| CN109488290A (zh) * | 2017-09-12 | 2019-03-19 | 中国石油天然气股份有限公司 | 钻井液对储层伤害程度的评价方法及装置 |
| CN107807078B (zh) * | 2017-09-13 | 2019-11-08 | 中国石油天然气股份有限公司 | 地层条件下岩心可动水饱和度在线检测方法 |
| CN108049863A (zh) * | 2018-01-11 | 2018-05-18 | 中国海洋石油集团有限公司 | 一种钻井液化学渗透压差测试装置及测试方法 |
| CN107965310A (zh) * | 2018-01-11 | 2018-04-27 | 中国海洋石油集团有限公司 | 一种钻井液侵入深度测试装置和测试方法 |
| CN108776093A (zh) * | 2018-05-18 | 2018-11-09 | 中国海洋石油集团有限公司 | 一种基于高温高压气藏特征的钻井液损害程度评价方法 |
| CN109061099B (zh) * | 2018-06-11 | 2021-01-05 | 西南石油大学 | 一种非均质致密岩石损伤程度的非破坏性实验评价方法 |
| CN109239311B (zh) * | 2018-10-23 | 2024-03-22 | 中国石油化工股份有限公司 | 堵剂充满度测试方法 |
| CN109254027A (zh) * | 2018-11-01 | 2019-01-22 | 西南石油大学 | 一种基于低场核磁共振的泥饼封堵评价装置及评价方法 |
| CN109932296B (zh) * | 2019-02-26 | 2021-09-28 | 中国石油天然气股份有限公司 | 一种定量表征贾敏效应动态变化的方法 |
| CN109973063B (zh) * | 2019-03-13 | 2020-05-08 | 中国石油大学(北京) | 用于确定碳化水驱油过程中碳化水对储层伤害程度的方法 |
| CN109916799B (zh) * | 2019-03-22 | 2021-04-20 | 西南石油大学 | 测量非常规致密气藏自发渗吸相对渗透率的实验方法 |
| CN109946437B (zh) * | 2019-04-15 | 2021-07-23 | 西南石油大学 | 一种兼顾基块与裂缝系统的裂缝性致密储层工作液损害评价方法 |
| CN109946436B (zh) * | 2019-04-15 | 2021-06-29 | 西南石油大学 | 一种兼顾基块与裂缝的裂缝性致密气层工作液损害评价方法 |
| CN110376113A (zh) * | 2019-07-29 | 2019-10-25 | 中国地质大学（北京） | 一种评价煤储层水锁伤害方法和水锁解除液 |
| CN110376112A (zh) * | 2019-07-29 | 2019-10-25 | 中国地质大学（北京） | 一种评价致密气储层水锁伤害程度的方法 |
| CN110487835B (zh) * | 2019-09-18 | 2022-04-05 | 西南石油大学 | 一种计算致密油气藏储层饱和度指数的新方法 |
| CN111337408B (zh) * | 2020-03-27 | 2022-02-08 | 西南石油大学 | 一种利用低场核磁共振设备测试岩石裂缝孔隙度的方法 |
| CN111827992B (zh) * | 2020-07-29 | 2023-07-25 | 中海石油（中国）有限公司天津分公司 | 一种海上油田生产井储层伤害类型和程度诊断方法 |
| CN111967162A (zh) * | 2020-08-20 | 2020-11-20 | 西南石油大学 | 一种致密砂岩气藏钻井液储层伤害评价方法 |
| CN111852465B (zh) * | 2020-08-27 | 2023-06-20 | 中国石油大学(华东) | 一种水基泥浆钻井获得页岩储层岩心原始含水量评价方法 |
| CN112213345A (zh) * | 2020-09-18 | 2021-01-12 | 中国石油大学(华东) | 一种模拟高温高压条件的核磁共振系统及方法 |
| CN114594033B (zh) * | 2020-12-04 | 2024-06-25 | 中国石油天然气股份有限公司 | 一种评价油气储层水敏、水锁和固相伤害的方法 |
| CN113218834A (zh) * | 2021-03-31 | 2021-08-06 | 长江大学 | 定量描述致密气压裂液与储层渗吸伤害的实验装置及方法 |
| CN113310867B (zh) * | 2021-05-14 | 2022-03-18 | 西南石油大学 | 模拟气藏高温高压地层反渗吸水锁伤害的实验测试方法 |
| CN113358683B (zh) * | 2021-06-11 | 2022-04-12 | 西南石油大学 | 一种研究岩心端面效应的水驱油实验装置及方法 |
| CN113655082B (zh) * | 2021-10-15 | 2021-12-21 | 西南石油大学 | 一种评价致密页岩储层的入井流体的优选方法 |
| CN114088602B (zh) * | 2021-11-19 | 2024-01-26 | 西南石油大学 | 一种基于油层钻屑的储层工作液损害评价方法 |
| CN114414428B (zh) * | 2022-01-13 | 2022-09-02 | 东北石油大学 | 一种评价侵入液对页岩油可动量影响的实验方法 |
| CN114152555B (zh) * | 2022-02-09 | 2022-04-15 | 成都理工大学 | 一种钻井液漏失损害模拟装置及方法 |
| CN114778402A (zh) * | 2022-03-09 | 2022-07-22 | 中国石油天然气股份有限公司 | 一种松散岩心储层伤害评价方法及装置 |
| CN114632368B (zh) * | 2022-04-22 | 2024-04-26 | 中国地质调查局水文地质环境地质调查中心 | 一种有机污染场地多滤层循环井解堵破阻方法 |
| CN115306373A (zh) * | 2022-09-16 | 2022-11-08 | 陕西延长石油（集团）有限责任公司 | 一种适用于不同粘土矿物含量储层伤害的预测方法 |
| CN115717526A (zh) * | 2022-12-05 | 2023-02-28 | 北京斯迪莱铂油气技术有限公司 | 一种用于预测致密气储层水力压裂伤害程度的方法及系统 |
| CN116050623B (zh) * | 2023-02-01 | 2023-05-30 | 西南石油大学 | 一种致密气储层供气能力的计算与评价方法 |
| CN115961950B (zh) * | 2023-02-03 | 2024-06-25 | 西南石油大学 | 基于岩石声学信息的水岩作用下钻井液密度下限预测方法 |
| CN117113884B (zh) * | 2023-10-17 | 2024-01-26 | 西南石油大学 | 一种有水气藏气井井周剩余含水饱和度的确定方法 |
| CN117269224B (zh) * | 2023-11-21 | 2024-02-09 | 天津朔程科技有限公司 | 一种气层水锁伤害的评价方法及系统 |
| CN117591777B (zh) * | 2024-01-18 | 2024-04-30 | 成都英沃信科技有限公司 | 一种水淹气井复产的井周临界含水饱和度确定方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120151998A1 (en) * | 2010-12-21 | 2012-06-21 | Schlumberger Technology Corporation | Wettability and matrix imbibition analysis |
| WO2015112449A1 (en) * | 2014-01-24 | 2015-07-30 | Schlumberger Canada Limited | Workflow for resaturation and multidimensional nmr analysis of unconventional core samples |
| CN104695945A (zh) * | 2015-02-06 | 2015-06-10 | 中国海洋石油总公司 | 一种钻井液滤液侵入储层深度测量方法和系统 |
-
2016
- 2016-06-02 CN CN201610390509.0A patent/CN106093299B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN106093299A (zh) | 2016-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106093299B (zh) | 一种致密气储层钻井液伤害评价实验方法 | |
| CN106501144B (zh) | 一种基于核磁共振双截止值的致密砂岩渗透率计算方法 | |
| CN109030311B (zh) | 基于核磁共振t2谱敏感参数的孔隙结构分类与识别方法 | |
| CA2730532C (en) | Formation evaluation instrument and method | |
| US5086643A (en) | System and method for determining multi-phase relative permeability of a subterranean reservoir | |
| Wan et al. | Pore-scale analysis of gas huff-n-puff enhanced oil recovery and waterflooding process | |
| Abuamarah et al. | A proposed classification for the reservoir quality assessment of hydrocarbon-bearing sandstone and carbonate reservoirs: A correlative study based on different assessment petrophysical procedures | |
| US7388373B2 (en) | Method of determining the permeability of an underground medium from NMR measurements of the permeability of rock fragments from the medium | |
| CN106248712B (zh) | 岩心中微裂缝的造缝方法、微裂缝密度的测量方法和模型建立方法及微裂缝岩心的制备方法 | |
| CN105604546B (zh) | 双重介质碳酸盐岩储层的定量分类方法 | |
| CN111007230A (zh) | 定量评价陆相湖盆低孔隙度致密油储层含油量的方法 | |
| CN109061099B (zh) | 一种非均质致密岩石损伤程度的非破坏性实验评价方法 | |
| CN110208487A (zh) | 一种基于ct扫描的页岩水化损伤测试方法 | |
| Lenormand et al. | Advances in measuring porosity and permeability from drill cuttings | |
| CN113916745B (zh) | 一种无损测量水驱气藏微观孔隙结构变化规律的实验方法 | |
| Ma et al. | Nuclear magnetic resonance experiment on the influence of confining pressure on spontaneous imbibition of water in coal | |
| CN112487620B (zh) | 一种页岩油可动资源量的评价方法 | |
| Ren et al. | Determination of microscopic waterflooding characteristics and influence factors in ultra-low permeability sandstone reservoir | |
| Ubani et al. | ADVANCES IN CORING AND CORE ANALYSIS FOR RESERVOIR FORMATION EVALUATION. | |
| Lynch et al. | Development of a laboratory technique for obtaining soil water retention curves under external loading in conjunction with high-capacity tensiometers | |
| CN110231268B (zh) | 一种致密储层岩石油水赋存孔径分布的分析方法 | |
| US11852576B2 (en) | Method for determining the pore size distribution in a reservoir | |
| Yiren et al. | Simulation of mud invasion and analysis of resistivity profile in sandstone formation module | |
| Fleury | FRIM: a fast resistivity index measurement method | |
| Odiachi et al. | Determination of pore fluid salinity in tight rocks without fluid extraction |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190611 Termination date: 20200602 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |