且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可不同于此处 的顺序执行所示出或描述的步骤。
[0053] 如图1所示，本发明实施例的技术思路为，构造响应方程（补偿声波、补偿中子、补 偿密度）的理论测井值，根据矿物质量含量结果和常规测井资料（补偿声波、补偿中子、补 偿密度）综合优化反演计算页岩地层矿物体积、干酪根体积和孔隙度大小。
[0054] 连施例一
[0055] 本实施例提供一种确定页岩地层干酪根含量的方法，W下结合图2对本实施例提 供的方法作详细说明。
[0056] 在步骤S201中，建立由骨架矿物、干酪根和孔隙组成的页岩地层体积模型。参考 图3,单位体积的页岩地层包括局部均匀的多种骨架矿物、干酪根和孔隙，其中孔隙包括气 体和水，满足
[0057]
[0058] 其中Vj为第j种骨架矿物体积含量，V|(为干酪根体积含量，为气体体积含量， 为水的体积含量，N为骨架矿物数量。Vi，…Vj，VK，Φβ，Φ?表示某一地层深度处的体积含 量。
[0059] 具体的，不同地区矿物成份并不相同，常见的骨架矿物包括；伊利石、绿泥石、蒙脱 石、石英、长石、云母方解石、白云岩、硬石膏和黄铁矿等。通常情况下，页岩地层中的娃质包 括石英、长石和云母，巧质包括方解石和白云岩，泥质包括蒙脱石、伊利石和绿泥石。
[0060] 在步骤S202中，基于所述体积模型构造页岩地层在不同体积含量参数条件下的 测井响应方程，确定响应方程误差和测量误差。
[0061] 页岩地层不同骨架矿物、干酪根和孔隙条件下的测井响应方程为：
[0062]
[006引其中为第i种测井曲线的响应结果；V,为第j种页岩地层矿物体积含量Λ为 干酪根体积含量；Φ,为气体体积含量，Φ?为水的体积含量;Ml,为第i种测井曲线对应第 j种矿物的骨架值；MAiK、Mig、分别为干酪根、气体和水对应的第i种曲线骨架值，N为 骨架矿物数量。
[0064] 在本实施例中，构造不同体积含量参数条件下的补偿声波测井、补偿中子测井和 补偿密度测井响应方程：
[0065] 补偿声波测井曲线响应方程
，
[006引其中MCj为第j种骨架矿物对应的声波骨架值;ACk、ACg、A(；分别是干酪根、气体 和水的声波骨架值；
[0067] 补偿中子测井曲线响应方程
[0068]
[0069] 其中MCNLj为第j种骨架矿物对应的中子骨架值；CNLk、CNLg、CNU分别是干酪根、 气体和水的中子骨架值；
[0070] 补偿密度测井曲线响应方程
[0071]
[007引其中MDENj为第j种矿物对应的密度骨架值；DENk、DENg、DEN"分别是干酪根、气体 和水的密度骨架值。
[0073] 由于测井仪器与数据采集系统均存在一定的测量误差，本实施例中某种测井响应 的测量误差取其实际测量数值大小的5%。
[0074] 任何测井方程都是根据对实际地层作一系列的数学物理简化之后建立的测井解 释模型得出的理论公式，响应方程中存在一定的响应方程误差τ1。
[00巧]在步骤S202中，响应方程误差包括：
[007引声波测井响应方程误差
，其中，δAtmf、δAti分别为 泥浆滤液和骨架矿物的补偿声波误差，V，为第j种骨架矿物体积含量，Φ为孔隙度，Φ= ，如上文所述，为气体体积含量，Φ"为水的体积含量；
[0077] 补偿中子测井响应方程误差：
[0078] 其中，δCNLmf、δCNLi分别为泥浆滤液和骨架矿物的补偿中子误差；
[0079] 补偿密度测井响应方程误差
[0080]其中，δPmf、δPi分别为泥浆滤液和骨架矿物的密度误差。
[0081] 在步骤S203中，采集实际测井曲线，W页岩地层矿物质量含量为约束条件，根据 所述测井响应方程，响应方程误差和测量误差建立测井解释的目标函数：
[0082]
[008引其中，
[0084]曰1为第i种测井曲线实际测量值，％;
[00财 X = (Vi,…Vj, V|(，4g，Φ J为页岩地层体积含量参数，^ ;
[008引Z为当前测井深度，m;
[0087] 01为第i种测井曲线的测量误差；
[008引 τ1为构造的第i种测井曲线响应方程误差；
[0089]fi(x，Z)为在Z深度时构造的第i种id井曲线响应；
[0090] m为测井曲线的种类；
[0091] g.，ω为X的第j种约束条件；
[0092]τi为第j种约束误差；
[0093]P为约束个数。
[0094] 在本实施例中，fi(x，Z)分别为补偿声波测井曲线响应AC，补偿中子测井曲线响应 C化，补偿密度测井曲线响应DEN。
[0095] 优选的，本实施例中将页岩地层的矿物质量含量作为约束条件：
[0096] 约束条件g,(x)表示为
其中，
[0097]Vj为第j种页岩地层骨架矿物的体积含量，Wj为页岩地层的矿物质量含量，Pb为 密度测井值，Pμα为矿物骨架密度。页岩地层矿物质量含量W，由全岩分析法或者ECS元素 俘获测井资料确定。
[009引相应的质量约束误差τ,取值为1。
[0099] 不限于此，还包括W下约束条件：
[0100] Vj>0,j为矿物种类；>0, >0。
[0101] 在步骤S204中，计算使所述目标函数取得最小值的最优体积含量参数，绘制最优 体积含量参数对应的测井响应的正演结果曲线。
[0102] 由目标参数F(x，a)可构建最优化测井解释的数学模型：
[0103]
[0104] 计算使所述目标函数F(x，a)取得最小值的最优体积含量参数x=(Vl，… V，，Vk，φ,，Φ?)，即可确定与实际体积含量最接近的体积含量参数。
[0105] 优选的，在本实施例中利用BFGS变尺度法确定搜索方向为馬。在搜索方 向苟上可寻找目标函数的极小值。
[0106] 其中正定对称矩阵馬为尺度矩阵，京=w持A)表示目标函数F(x，a)在第k次迭 代点处的梯度，在第k次搜索方向上任意一点表示为χ= 目标函数F(x，a)转换 为单变量函数Φ(t)。
[0107] 利用DSP0W抛物线插值法求取单变量函数Φ(t)的极小值点t。，根据t。确定页岩 地层最优体积含量参数X= (Vi，…Vj，Vi(，4g，Φ?)。
[0108] 上文中BFGS变尺度法和DSPOW抛物线插值法是本领域技术人员公知的测井解释 最优化方法，可参考《高等学校教材；测井数据处理与综合解释》（雍世和著，2007年8月1 日出版）。
[0109] 在步骤S205中，检验所述测井响应的正演结果曲线是否与实际测井曲线匹配，将 匹配的正演结果曲线对应的体积含量参数中的干酪根体积含量作为确定结果。
[0110] 具体的，根据实际测井曲线、响应方程误差和测量误差绘制正演曲线的的置信区 间
，其中为第i种测井曲线实际测量值；0 1和τ1对应为 测量误差和测井曲线响应方程误差；
[0111] 当所述测井响应的正演结果曲线处于正演曲线的置信区间内时，所述检验所述测 井响应的正演结果曲线与实际测井曲线匹配，即说明反演得到的最优体积含量参数X= (Vi，...Vj，VK，4g，Φ")充分反应了页岩地层中的实际体积含量。
[0112] 此外，可根据目标函数最优值的大小检验测井响应的正演结果曲线是否与实 际测井曲线匹配。对于利用最优化方法计算的最优体积含量参数X%目标函数最优值 F*(x*，a)w〇。F*(x*，a)是衡量正演测J井曲线是否逼近实际测I井曲线的总标志。F*(x*，a)值 较小，说明测井响应的正演结果曲线与实际测井曲线匹配程度较好。
[0113] 综上所述，本发明实施例提供的方法可根据页岩地层矿物质量解释结果，结合常 规测井资料中的补偿声波、补偿中子和密度测井资料反演计算出页岩地层