测井资料处理与综合解释实验

1、 测井资料处理与解释测井资料处理与综合解释实验专业名称：勘查技术与工程 学生姓名：张奎涛学生学号：201326020109指导老师：岳崇旺、王飞完成日期：2017-1-2目 录实验一 定性划分储集层并定量解释.1实验二 利用综合方法估计地层泥质含量.7实验三 含泥质复杂岩性地层综合测井处理.18实验一 定性划分储集层并定量解释一、 实验目的：通过对测井曲线特征的分析和认识，掌握定性划分砂泥岩剖面储集层的基本方法，并应用阿尔奇公式，进行储层参数的计算，巩固已经学过的钻井地球物理课程的主要内容与应用。二、 实验要求正确划分出储集层和非储集层，对砂泥岩剖面能区分开较明显的油水层。进行测井曲线读数，简

2、单地计算出孔隙度、饱和度等参数。三、 实验场地、用具与设备测井实验室或一般的教室，长直尺、铅笔、像皮和计算器四、 实验内容：1 测井曲线图的认识；图1是某井的综合测井曲线图。图中共有5道，第一道主要为反映岩性的测井曲线道，包括：自然电位测井曲线曲线符号为SP、记录单位mv；自然伽马测井曲线曲线符号为GR、记录单位API；井径测井曲线曲线符号为CAL，记录单位in或cm；岩性密度测井曲线（光电吸收界面指数）曲线符号为PE；第二道是深度道；通常的深度比例尺为1：200 或1：500第三道是反映含油性的测井曲线道，包括深中浅三条电阻率测井曲线，分别是：深侧向测井曲线曲线符号为LLD、记录单位m；浅侧

3、向测井曲线曲线符号为LLS、记录单位m；微球形聚焦测井曲线曲线符号为MSFL、记录单位m；电阻率测井曲线通常为对数刻度。第四道为反映孔隙度的测井曲线道，包括：密度测井曲线曲线符号为DEN或RHOB，记录单位g/cm3；中子测井曲线曲线符号为CNL或PHIN，记录单位%，有时为v/v。声波测井曲线曲线符号为AC或DT，记录单位us/ft，有时为us/m。中子和密度测井曲线的刻度的特点是保证在含水砂岩层上两条曲线重迭，在含气层上，密度孔隙度大于中子孔隙度，在泥岩层上，中子孔隙度大于密度孔隙度；第五道是反映粘土矿物类型的测井曲线道，包括自然伽马能谱测井中的三条曲线：放射性钍测井曲线曲线符号为Th或T

4、HOR，记录单位是ppm；放射性铀测井曲线曲线符号为U或URAN，记录单位ppm；放射性钾测井曲线曲线符号为K或POTA，记录单位%，有时为v/v。2 测井曲线特征（1）砂泥岩剖面的测井曲线特征砂泥岩剖面储集层（砂岩）的典型特征是，一般自然电位有明显的异常，异常的方向和幅度取决于泥浆滤液电阻率（Rmf）和地层水的电阻率（Rw），或者说与Rmf与Rw的比值有关,如果Rmf> Rw，则为负异常，否则为正异常。如果砂层中不含放射性矿物，自然伽玛曲线亦显示低值。微电极曲线一般在砂岩层幅值高，并出现正幅差。而泥岩的幅度和幅差均较低，当井眼条件不好时，可能会出现曲线跳动现象。砂岩中含灰质较多的夹层，

5、因为致密电阻率异常高，幅度差很小或没有。一般幅度差的大小标明了储集层渗透性的好坏。普通电阻率测井曲线在泥岩处显示为低值。砂岩处显示为高值，含油砂岩幅值就更高，如有两条探测深度不同的Ra 曲线，幅值的差别显示着低侵、高侵。通常在油层上为低侵，水层上为高侵。井径在泥岩层扩大，砂岩层缩小（略小于钻头直径）。具体特征总结见表1表1 砂泥岩剖面测井曲线特征 地层测井曲线储集层砂岩非储集层泥岩自然电位负异常（Rw<Rmf）正异常（Rw>Rmf）泥岩基线自然伽马低高井径缩径扩径深中浅电阻率高阻低阻声波<300us/m>300us/m钍低高铀低高钾低高（2）碳酸盐岩剖面的测井曲线特征碳

6、酸盐岩剖面的测井解释任务，就是从致密的围岩中找出孔隙性、裂缝性的储集层，并判断其含油性。碳酸盐岩剖面电阻率一般较高，自然电位效果不好。为了区分岩性和划分储层，一般使用自然伽马测井曲线。 储集层相对于致密的围岩具有低阻、低自然伽马以及孔隙度测井反映孔隙度较大的特点。3 划分储集层的基本方法与原则基本要求：凡一切可能含油气的地层都要划出来，要适当地划分明显的水层。具体要求为：（1） 估计为油层、气层、油水同层和含油水层的储集层都必须分层解释。（2） 厚度半米以上的电性（测井曲线）可疑层（即指从测井曲线上看有油气的地层）或录井显示为微含油级别以上的储集层必须分出。（3） 选择出作为确定地层水电阻率R

7、w 的标准水层（厚度大、岩性纯、不含油）要划分出来。（4） 录井、气测有大段油气显示而测井曲线显示不好的储集层，应选取一定层位，尤其是该组储层的顶部层位，进行分层。（5） 当有多套油水系统，油层组包括若干水层时，只解释最靠近油层的水层。（6） 对于新区探井，应做细致工作，对各个储层均应酌情选层解释，以使不漏掉可能有油气的地层。4 正确划分出储集层的方法（1）砂泥岩剖面通常是自然电位（SP）曲线的异常确定渗透层的位置，用微电极曲线确定分层界面，分层前，应将井场收集的井壁取芯、气测显示等有关油气显示的资料标注在综合测井曲线图上，并根据邻井的测井和试油等资料对本井的油水关系作出初步估计。分层时应注意

8、：l 确定分层的界面深度时，应左右环顾，照顾到分层线对每条测井曲线的合理性。l 分层的深度误差不应大于0.1m。l 渗透层中，凡是0.5m以上的非渗透性夹层（泥岩或致密层），应将夹层上下的渗透层分两层解释。l 岩性渐变层顶界（顶部渐变层）或底界（底部渐变层）分层深度应在岩性渐变结束处。l 一个厚度较大的渗透层，如有两个以上解释结论，应按解释结论分层。l 在同一解释井段，如果油气层与水层岩性、地层结构和孔隙度基本相同，则油气层是纯水层的电阻率的3-5倍。纯水层的自然电位异常最大，油气层异常明显偏小，油水同层介于油、水层之间。并且厚度较大的油水同层，自上而下电阻率有明显减小的趋势。（2）碳酸盐岩剖

9、面碳酸岩盐剖面划分渗透性地层的5 测井曲线读数分层以后，要从有关的主要测井曲线将代表该储层的测井曲线读数，以便计算孔隙度、饱和度等地质参数，在厚度较大的储集层中按测井曲线变化确定几个取值区，对每个取值区对应读数计算，几种主要测井曲线取值区的最小厚度如下：各种孔隙度测井0.6m。侧向测井0.6m感应测井，低阻0.6m，高阻层1.5m。每种测井曲线分层和取值要符合其方法特点，例如声波测井扣除致密夹层，选用与渗透层相对应部分的平均值。电阻率测井曲线则扣除致密夹层，选用与渗透层相对应部分的极大值的平均值。另外注意孔隙度与电阻率测井曲线对应取值的原则。因为要用两者结合计算地层的含水饱和度，两者当然应该是

10、对应深度上同一地层或同一取值区的读数。岩层含油性的定性判断，主要依据井曲线的测井曲线特征，而电性特征是岩石物性、岩性和含油性的综合反映。因此在判断地层的含油性时，一般应将测量井段首先按照地层水矿化度的不同分为不同的解释井段，然后才有可能对每一个解释井段在充分考虑其岩性特点的前提下进行含油性解释。由于地下地层复杂性，仪器的局限性，上述原则是一般性的。要做到正确地解释，一方面应多收集资料，认真分析曲线，另一方面还要了解区域性特点和规律，要积累经验。6 计算出孔隙度、饱和度等参数。读数以后，还要做一些定量计算，常用的公式：孔隙度：含水饱和度：上式中为当前层的声波时差，为地层水的声波时差,189us/

11、ft(623us/m)，为固体骨架的声波时差，对于砂岩骨架，主要矿物为石英，其声波时差为55.5us/ft。a是常数，对于砂岩地层通常取1.0，为当前层的电阻率，m为胶结指数。五、 实验结果处理结果如图2所示 具体分层如下表所示：储层序号顶部深度（m）底部深度（m）厚度（m）测井曲线读数孔隙度含水饱和度SPGRDEN11122.4114219.6-6.49100.012.31413.84586.4421154.2116216.8-8.023103.5732.37311.61586.33231218.0512267.95-10.504115.0922.28913.94883.91441236.2

12、1252.616.4-14.695104.6862.33013.22978.80751272.41286.1513.75-17.636106.9982.34813.90163.182由自然电位的划分原则可知，在图2的自然电位的曲线上，有五个明显的偏离泥基线的负异常，故可以划分出五个储集层。然后就可以求出泥质含量，进而求出泥质校正后的孔隙度，然后根据阿尔奇公式，就可求出含水饱和度、渗透率，进而划分油气水层。 图1 某井的综合测井曲线图23451图2 SN183井测井解释综合图实验二利用综合方法估计地层泥质含量一、 实验目的：通过实际计算，巩固掌握利用多种测井资料确定泥质含量的方法。二、 实验要求

13、自编程序，在计算机上运算出地层泥质含量。三、 实验场地、用具与设备计算中心，尺子、像皮和计算机；四、 实验内容： 什么是泥质含量：泥质是指颗粒直径小于0.01mm的碎屑物质，泥质含量，也叫做泥质体积，是指泥质的体积占岩石总体积的比：确定Vsh的重要性泥质含量的确定，在泥质砂岩储集层的定量解释中具有重要意义。多年来人们提出许多计算泥质含量的理论和方法。目前求取泥质含量的方法大致可分为两类，一类是用每种测点各求出一个泥质含量，然后求出最佳值。当岩石含有泥质时，各种测井曲线均或多或少地受到泥质的影响，其影响的程度受Vsh的决定，评价岩石的特性时，只有已知Vsh，才知道由于泥质带来的影响，从而将泥质的

14、影响校正掉。一般而言，用自然伽马或自然伽马能谱或自然电位来求取泥质含量效果最好，但自然伽马要求储层中除了泥质外，其他物质不含放射性矿物。自然电位要求地层水电阻率保持不变，且储层中的泥质与相邻泥岩的的成分相同。用其他方法计算泥质含量则要求更为苛刻的条件：如电阻率方法要求储层的孔隙度和含水饱和度均要很小。中子和声波方法则要求孔隙度很小。3. 确定Vsh的方法：（） 自然伽玛法式中，分别是砂岩和泥岩层的自然伽马值，GCUR是与地层有关的经验系数，新地层（第三系地层）GCUR=3.7，老地层GCUR=2.0.（） 自然电位法式中，是当前层的自然电位读数，和分别是纯地层和泥质地层的自然电位读数（）电阻率

15、（b=1.5）（）中子法式中，是当前层的视中子孔隙度读数，是泥岩层的视中子孔隙度读数。（） 交会图法以中子密度测井交会图为例，通过对图2所示的石英点（Q）、水点（W）和泥岩点（SH）构成的三角形进行分解，依据资料点所落入三角形中的位置，可以推测出来泥质含量。或者利用下式进行计算（依据点到直线的距离计算方法）：式中，SHWQ0是石英点（Q）和水点（W）连线的直线方程。依据任意两点的直线，用石英点()和水点()两个点的参数可以推出：，当然，也可以用中子声波、声波密度交会图的类似方法求Vsh 。如果解释层段上没有纯泥岩层时，上述交会图法所定出的泥岩点位置并不代表实际的泥岩参数，导致用交会图方法估计的

16、泥质含量比实际值偏大。上述几种方法计算的泥质含量往往都有一定的条件，当条件满足时， 泥质含量的近似结果，当条件不满足时，计算的泥质含量均可能偏高，所以在实际处理时，最后选取其中的最小值作为接近实际的泥质含量。五、 实验结果1.已知条件：GRmin=50 GRmax=110 SPmin=0 SPmax=100 b=1.5 Rsh=5 2.测井数据depthspmsflllslldrhobnphigrdtcal2680.9888.504.809.0910.842.5517.00103.1773.1311.992682.0983.402.026.588.212.0720.00108.2577.751

17、4.192683.0170.001.585.717.171.5924.00107.8181.8115.592684.0280.001.925.967.481.9623.00108.7581.2114.792685.0484.302.986.818.282.3020.00105.0878.3812.032686.0578.509.607.158.222.5919.00112.1977.1310.592687.0736.0014.1810.9311.732.6417.0098.4472.409.822688.0926.5019.3818.1319.142.5817.0085.5470.739.49

18、2689.0025.2068.9450.0651.542.5011.0053.4865.949.262690.0227.0055.8145.6645.832.5210.0064.0665.319.163.计算结果spmsflllslldrhob-nphinphigr泥质含量0.885001.027590.671330.596980.542990.566670.805350.542990.834001.829840.832710.718490.463950.666670.947160.463950.700002.155470.915280.786381.433950.800000.934210.

19、700000.800001.892830.889510.764500.609270.766670.962040.609270.843001.412000.813860.714430.071650.666670.856730.071650.785000.647340.787850.717900.710030.633331.069200.633330.360000.499110.593700.566390.752580.566670.687470.360000.265000.405270.423690.408650.612850.566670.424370.265000.252000.173910

20、.215270.211130.204880.366670.027910.027910.270000.200220.228890.228320.214510.333331.279430.20022附录：程序源代码! 功能:利用测井数据求泥质含量! 时间:2016/12/22!-! 参数说明:! cmdfile:参数文件! file_data:输入测井数据文件! file_contant:常数参量文件! file_sp:输出自然电位估计的泥质含量文件! file_msfl:输出微球型聚焦估计的泥质含量文件! file_lls:输出浅侧向估计的泥质含量文件! file_lld:输出深侧向估计的泥质含

21、量文件! file_nphi:输出补偿中子估计的泥质含量文件! file_gr:输出自然伽马估计的泥质含量文件! file_n_d:输出中子-密度交会图估计的泥质含量文件! n:测井数据的个数! grmin,grmax:自然伽马的最小值与最大值! spmin,spmax:自然电位的最小值与最大值! b,r_sh:已知常数! fai_nsh:已知常数! p_sh,p_ma,p_f:已知常数! fai_nf,fai_nma:已知常数! sp(:):自然电位! msfl(:):微球型聚焦! lls(:):浅侧向! lld(:):深侧向! nphi(:):补偿中子! gr(:):自然伽马! rhob(

22、:):岩性密度!-program cjsyimplicit nonecharacter*80 cmdfilecharacter*80 file_data,file_contantcharacter*80 file_sp,file_msfl,file_lls,file_lld,file_nphi,file_gr,file_n_dinteger nreal grmin,grmax,spmin,spmax,b,r_sh,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nmareal,allocatable:sp(:),msfl(:),lls(:),lld(:),nphi(:),gr

23、(:),rhob(:) cmdfile='cmd.txt'call read_cmd(cmdfile,file_data,file_contant,file_sp,file_msfl,file_lls,file_lld,file_nphi,file_gr,file_n_d)call read_number(file_data,n)allocate(sp(1:n),msfl(1:n),lls(1:n),lld(1:n),nphi(1:n),gr(1:n),rhob(1:n)call read_data(file_data,n,sp,msfl,lls,lld,nphi,gr,rho

24、b)call read_contant(file_contant,grmin,grmax,spmin,spmax,b,r_sh,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nma)call get_sp(file_sp,n,sp,spmin,spmax)call get_r(file_msfl,n,msfl,b,r_sh)call get_r(file_lls,n,lls,b,r_sh)call get_r(file_lld,n,lld,b,r_sh)call get_nphi(file_nphi,n,nphi,fai_nsh)call get_gr(file_gr,n,

25、gr,grmin,grmax)call get_n_d(file_n_d,n,nphi,rhob,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nma)deallocate(sp,msfl,lls,lld,nphi,gr,rhob)end!-! ! 功能：读取参数文件! 输入参数说明：! cmdfile:参数文件! 输出参数说明：! file_data:输入测井数据文件! file_contant:常数参量文件! file_sp:输出自然电位估计的泥质含量文件! file_msfl:输出微球型聚焦估计的泥质含量文件! file_lls:输出浅侧向估计的泥质含量文件! fil

26、e_lld:输出深侧向估计的泥质含量文件! file_nphi:输出补偿中子估计的泥质含量文件! file_gr:输出自然伽马估计的泥质含量文件! file_n_d:输出中子-密度交会图估计的泥质含量文件!-subroutine read_cmd(cmdfile,file_data,file_contant,file_sp,file_msfl,file_lls,file_lld,file_nphi,file_gr,file_n_d)implicit nonecharacter*(*) cmdfilecharacter*(*) file_data,file_contantcharacter*(*

27、) file_sp,file_msfl,file_lls,file_lld,file_nphi,file_gr,file_n_dopen(10,file=cmdfile,status='old')read(10,*) file_dataread(10,*) file_contantread(10,*) file_spread(10,*) file_msflread(10,*) file_llsread(10,*) file_lldread(10,*) file_nphiread(10,*) file_grread(10,*) file_n_dclose(10)endsubrou

28、tine read_cmd!-! ! 功能：读取测井数据的个数! 输入参数说明：! filename:输入测井数据文件! 输出参数说明：! number:测井数据的个数!-subroutine read_number(filename,number)implicit nonecharacter*(*) filenameinteger numberopen(11,file=filename,status='old')number=0 do while(.not.eof(11) read (11,*,end=100,ERR=100) number=number+1100 end d

29、oclose(11)number=number-1endsubroutine read_number!-! ! 功能：读取测井数据! 输入参数说明：! file_data:输入测井数据文件! n:测井数据的个数! 输出参数说明：! sp(:):自然电位! msfl(:):微球型聚焦! lls(:):浅侧向! lld(:):深侧向! nphi(:):补偿中子! gr(:):自然伽马! rhob(:):岩性密度!-subroutine read_data(file_data,n,sp,msfl,lls,lld,nphi,gr,rhob)implicit nonecharacter*(*) file

30、_datainteger nreal sp(1:n),msfl(1:n),lls(1:n),lld(1:n),nphi(1:n),gr(1:n),rhob(1:n)integer ireal hopen(10,file=file_data,status='old')read(10,*)do i=1,n,1 read(10,*) h,sp(i),msfl(i),lls(i),lld(i),rhob(i),nphi(i),gr(i)enddoclose(10)endsubroutine read_data!-! ! 功能：读取常数参量! 输入参数说明：! file_contant:

31、常数参量文件! 输出参数说明：! grmin,grmax:自然伽马的最小值与最大值! spmin,spmax:自然电位的最小值与最大值! b,r_sh:已知常数! fai_nsh:已知常数! p_sh,p_ma,p_f:已知常数! fai_nf,fai_nma:已知常数!-subroutine read_contant(file_contant,grmin,grmax,spmin,spmax,b,r_sh,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nma)implicit nonecharacter*(*) file_contantreal grmin,grmax,sp

32、min,spmax,b,r_sh,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nmaopen(10,file=file_contant,status='old')read(10,*) grmin,grmaxread(10,*) spmin,spmaxread(10,*) bread(10,*) r_shread(10,*) fai_nshread(10,*) p_sh,p_ma,p_fread(10,*) fai_nf,fai_nmaclose(10)end subroutine read_contant!-! ! 功能：用自然电位的估计泥质含量! 输入参

33、数说明：! n:测井数据的个数! sp(:):自然电位! spmin,spmax:自然电位的最小值与最大值! 输出参数说明：! filename:输出自然电位估计的泥质含量文件!-subroutine get_sp(filename,n,sp,spmin,spmax)implicit nonecharacter*(*) filenameinteger nreal sp(1:n)real spmin,spmaxinteger ireal v_shopen(10,file=filename,status='unknown')do i=1,n,1 v_sh=(sp(i)-spmin)

34、/(spmax-spmin) write(10,*) v_shenddoclose(10)end subroutine get_sp!-! ! 功能：用电阻率的估计泥质含量! 输入参数说明：! n:测井数据的个数! r(:):电阻率! b,r_sh:已知常数! 输出参数说明：! filename:输出电阻率估计的泥质含量文件!-subroutine get_r(filename,n,r,b,r_sh)implicit nonecharacter*(*) filenameinteger nreal r(1:n)real b,r_shinteger ireal v_shopen(10,file=f

35、ilename,status='unknown')do i=1,n,1 v_sh=(r_sh/r(i)*(1.0/b) write(10,*) v_shenddoclose(10)end subroutine get_r!-! ! 功能：用中子的估计泥质含量! 输入参数说明：! n:测井数据的个数! nphi(:):补偿中子! fai_nsh:已知常数! 输出参数说明：! filename:输出中子估计的泥质含量文件!-subroutine get_nphi(filename,n,nphi,fai_nsh)implicit nonecharacter*(*) filenamei

36、nteger nreal nphi(1:n)real fai_nshinteger ireal v_shopen(10,file=filename,status='unknown')do i=1,n,1 v_sh=nphi(i)/fai_nsh write(10,*) v_shenddoclose(10)end subroutine get_nphi!-! ! 功能：用中子的估计泥质含量! 输入参数说明：! n:测井数据的个数! gr(:):自然伽马! grmin,grmax:自然伽马的最小值与最大值! 输出参数说明：! filename:输出自然伽马估计的泥质含量文件!-su

37、broutine get_gr(filename,n,gr,grmin,grmax)implicit nonecharacter*(*) filenameinteger nreal gr(1:n)real grmin,grmaxinteger ireal v_sh,gcurgcur=2.0open(10,file=filename,status='unknown')do i=1,n,1 v_sh=(gr(i)-grmin)/(grmax-grmin) v_sh=(2*(gcur*v_sh)-1)/(2*gcur)-1) write(10,*) v_shenddoclose(10

38、)end subroutine get_gr!-! ! 功能：用中子-密度交会图的估计泥质含量! 输入参数说明：! n:测井数据的个数! nphi(:):补偿中子! rhob(:):岩性密度! fai_nsh:已知常数! p_sh,p_ma,p_f:已知常数! fai_nf,fai_nma:已知常数! 输出参数说明：! filename:输出中子-密度交会图估计的泥质含量文件!-subroutine get_n_d(filename,n,nphi,rhob,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nma)implicit nonecharacter*(*) filen

39、ameinteger nreal nphi(1:n),rhob(1:n)real fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nmainteger ireal v_sh,a,b,ca=p_ma-p_fb=fai_nf-fai_nmac=fai_nma*p_f-fai_nf*p_maopen(10,file=filename,status='unknown')do i=1,n,1 v_sh=(abs(a*nphi(i)+b*rhob(i)+c)/(abs(a*fai_nsh+b*p_sh+c) write(10,*) v_shenddoclose(10)en

40、d subroutine get_n_d实验三 含泥质复杂岩性地层综合测井处理一、实验目的：综合所学习的测井解释与数据处理知识，掌握一般含泥质复杂岩性地层测井数据处理的步骤和方法。二、实验要求用C语言编程，调试通过，并能计算出正确的结果。三、实验场地、用具与设备计算中心，计算机；四、实验内容：（一） 原理：测井数据处理与解释的主要任务就是依据“四性”关系，用测井所获得的地层的各种“物理参数”，定量地计算出地下地层的岩性（）、储集性(孔隙度和渗透率)以及含油性（饱和度Sw）。当储层含有泥质时，还要求取泥质含量，对测井数据进行泥质校正。对于含有泥质骨架由两种矿物构成,且矿物含量随深度变化的复杂岩性

41、地层，其体积模型如图3所示。这种泥质储层测井处理的过程见图4所示。其基本思路就是用几种方法分别计算出独立的泥质含量Vsh ，然后取最极小值作为该深度的泥质含量。对测井数据进行泥质含量校正，两种孔隙度组合求出孔隙度和两种矿物的含量，在恢复为含泥质地层，计算出实际地层的孔隙度和矿物含量。在计算过程中所设计的参量和模型的变化见图3所示。VshV1V2恢复为泥质岩石求Sw K，VshNwWV1V2 VshV1V2Vma1Vma2纯岩 石 图3 体积模型与计算过程中的变化 实际上求上述四种体积含量的方法有两种：（） 三种孔隙度测井联立解方程（矩阵方法，计算的核心是求逆阵）接出，Vsh , V1 , V2

42、 (2) 独立方法计算泥质含量，二种孔隙度测井组合，解出纯岩石模型的各个部分的体积。读入一个深度的测井数据孔隙度，矿物1和矿物2构成的体积模型泥质含量Vsh GR 取最小值作为泥质含量Vsh SPRtTHCNL-DEN泥质含量Vsh泥质含量Vsh泥质含量Vsh泥质含量Vsh 中子、密度测井曲线泥质校正Vsh计算含水饱和度Sw矩阵方法计算纯地层孔隙度，矿物1和矿物2的含量计算渗透率K读入测井解释参数最后一个深度点结束是否图3 计算步骤图先求出泥质含量，再校正泥质的影响后的测井数据联立求解, Vma1 , Vma2 ，相当按纯岩石计算，然后再根据泥质的多少换算出实际岩石的孔隙度和矿物含量。， （二）计算过程。（） 根据GR，SP求Vsh 当Vsh 50%时，按泥岩处理 =0，Sw =1.0, Vsh =100, V1 =0, V2 =0.K=0当Vsh <50%时，按下述方法处理。（2）进行泥质校正 （3）求纯岩石的孔隙度和矿物含量用矩阵表示为：CX=L其中：C= 则孔隙度和两种矿物的含量可以用解常数阵的逆的方法得出：X=C-1L（3）恢复成泥质岩石（） 计算含水饱和度（混合泥质公式）利用1963年Simandoux的（混合泥质公式）计算含水饱和度。可以看成是关于Sw的一元二次方程。对于A、B、C分别上述一元二次方程的系数。当a=1.0,