一种地球物理勘探利用测井曲线自动求取岩性曲线的方法
【技术领域】
[0001] 本发明属于石油地球物理勘探领域，具体涉及一种地球物理勘探利用测井曲线自 动求取岩性曲线的方法，针对测井约束地震反演问题，涉及应用测井资料进行岩性分析识 别技术，可应用于石油地球物理勘探中的岩性反演。
【背景技术】
[0002] 随着油田勘探开发的深入以及储层描述技术的发展，大的整装的构造油气藏很难 发现，目前重点研究的是复杂隐蔽性的岩性油气藏，地震反演技术在岩性油气藏的发现中 发挥了较大的作用。随着隐蔽岩性油气藏开发的进程，对地震反演的需求也越来越高，常规 的波阻抗反演技术在勘探时期能够满足简单的岩性识别的要求，但针对精细储层预测的需 求，波阻抗反演技术很存在一定的多解性，需要开展岩性反演，直接反演岩性数据体，为地 质人员提供直观的、多解性小的岩性数据体，提高岩性识别的精度。岩性反演效果非常好， 但进行岩性反演时，要求井中有岩性曲线，实际的生产中，只有小部分的探井才进行岩性录 井和取芯工作，才能有岩性曲线，其他的开发井等都没有没有岩性曲线，因此，需要对应用 测井曲线进行岩性曲线求取的技术开展研究工作。
[0003] 现有的应用测井曲线划分岩性的方法主要是应用伽马曲线来识别砂泥岩，伽马低 的认为是砂岩，伽马高的认为是泥岩，更准确的做法是应用伽马曲线重新构建一条泥质含 量曲线（V sh)，采用公式如下：（公式1)
[0004]
[0005] Bq :纯地层背景值
[0006] B0 = P sd · GRsd (或者 P ls · GRls)
[0007] P b、P sh、P sd、P ls分别是目的层、泥岩层、纯砂岩、纯石灰岩的体积密度；
[0008] GR、GRsh、GRsd、GRls分别是目的层、泥岩层、纯砂岩、纯石灰岩的自然伽马测井。
[0009] 求得泥质含量后，在通过取值来判断砂岩和泥岩等，该取值一般根据经验，比如泥 质含量小于50为砂岩，大于50为泥岩。
[0010] 该方法主要用伽马和密度曲线来求取泥质含量曲线，然后利用泥质含量曲线求取 岩性曲线，经过与实际的岩性录井资料对比分析，发现存在一些误差，包括这些纯泥岩和纯 砂岩等的自然伽马值的确定也有很大的人为性，因此求取的泥质含量曲线不是很精确，所 以有必要开展泥质含量曲线方面的精确求取研究，需要同时结合其他多条曲线进行岩性曲 线的求取
【发明内容】
[0011] 本发明的目的在于解决上述现有技术中存在的难题，提供一种地球物理勘探利用 测井曲线自动求取岩性曲线的方法，提高岩性曲线的精度。
[0012] 本发明是通过以下技术方案实现的：
[0013] 一种地球物理勘探利用测井曲线自动求取岩性曲线的方法，包括：
[0014] (1)加载测井曲线，包括伽马、自然电位、电阻率、中子和中子寿命曲线；
[0015] (2)分别对这5条曲线进行校正，得到校正好的曲线；
[0016] (3)从校正好的曲线上分别读取5条曲线的最大值GMAX1和最小值GMIN 1 ;
[0017] (4)求取曲线泥质含量相对值SH1 :
[0018]
[0019] 式2中的SHLG1为测井曲线，即为步骤（1)中的5条曲线，i = 1，2,3,4,5分别代 表伽马、自然电位、电阻率、中子和中子寿命曲线，即采用同一个公式计算各曲线的泥质含 量；
[0020] (5)计算泥质含量：
[0021]
[0022] W > > J
[0023] G⑶R为与地层年代有关的经验参数；
[0024] (6)求得最终的泥质含量 Vsh : Vsh = min (Vshl，Vsh2，…，Vsh5);
[0025] (7)应用步骤（6)求得的最终泥质含量Vsh和实测密度曲线Density联合求取岩 性曲线；
[0026] (8)输出岩性曲线。
[0027] 所述步骤⑵中的校正包括：
[0028] 奇异值消除校正和或深度校正和或环境校正。
[0029] 所述步骤（5)中，如果是新地层，则GCUR取3. 7,如果是老地层，则GCUR取2。
[0030] 所述步骤（7)是这样实现的：
[0031] 采用逻辑判断：
[0032] 当Vsh < 50且Density〈地层密度值（KG/M3)时候，岩性曲线的取值为数值1，代 表的岩性为砂岩；当V sh > 50且Density >地层密度值（KG/M3)时候，岩性曲线的取值为数 值〇,代表的岩性为泥岩；这样就构建了一条由数值1和〇表不的岩性曲线，1代表砂岩，〇代 表泥岩；
[0033] 所述砂泥岩密度门槛值是根据具体地层综合分析后确定的值。
[0034] 与现有技术相比，本发明的有益效果是：常规泥质含量曲线的计算方法，应用的测 井曲线种类少，容易受到特殊岩性段的影响，同时，纯的岩性地层密度值得读取存在人为因 素，导致求得的泥质含量曲线与实际岩性录井资料存在较大误差。
[0035] 本发明能够应用多条测井曲线，包括伽马、中子、自然电位、电阻率、热中子等曲线 求取泥质含量曲线，然后再应用泥质含量曲线和密度曲线联合求取岩性曲线，能够同时综 合6条测井曲线所反映的岩性信息，同时考虑了地层年代信息，因此该方法求得的岩性曲 线在精度方面得到了较大提高，能够在没有岩性录井资料的井中开展应用，尤其在大量没 有做录井的开发井中应用，可以快速准确的求得岩性曲线。
【附图说明】
[0036] 图1某井用公式1方法求得的泥质含量曲线与伽马对比图，可以看到由于特殊岩 性段的存在，导致整体求得的泥质含量曲线偏低，也就是说砂岩成分偏高了，该方法计算的 泥质含量曲线存在较大的误差。
[0037] 图2某井多条测井曲线及用本发明方法求得的泥质含量曲线，该方法求得的泥质 含量曲线比图1求得的泥质含量在精度上得到了提高。
[0038] 图3某井用常规方法和本发明所用方法求得泥质含量曲线对比图，泥质含量1为 常规方法求得，泥质含量2为本发明方法求得，可以看到泥质含量2曲线更精确，受特殊岩 性段的影响很少。
[0039] 图4某井用步骤⑦求得的岩性曲线（右1列，黑色为砂岩），最终求得的砂岩与伽 马曲线和泥质含量曲线有很好的对应关系（对应低伽马和低泥质含量），同时也受密度曲 线的影响，能