0　引　言

不同沉积环境下沉积的岩石矿物种类及其含量不同,其所含元素的种类及其含量也会有较大的变化[1]。地层中化学元素的种类及丰度与矿物种类及含量有着密切的关系,通过确定元素含量和矿物含量之间的转换关系,可以把元素含量转换成矿物含量[2]。元素俘获能谱测井(ECS,Elemental Capture Spectroscopy)是斯伦贝谢公司推出的新型地层元素测井仪器,它利用中子与地层各元素作用发生辐射俘获核反应时瞬发的伽马射线能量和数量不同,这些不同取决于特定的核,每一种核具有其特征的俘获伽马射线谱,因此,通过测量和分析俘获伽马能谱,确定地层的主要元素和含量[3-6]。ECS测井能从岩石成分角度解决岩性识别问题,对识别成分差异较大而颜色、结构、构造差异不明显的复杂岩性具有极其重要的意义[7-11]。

对用户请求进行全面深入的分析，必须要寻找并利用那些包含以下信息或知识的知识源：灾区相关信息（对灾区灾情信息、灾区地理位置、灾区天气情况及所需救灾物资等信息的约束）、救灾物资供应信息（对供应方位置、可供应物资等信息的约束）、应急物资运输信息（对运输工具、最优的运送路径、运送时间等信息的约束）。

目前,对ECS资料的处理主要采用斯伦贝谢公司提供的元素与矿物之间的经验转换系数,但在对研究区的ECS资料处理过程中,发现利用该转换系数的处理结果与薄片分析数据存在较大的误差,说明这套转换系数已不适应于该地区的ECS资料处理。针对这一问题,利用X衍射荧光和X衍射全岩实验分析数据建立新的氧化物闭合模型,通过借助最小二乘方法和广义逆矩阵求解线性方程组的方法求解,得到了该地区的经验转换系数,利用该套经验转换系数有效提高了该区ECS测井资料的处理精度,满足了现场的应用需要。

1　矿物含量的计算方法

1.1　元素含量的计算

根据X衍射和X荧光分析,确定研究区具有8种元素(Si、Al、K、Fe、Ca、Ti、Mg、S)、8种矿物(黏土矿物、菱铁矿、黄铁矿、石英、正长石、斜长石、方解石、白云石);而ECS测井测量得到的仅有Si、Ca、Fe、S、Ti等元素的产额曲线,因此,需要建立氧化物闭合模型确定8种元素的含量曲线[13]。

通过每个元素的产额Yi除以相对灵敏度因子Si,可以比较容易地求得元素的相对含量值,Si是可以在实验室确定的仪器常数。元素的相对含量与所需的绝对元素含量wi之间可以通过随深度变化的归一化因子F值相联系,即

(1)

式中,wi为地层中第i种元素的重量百分含量;F为随深度变化的归一化因子;Yi为地层中第i种元素的相对产额;Si为第i种元素的相对灵敏度因子。归一化因子F应满足闭合条件,即所有元素的重量百分含量之和为1。但是,因子F是一个非常复杂的函数,直接准确计算F值几乎是不可能的[2,12]。同时,只用俘获伽马射线谱较难确定碳、氧、钠和镁元素含量,因此,元素的闭合模型条件难以满足。所以,在没有碳和氧元素的情况下,采用了一个近似闭合模型,即将俘获伽马射线谱所能确定的元素转换成氧化物或碳酸盐矿物,使这些元素的氧化物和碳酸盐的质量百分数之和为1,每个深度点有特定方程

(2)

式中,Xi为元素i的氧化物或碳酸盐的质量与第i种元素的质量比,定义为元素i的氧化物指数;WK为利用自然伽马能谱测井确定的K的质量百分含量;WAl通过式(3)计算得到

WAl=wAl[1-XSiWSi-XCaWCa-

XMgWMg-1.99WFe]

(3)

利用式(3)变换可得转换系数

研究发现，D-二聚体在下肢深静脉血栓患者血液中所占比例明显高于健康人群[12]。这是由于其作为凝血酶激活与继发性纤溶的特异性分子标志物[13-14]。本研究结果显示，两组患者术后血液流变学指标均发生了改变，D-二聚体也异常升高，均出现了血栓形成倾向，但A组上述变化更明显。分析其原因可总结为以下两个方面：(1)血流动力学改变。对重度股骨粗隆间骨折患者，采用持续硬膜外麻醉能在一定程度上促进动脉血流，提高静脉排空率，对血栓的形成产生抑制作用[15]。(2)血液的高凝状态改变。采用硬膜外麻醉能有效降低术后患者骨折处血窦渗血的发生率，抑制局部小血管血流，这在一定程度上避免了术中及术后过多失血。

1.2　矿物含量求解

通过求逆矩阵,就可以得到用元素百分含量表示的矿物重量百分含量,即

E=C M

(4)

式中,E为元素重量百分含量构成的矩阵;M为矿物重量百分含量构成的矩阵;C为转换系数。

Herron等[15]采用数理统计中的因子分析法得出元素含量与矿物的转换关系为

本实验保留了运用工具药分析两种作用于心脏药物的实验方法，通过选用适宜的药物浓度使维拉帕米和普萘洛尔对离体心脏的抑制作用程度相近。在此基础上分别观察两药对异丙肾上腺素作用的影响，并得出两药机制不同的结论。该实验成功率高，可重复性强，适合在学生实验中进行。

M=C-1E

(5)

式中,C-1为转换系数矩阵的逆矩阵。

“无规矩不可成方圆”。党中央的“八项规定”“三严三实”及“自律准则”“处分条例”等就是最好的规矩。千里之堤毁于蚁穴，一些巨贪，往往就是从贪小便宜步步为“盈”“炼”成的。“易鱼”类似事件的严肃处理，体现的是“抓早抓小”，防微杜渐，以遏制由量变到质变。“好规矩”需善始善终，严要求产生好成效，严管理养成好习惯，严标准才能得民心。杜绝“易鱼”类似案例发生，打虎学武松拳头要硬的同时，拍蝇力度亦不能松。

(6)

式中,aij为第j种矿物中第i种氧化物的质量分数;xj为岩石中第j种矿物的质量分数;yi为岩石中第i种氧化物的质量分数。

利用求解矿物含量模型的解法[16]对转换系数矩阵的方程进行求解,得到研究区的转换系数见表3。

(7)

2　利用X衍射数据刻度转换系数

显效：常规治疗或者康复治疗之后，基本恢复了各项项肢体功能和语言功能；有效：各项肢体功能有明显改善且可以简单对话；无效：肢体功能和语言能力没有明显的效果。通过对两组患者的治疗调查发现，护理满意度高达78%，护理满意度=（非常满意+满意）/总例数×100%。

在上述求解矿物含量矩阵方程的过程中,需要已知元素含量以及转换系数。元素含量可根据ECS测井资料获取,而转换系数矩阵的准确获取则为ESC资料求取矿物含量的关键,通常不同地区具有适合于该地区的地区转换系数。为此,需利用X衍射数据反演该地区的转换系数。

对研究区的井取样进行X衍射、X荧光测试,根据实验分析资料和现场资料确定主要矿物种类及其含量,结合测井得到的元素资料,确定主要元素种类,并整理计算得到所需的元素含量。对所选7个样品进行X衍射全岩分析,得到该地区的矿物成分及含量,主要由石英、斜长石、正长石、方解石、白云石、黏土矿物、黄铁矿、菱铁矿8种矿物组成(见表1);7个样品中的5个样品的元素种类及含量则通过X-射线荧光光谱分析确定(见表2)。

闭合模型中不考虑镁的影响,可以利用光电吸收截面指数确定镁的含量[14],利用元素俘获测井得到的元素产额结合建立的氧化物闭合模型利用式(2)求得F值,再利用式(1)即可求取其他元素的含量曲线,最终得到8种元素的相对含量曲线。

假定岩石中有n种矿物,每种矿物包含m种氧化物成分,则求矿物含量xj(j=1,…,n)的问题,可转化为解约束方程组问题

C=E-1M

(8)

记x=(x1,…,xn),为求得带约束方程组(6)的解,可将式(6)转化成极值函数的最优化求解问题

根据传统CORDIC算法的原理，把初始相位角转换到[0,π/4]范围内，则可以跳过n=0级的迭代.设初始向量经过CORDIC算法N次迭代后，此时向量的坐标为(XN,YN)，则幅度可表示为相位可表示为θ=ZN+arctan(YN/XN).因此，可定义幅度修正因子剩余未旋转角度Δθ=arctan(YN/XN).由于剩余未旋转角度可表示为易知其最大值为：

表1　研究区X衍射矿物成分样本点分析数据表

样本序号矿物百分含量/%黏土矿物菱铁矿黄铁矿石英正长石斜长石方解石白云石1228880000006843617000000590638722949643310006216625060332117000032347531630007067426880000000004650700016953670310375387984329159351161000000501501384128428455118686178800005815513013038215426500142072142000072546332540333647104191332

表2　研究区X荧光元素成分样本点分析数据表

本序号元素相对百分含量/%SiAlKFeCaTiMgS2436924303202700600100700033370963311232011033079000431577082341602430160961845359954318106922101102300073452537177121513020021139

表3　研究区X衍射标定转换反演参数表

矿物元素/%AlSiFeTiCaMgKS石英01468501正长石103210斜长石123042方解石0240037白云石21129黏土矿物1322240606050716黄铁矿46555345菱铁矿482

图1　X-T4井矿物含量处理成果图 *非法定计量单位,1 ft=12 in=0.304 8 m,下同

3　处理结果分析

图2　标定前后计算石英-长石含量误差分析图

图3　标定前后计算黏土矿物含量误差分析图

图4　标定前后计算碳酸盐岩含量误差分析图

利用X衍射数据得到的转换系数对该地区的ECS测井资料进行处理。图1为X-T4井的矿物含量处理成果图。从图1中可以看出,利用岩心X衍射资料确定的转换系数计算的矿物含量与斯伦贝谢公司的转换系数确定的矿物含量曲线的变化趋势基本一致,但是前者与薄片鉴定的矿物含量数据更加吻合。分别将两者计算的矿物含量进行交会图分析(见图2、图3、图4)。可以看出,使用X衍射数据反演的转换系数求得的矿物含量与薄片分析矿物含量的数据点较均匀地分布在对角线附近,而使用斯伦贝谢的经验转换系数求得的矿物含量与薄片分析矿物含量的数据点则较分散。总体上,由斯伦贝谢公司计算的石英、长石质矿物含量偏大,而碳酸盐岩和黏土矿物含量偏小。利用X衍射数据得到的转换系数求解得到的矿物含量值则更加接近地层实际矿物含量值,准确度更高。

4　结　论

(1) 元素俘获测井资料的处理中所利用的元素与矿物之间的转换系数具有较强的地区经验性,不同地区应建立适合本地区的转换系数以适应该地区的ECS资料处理。

(2) 利用岩心X衍射数据,通过最小二乘法求解,可以建立地区性元素与矿物的转换系数。

(3) 利用该方法得到的转换系数经实际资料的处理,并经检验后表明,该转换系数更加适合于该地区元素俘获测井资料的处理,处理精度明显提高。

参考文献:

[1]　章海宁, 张金功, 岳爱忠, 等. 利用地层元素识别沉积岩岩性和矿物含量计算方法 [J]. 测井技术, 2016, 40(6): 683-688.

[2]　孙建孟, 姜东, 尹璐. 地层元素测井确定矿物含量的新方法 [J]. 天然气工业, 2014, 34(2): 42-47.

[3]　袁祖贵, 楚泽涵. 一种新的测井方法(ECS)在王庄稠油油藏中的应用 [J]. 核电子学与探测技术, 2003, 23(5): 417-423.

[4]　袁祖贵, 成晓宁, 孙娟. 地层元素测井(ECS)——一种全面评价储层的测井新技术 [C]∥北京核学会核应用技术学术交流会, 2004: 208-213.

[5]　刘绪纲, 孙建孟, 郭云峰. 元素俘获谱测井在储层综合评价中的应用 [J]. 测井技术, 2005, 29(3): 236-239.

[6]　程华国, 袁祖贵. 用地层元素测井(ECS)资料评价复杂地层岩性变化 [J]. 核电子学与探测技术, 2005, 25(3): 233-238.

[7]　龚劲松, 杨鸣宇, 王静, 等. ECS元素测井技术在非常规储层评价中的应用 [J]. 油气藏评价与开发, 2014(2): 76-80.

[8]　程华国, 袁祖贵, 刘宁. 用地层元素测井资料确定储层黏土含量 [J]. 中国石油大学学报(自然科学版), 2004, 28(2): 28-30.

[9]　程希, 任战利. 利用地球化学元素与矿物关系识别GS油田下干柴沟组地层岩性 [J]. 地球物理学进展, 2008, 23(6): 1903-1908.

[10] 魏国, 赵佐安. 元素俘获谱(ECS)测井在碳酸盐岩中的应用探讨 [J]. 测井技术, 2008, 32(3): 285-288.

[11] 张锋, 刘军涛, 冀秀文, 等. 地层元素测井技术最新进展及其应用 [J]. 同位素, 2011, 24(b12): 21-28.

[12] 孙中春, 王振林, 王伟. 复杂岩性储层元素俘获测井评价方法研究 [J]. 吉林大学学报(地球科学版), 2015(增刊1).

[13] 庞巨丰. 地层元素中子俘获伽马能谱测井解释理论和方法 [J]. 测井技术, 1998, 22(2): 116-119.

[14] 李舟波, 王祝文. 科学钻探中的元素测井技术 [J]. 地学前缘, 1998, 5(1-2): 119-130.

[15] HERRON M M. Mineralogy from Grochemical Well Logging [J]. Clay and Clay Minerals, 1986, 34(2): 204-213.

[16] FENG Z, LI X T, WU H L, et al. Multimineral Optimization Processing Method Based on Elemental Capture Spectroscopy Logging [J]. Applied Geophysics, 2014, 11(1): 41-49.