本发明专利技术涉及求取渗透率的岩石激发极化谱分析方法。解决了现有方法测井成本高、信号弱及误差大精确度低的问题。其特征在于：在井下采集激发极化衰减谱，经过多指数衰减拟合得到激发极化驰豫时间谱，然后求出驰豫时间的几何平均值和算术平均值，结合孔隙度测井资料和实验常数，计算得到地层的渗透率曲线。具有测井成本相对较低、信号幅度大、信噪比高、衰减时间足够长及测量渗透率精度高的特点。
【技术实现步骤摘要】
本专利技术涉及一种激发极化谱测井方法及其解释方法，特别是。
技术介绍
储层孔隙度的测量和评价有比较成熟的测井技术和解释方法，但在确定渗透率方面，由于决定渗透率大小的主导因素是储层孔隙结构，而孔隙结构的复杂性及其严重的非均质性，使得渗透率测井技术和解释方法一直没有很大进展，目前主要采用实验室测量和现场核磁共振测井，实验室测量难以再现地层条件，在取芯过程会使岩心产生不必要的裂纹，这些都对渗透率的测量产生极大的误差。而核磁共振测井方法存在的缺陷是测井成本非常高，信号微弱，探测深度浅，对仪器的信噪比要求也非常高，同时该技术测量了不连通的孔隙，这些孔隙对渗透率没有贡献。传统的激发极化测井主要是利用极化率参数，建立极化率与地层水矿化度以及泥质含量之间的关系，但是极化率参数的影响因素非常复杂，比如围岩、井眼，供电时间、泥质分布等，尤其重要的是极化率参数与测量方式以及测量时间密切相关，对仪器的要求非常高。
技术实现思路
为了克服现有方法测井成本高、信号弱及误差大、精确度低的不足，本专利技术提供一种，该方法首先由恒流源通过控制单元给供电电极供电，使岩心产生极化场，然后断电，极化场随时间逐渐衰减，测量衰减弱，对衰减谱进行多指数衰减反演，得到弛豫时间谱，然后求出驰豫时间的几何平均值和算术平均值，结合孔隙度测井资料和实验常数，计算得到地层的渗透率曲线。本专利技术的技术方案是一种，a、在井下采集激发极化衰减谱，经过多指数衰减拟合 V(t)V(0)=Σi=1nfiexp(-t/Ti)---(1)]]>fi代表时间常数为Ti的孔体积百分比，且fi≥0；b、将方程(1)改写成如下形式y＝M·f (2)式中y＝(y1，y2，...，ym)T是归一化的激发极化衰减电位，即yi≡V(ti)V(0),]]>t(i＝1，2，...，m)为衰减时间，V(ti)为衰减时间为ti的衰减电位，V(0)为断电瞬间(即t＝0)时的衰减电位，M＝m×n，，(m＞n)，f＝(f1，f2，...，fn)T衰减项的系数，Tj(j＝1，2，...，n)为预先指定的驰豫时间时间分布，其典型取法为在某个合适的区间(Tmin，Tmax)按照对数规律均匀地选取n个点，其中Tmin为选择的最小T，Tmax为选择的最大的T，c、求解方程(2)得到解f，且f对T作图可以得到激发极化弛豫时间谱。d、时间常数Ti和孔径Ri之间有如下的关系Ri=D×Ti---(3)]]>其中Ri为第i个孔的半径，D为溶质的扩散系数；每一个时间常数代表一种孔径，而前面的系数代表这种孔径的体积百分比；e、渗透率模型为K＝c×(φm×T)b(4)式中c、m和b为常数，K为渗透率，φ为孔隙度；f、求出驰豫时间的几何平均值Tg和算术平均值TaTa=Σi=1nfi×TiΣi=1nfi---(5)]]>Tg=(Πi=1nTifi)1Σi=1nfi---(6)]]> g、取m＝4，且用几何平均值和算术平均值代替渗透率模型中的弛豫时间平均值，则式(4)可以写成如下形式K＝c×(φ4×Tg)b(7)K＝x×(φ4×Ta)b(8)h、渗透率模型(7)和(8)中的系数c和b是通过每个区块的岩心的孔隙度和渗透率参数进行拟合而得到的；i、结合相应区块现场测量的孔隙度和弛豫时间常数平均值，利用渗透率模型(7)和(8)得到相应的渗透率。本专利技术具有的有益效果是岩石的激发极化电位衰减谱中包含许多有用的地质信息，包括孔隙结构、渗透率、地层水矿化度以及泥质含量，衰减谱的形状与地层的孔隙结构相关，从而可以反映渗透率参数。因此可以利用激发极化衰减谱求取渗透率，这种技术信号幅度大，信噪比高，探测深度深，只测量含水的孔隙，不测量含油孔隙，因此按上述模型求取的渗透率误差小，精确度高。利用激发极化衰减谱，经过多指数衰减拟和得到驰豫时间谱，建立了渗透率与孔隙度、驰豫时间几何平均值之间的定量关系。我们可以利用测量得到的一段衰减谱，经过拟合也能得到较好的结果，影响因素少，由于测量的是衰减谱，因此与传统的激发极化测井相比，对仪器的要求要低。基于激发极化谱分析技术的激发极化谱测井属于电测井范畴，具有测井成本相对较低、信号幅度大、信噪比高、衰减时间足够长及测量渗透率精度高的特点，因此激发极化谱分析技术具有非常好的应用前景。附图说明图1是本专利技术典型岩心的激发极化衰减谱及其相应的弛豫时间谱；图2岩心渗透率与孔隙度之间的关系图；图3岩心渗透率与孔隙度、驰豫时间几何平均值之间的关系图；图4岩心渗透率与孔隙度、驰豫时间算术平均值之间的关系图。图中a-激发极化衰减谱，b-反演得到的驰豫时间谱，c-区块1，d-区块2，e-区块3，f-区块1，g-区块2，h-区块3，i-区块1，j-区块2，k-区块3。具体实施例方式下面结合附图将对本专利技术作进一步说明该方法采用激发极化理论建立在双电层基础之上，即处于深液中的岩石颗粒，尤其是粘土颗粒因为带有多余的负离子而吸引溶液中的阳离子，在其周围形成双电层，当在岩石的两端施加一恒定外电场后，则会使得双电层发生形变或局部浓度的变化，即产生极化电场，当外电场去掉后，由于扩散电位，即二次场衰减电位，这就是激发极化电位。激发极化测井技术通过探测孔隙内流体的电驰豫特征，能够提供储层渗透率、地层水矿化度、孔径大小和分布等信息。岩石在完全饱和水的情况下，孔隙内水中离子的电驰豫特性受到岩石孔隙结构的影响，表现出不同的驰豫特性，激发极化测井测量的信号是由不同大小孔隙内水的离子极化信号的叠加，经过复杂的数学拟合得到激发极化弛豫时间谱，因此激发极化弛豫时间谱的分布反映了孔隙大小的分布，大孔隙内的组分对应长的激发极化弛豫时间值，小孔隙组分对应短的激发极化弛豫时间值。既然激发极化弛豫时间分布谱和压汞曲线都反映了岩石的孔隙结构特征，那么两者之间必然存在密切的关系，可以应用激发极化测井资料对孔隙结构进行评价，确定毛管压力曲线以及渗透率值。理论分析表明，单个孔隙中的激发极化二次场电位的衰减满足单指数衰减规律，但是岩石内部是由一系列大小不等的孔隙群体组成的，所以在岩石激发极化二次场电位V(t)是一系列单个孔隙极化电位的叠加，二次场V(t)按照t＝0时的二次场V(0)归一化，然后进行多指数拟合V(t)V(0)=Σi=1Nfiexp(-t/Ti)---(1)]]>fi代表时间常数为Ti的孔体积百分比，其中fi≥0。该方程可以写成如下形式 y＝M·f(2)式中y＝(y1，y2，...，ym)T是归一化的激发极化衰减电位，即yi≡V(ti)V(0),]]>ti(i＝1，2，...，m)为衰减时间，V(ti)为衰减时间为ti的衰减电位，V(0)为断电瞬间(即t＝0)时的衰减电位，M＝m×n，，(m＞n)，f＝(1，f2，...，fn)T衰减项的系数，Tj(j＝1，2，...，n)为预先指定的驰豫时间时间分布；T的典型取法为在某个合适的区间(Tmin，Tmax)按照对数规律均匀地选取n个点，其中Tmin为选择的最小T，Tmax为选择的最大的T，求解方本文档来自技高网...
【技术保护点】
一种求取渗透率的岩石激发极化谱分析方法，其特征在于： ａ、在井下采集激发极化衰减谱，经过多指数衰减拟合： Ｖ（ｔ）／Ｖ（０）＝＊ｆ↓［ｉ］ｅｘｐ（－ｔ／Ｔ↓［ｉ］） （１） ｆ↓［ｉ］代表时间常数为Ｔ↓［ｉ］的孔体积百分比，其中ｆ↓［ｉ］≥０； ｂ、多指数衰减方程可以写成如下形式： ｙ＝Ｍ.ｆ （２） 式中ｙ＝（ｙ↓［１］，ｙ↓［２］，…，ｙ↓［ｍ］）↑［Ｔ］是归一化的激发极化衰减电位，即ｙ↓［ｉ］≡Ｖ（ｔ↓［ｉ］）／Ｖ（０），ｔ↓［ｉ］（ｉ＝１，２，…，ｍ）为衰减时间，Ｖ（ｔ↓［ｉ］）为衰减时间为ｔ↓［ｉ］的衰减电位，Ｖ（０）为断电瞬间（即ｔ＝０）时的衰减电位，Ｍ＝［ｅｘｐ（－ｔ↓［ｉ］／Ｔ↓［ｊ］）］↓［ｍ×ｎ］，（ｍ＞ｎ），ｆ＝（ｆ↓［１］，ｆ↓［２］，…，ｆ↓［ｎ］）↑［Ｔ］衰减项的系数，Ｔ↓［ｊ］（ｊ＝１，２，…，ｎ）为预先指定的驰豫时间时间分布； ｃ、Ｔ的典型取法为在某个合适的区间（Ｔ↓［ｍｉｎ］，Ｔ↓［ｍａｘ］）按照对数规律均匀地选取ｎ个点，其中Ｔ↓［ｍｉｎ］为选择的最小Ｔ，Ｔ↓［ｍａｘ］为选择的最大的Ｔ； ｄ、求解方程（２）得到解ｆ； ｅ、ｆ对Ｔ作图可以得到激发极化弛豫时间谱。...
【技术特征摘要】
【专利技术属性】
技术研发人员：童茂松，王伟男，姜亦忠，范清华，
申请(专利权)人：大庆石油管理局，
类型：发明
国别省市：23[中国|黑龙江]
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。