镜质组反射率反演的EASY_Ro数值模拟新方法
鄂尔多斯盆地渭北隆起奥陶系构造_热演化史恢复
第８ ８卷 第１ １期２ ０ １ ４年１ １月 地 质 学 报 ＡＣＴＡ ＧＥＯＬＯＧＩＣＡ ＳＩＮＩＣＡ Ｖｏｌ．８８ Ｎｏ．１１Ｎｏｖ． ２ ０ １ ４注：本文为国家自然科学基金项目（编号４１３７２１２８）、西北大学大陆动力学国家重点实验室（编号ＢＪ０８１３３－１）、国家重大专项（编号２０１１ＺＸ０５００５－００４－００７ＨＺ）和中国地质调查局科研项目（编号１２１２０１１３０４０３００－０１）资助的成果。
收稿日期：２０１４－０７－２９；改回日期：２０１４－０９－２５；责任编辑：周健。
作者简介：任战利，男，１９６１年生。
博士后，西北大学教授、博士生导师，主要从事盆地热史与油气成藏及油气评价研究工作。
通讯地址：７１００６９，陕西省西安市太白北路２２９号，西北大学地质学系；Ｅｍａｉｌ：ｒｅｎｚｈａｎｌ＠ｎｗｕ．ｅｄｕ．ｃｎ。
鄂尔多斯盆地渭北隆起奥陶系构造－热演化史恢复任战利１，２），崔军平１，２），李进步３），王继平３），郭科２），王维２），田涛２），李浩２），曹展鹏２），杨鹏２）１）西北大学大陆动力学国家重点实验室，西安，７１００６９；２）西北大学地质学系，西安，７１００６９； ３）苏里格气田研究中心，西安，７１００１８内容提要：鄂尔多斯盆地渭北隆起区构造位置独特，演化历史复杂。
该区下古生界奥陶系碳酸盐岩有机质丰度较高，是寻找天然气的有利地区。
奥陶系碳酸盐岩由于缺乏有效的古温标，热演化程度的确定及热演化历史的恢复一直是研究的难题。
本文利用渭北隆起奥陶系碳酸盐岩大量的沥青反射率测试资料，结合上覆晚古生代、中生代地层的镜质组反射率资料及磷灰石和锆石裂变径迹等古温标，恢复了渭北隆起的构造热－演化史。
研究结果表明古生界奥陶系热演化程度具有北高南低的特点。
奥陶系等效镜质组反射率普遍大于２．００％，处于过成熟干气阶段。
磷灰石裂变径迹资料表明渭北隆起抬升冷却具有南早北晚的特点。
南部奥陶系—下二叠统抬升早，约为１０２～１０７Ｍａ，北部自６５Ｍａ以来抬升，主要抬升时期为４０Ｍａ以来。
计算烃类成熟度史的实地TTI_R_o法
计算烃类成熟度史的实地 T T I - R o 法石 广 仁中国石油勘探开发研究院目前在盆地模拟中一般用镜质组反射率和干酪根 降解率表征烃类成熟度 ,它们均可作为生烃量史计算 的重要参数 。
此外的烃类成熟度指标不能用于计算生 应用 (1) 式时必须注意 ,其积分路径是沿着地史模 拟所得的关于某地层底界的埋藏史 (见图 1) ;还必须注 意 ,如果某地质时期地温不变或下降 ,干酪根热降解作 用暂停 ,就不计算该时期的 T T I ,直到地温回升至曾有 过的最高温度时再计算 。
烃量 1 ,2 ,故除非对专门研究热史有特殊价值 ,否则盆 地模拟一般不采用 。
目前重建烃类成熟度史主要有 3 种方法3 , 各自 适用范围为 : ①T T I (时间温度指数) - R o ( 镜质组反射 率) 法 。
需用 R o - 深度曲线求取 R o ,适用于勘探程度 较高的地区 , 但目前一般只制作单个 T T I - R o 公式 , 不仅使 R o 史的误差较大 ,还有明显的概念错误 。
②化 学动力学法 。
不需用 R o - 深度曲线求取 R o , 适用于 勘探程度中等的地区 ,但需要有热裂解实验资料 ,影响 了其应用范围 。
③E asy R o 法 。
不需用 R o - 深度曲线 求取 R o ,适用于勘探程度较低的地区 , 由于只采用通 用的化学动力学参数 ,不考虑实地差异 , 也影响 R o 史 的精度 。
针对 T T I - R o 法存在的问题 ,本文提出实地 T T I 图 1 模型井的埋藏史和 T T I 史及 R o 史剖面图2 实地 T T I - R o 公式的制作根据单井 m 个地层底界的现今 T T I 值以及 R o - 深度曲线上对应于各底界埋深最大时的 R o 值 ,得到的T T I - R o 曲线分为 m 段 , 每段 的 对 数 表 达 式 系 数 不同 。
- R o 法 ,即建立不同地质时期的 T T I - 概念正确 ,而且确保了 R o 史的精度 。
SYT5124-95沉积岩中镜质组反射率测定方法
SY 中华人民共和国石油天然气行业标准SY/T 5124—1995 ───────────────────────────────沉积岩中镜质组反射率测定方法1995—12—25发布 1996—06—30实施───────────────────────────────中国石油天然气总公司发布前言镜质组反射率是石油地质勘探中研究生油岩成熟度及古地温变化的常规分析参数，为了适应90年代专业技术的新发展，便于国内外技术交流，尽快与国际通用标准接轨，对SY 5124—86《有机质中镜质组反射率测定方法》进行了修订。
将原标准名称改为《沉积岩中的镜质组反射率测定方法》，适应范围扩大至对干酪根以及碎屑岩、碳酸盐岩和煤岩等全岩中的镜质组反射率测定，并对其中的光片制备、仪器调节、组分识别以及数据的精密度等内容进行了补充或修改，提高了标准的实用性、科学性和先进性。
本标准自生效之日起，同时代替SY 5124—86。
本标准由石油地质勘探专业标准化委员会提出并归口。
本标准起草单位：胜利石油管理局地质科学研究院。
本标准主要起草人：李佩珍本标准参加起草人：熊玉文王可仁张学军张敏尹玲张敏锋目次前言1 范围 (1)2 原理 (1)3 试剂材料及标样 (1)4 仪器设备 (1)5送样要求及光片制备 (2)6 测定对象及检测环境 (2)7 测定步骤 (2)8 数据处理及报告内容 (3)9 精密度 (3)中华人民共和国石油天然气行业标准SY/T 5124—1995代替SY 5124-86沉积岩中镜质组反射率测定方法───────────────────────────────1范围本标准规定了镜质组反射率的测定及数据统计处理方法。
本标准适用于岩石中富集的干酪根以及碎屑岩、碳酸盐岩和煤岩等全岩中镜质组反射率的测定。
2 原理镜质组反射率是指在波长546nm±5nm(绿光)处，镜质组抛光面的反射光强度对垂直入射光强度的百分比。
它是利用光电效应原理，通过光电倍增管将反射光强度转变为电流强度，并与相同条件下已知反射率的标样产生的电流强度相比较而得出。
沉积盆地热演化史研究方法
沉积盆地热演化史研究方法盆地热演化史研究方法很多，主要有地球动力学模型法及古温标法两类。
一、地球动力学模型法地球动力学模型法是通过对盆地形成和发展过程中岩石圈构造（伸展、减薄、均衡调整、挠曲形变等）及相应热效应的模拟（盆地定量模型），获得岩石圈热演化史（温度和热流的时空变化）。
不同类型的盆地，具有不同的热史模型，根据已知或假定的初始边界条件，通过调整模型参数，使得模型计算结果与实际观测的盆地构造沉降史相拟合，从而确定盆地底部热流史；进而结合盆地埋藏史，恢复盆地内地层的热演化历史。
不同类型的盆地由于其形成的地球动力学背景和成因机制的差异，导致盆地演化过程的不同。
因而描述其构造热演化过程的数学模型也是不同的，P．A．Allen和J．R．Allen（1990）在其论著中对岩石圈伸展作用形成的盆地、挠曲盆地及与走滑变形有关的盆地的热史模型都作过详细地论述。
（一）伸展盆地伸展盆地是目前研究较广泛、研究程度较高的盆地类型，裂谷、拗陷、拗拉槽和被动大陆边缘是其基本样式。
在地壳和岩石圈伸展、减薄作用下形成，其主要的构造热作用过程包括：岩石圈的伸展减薄、地幔侵位、与热膨胀和冷却收缩以及沉积负载相关的均衡调整。
裂谷是地壳中的拉张区，现代裂谷具有负的重力异常、高热流值和火山活动等特征，表明在深部存在某种热异常。
裂谷分主动裂谷与被动裂谷两种类型。
1978年McKenzie研究了被动裂谷或机械伸展模型的定量结论后，提出了瞬时均匀伸展模型。
该模型假定地壳和岩石圈的伸展量是相同的（即均匀伸展）；伸展作用是对称的，不发生固体岩块的旋转作用。
因此，这是纯剪切状态。
构造沉降主要取决于伸展量、伸展系数（β）以及初期地壳与岩石圈的厚度比值。
该模型可概括如下：①拉张盆地的总沉降量由两部分组成：其一是由初始断层控制的沉降，称为初始沉降，它取决于地壳的初始厚度及伸展系数β；其二是岩石圈等温面向着拉张前的位置松驰，从而引起的热沉降，热沉降只取决于伸展量的大小；②模拟结果表明，断层控制的沉降是瞬时性的，而热沉降的速率随时间呈指数减小，这是由于热流随时间减小的结果。
镜质组反射率反演的easy%ro数值模拟新方法
镜质组反射率反演的easy%ro数值模拟新方法近年来，随着激光测量技术的发展，镜质组反射率反演技术已经成为一种有效的测量方法。
为了更好的使用激光测量技术，研究人员们一直在寻求新的，可靠的方法。
最近，一项新的技术easy%ro数值模拟，被提出来被用于镜质组反射率反演。
它根据数值模拟技术，采用有限元方法，对镜质组反射率进行大量的计算，以准确估计其真实反射率。
easy%ro数值模拟的优点在于它的高度精确性，可以更准确的反演出镜质组反射率的真实反射率。
除此之外，easy%ro数值模拟还能通过有限元方法分析整个计算系统，全面把握每个元素对反射率反演的影响。
另外，easy%ro数值模拟还可以针对不同的光谱特性，设计出最佳的参数，从而有效改善反射率反演结果。
easy%ro数值模拟方法显著提高了镜质组反射率反演的精度和可靠性。
在实际使用中，它能够进行更为精确的反射率反演，从而获得更精确的结果，提高测量精度。
除此之外，easy%ro数值模拟方法还可以用于其他技术领域，如电磁场测量技术、红外高光谱技术、宽带可调谐激光器应用及其它光学技术等等。
例如，它可以用于可见光区域的液体成像技术，以及可见光谱参数定量分析技术，增强测量结果的精度和准确性。
总之，easy%ro数值模拟方法是一种新兴的激光测量技术，在提高镜质组反射率反演精度方面具有重大意义。
它可以用于现有的技术领域，提供准确、可靠的测量结果，从而大大增强测量精度和改善数
据质量。
同时，easy%ro数值模拟方法的计算复杂度和精度也会给研究者带来更多的自由，推动技术发展。
镜质组反射率反演的easy%ro数值模拟新方法
镜质组反射率反演的easy%ro数值模拟新方法由于现有的反射率反演算法仅能够有限地处理镜质组反射率反演，为此，本文提出一种基于easy%ro数值模拟的新方法来改善镜质组反射率反演。
这种新方法将easy%ro数值模拟方法应用于反射率反演算法，以有效地解决上述问题。
首先，直接建立反射率数据库，并使用一种有效的投影映射方法将空间中的物理量映射到空间中的反射率值。
其次，使用四次spline插值方法，以构建一个二元的查找表，表中的变量代表材料的参数（如温度、湿度等），以提高解算精度。
最后，利用多种场景下的反射率变化模式，对最终结果进行校验和优化，确保解决方案的稳定性和准确性。
该方法的主要组成部分是easy%ro数值模拟方法。
easy%ro方法利用线性反射率的估计和四次样条插值，以提高反演精度和稳定性。
通过对四次样条插值技术的应用，可以使查找表中的变量与反射率值之间的关系更加精确、灵活，从而提高解算准确性。
另外，线性反射率估计技术可以有效抑制反射率变化产生的误差，从而提高反射率反演精度。
此外，本文还提出了一种有效的反射率反演模型，该模型能够有效利用多种场景下的反射率变化模式，可以有效地优化反射率反演程序，使反射率反演结果更加准确。
为了验证该新方法的有效性，我们进行了一系列的实验，取得了较好的结果。
首先，我们在模拟实验中构建了一种新的镜质组反射率数据库，并用投影映射方法将物理量映射到反射率。
接着，四次样条插值技术被用来构建查找表，从而提高了解算精度。
然后，在真实场景测试中，本文提出的模型能够有效地利用多种场景的反射率变化模式，从而有效地优化反射率反演解算过程，使反射率反演解决方案更加准确。
综上所述，本文提出了一种基于easy%ro数值模拟方法的新镜质组反射率反演算法，能够有效地解决镜质组反射率反演的问题，且经过实验验证，该方法的有效性得到了显著的改善。
因此，本文提出的新方法有望成为一种有效的反射率反演解决方案，对改善反射率反演的精度具有重要意义。
实验一：镜质组反射率测定方法
镜质组反射率测定方法1.原理煤的镜质组反射率，是镜质组(在绿光中546nm)的反射光强对垂直人射光强的百分比。
测定时，是根据CCD所接收的反射光强与其光电信号成正比的原理，在显微镜下一定强度的人射光中，对比镜质组和已知反射率的标准片的光电信号值而确定的。
2.适用范围实验所用样品为粉煤光片(或块煤光片)，适用于单煤或混合煤，也基本上适用于沉积岩中分散有机质(镜煤色体和其他固体有机质)的反射率测定。
3.使用仪器和材料3.1 偏光显微镜3.1.1 测量光源功率不小于30W的钨卤素灯或钨丝白炽灯。
3.1.2 起偏器和检偏器应能装卸和旋转。
3.1.3 孔径光圈和视域光圈其中心和大小能调节，并能调节到同一水平光轴上。
3.1.4 物镜，无应变的油浸物镜。
一般放大倍数为x25至x60。
当测定特别微小的颗粒时，可采用倍数更高的油浸物镜。
3.1.5 目镜，观察目镜与测试目镜的放大倍数一般为x10。
观察目镜中应装有十字线和测微尺。
3.1.6 载物台垂直于显微镜竖轴，转动360度时，对中倍物镜的焦距无影响。
载物台上应装有机械推动尺，其X, Y轴的最大移动范围不小于20mm，物台测微尺最小刻度为1/100mm。
3.2 分光光度计3.3 标准片和浸油3.3.1 反射率标准选用与煤的反射率接近的一组合格的反射率标准片。
应保持标准片的表面光洁。
经常检查其反射率值，一般用一系列标准片相互比较检查。
3.3.2 零标准片它在浸油中的反射率小于10-6% .33.3.3 油浸液最好采用在23℃时折射率(546n m光中)为1.518。
4.实验步骤4.1 制样用试样压平器、胶泥、载片将光片固定在载片上。
4.2 调整仪器4.2.1仪器启动打开电源、灯和仪器有关的电器部件，并调到规定的数值上。
4.2.2调节显微镜光学系统4.2.2.1校正物镜中心，使其与显微镜竖轴一致。
4.2.2.2调节照明系统，其步骤如下－在已经压平的标准片(或试样)上滴以浸油，置于载物台上，并准焦。
镜质体反射率的研究现状
镜质体反射率的研究现状孔伟思;方石;袁魏;王凯;聂镜奇;侯继盛【摘要】镜质体的反射率与其受热升温正相关，该过程不可逆且十分稳定。
镜质体反射率最早用于煤岩学，其后作为重要的数字标尺被用于分析烃源岩的生烃模式。
本文详述了该方法的基本原理、实验过程、数据校正以及在盆地热历史等方面应用，特别论述了20世纪以来的发展过程。
其中该方法在热史方面的应用，由早期简单图版法，发展至现今的复杂计算模拟，其中经典的是 Lerche 计算法和 Easy％Ro 模拟法，已经可以较为准确的恢复盆地的热史和经历的最高古地温。
此外利用地层中 Ro值与深度之间的关系恢复剥蚀量，以及分析盆地的类型和排驱压力。
目前，镜质体反射率正在朝微观化和系统化两个方向发展，一是通过实测数据点的双峰或多峰直方图了解芳香苯环的区别和不同来源，定量解释其受热过程，二是发挥自身不可逆性及稳定记录的特点，与其他古温标法配合，相互约束揭示其受热过程。
%Vitrinite reflectivity varies directly as the heat it receives, which is irreversible and stable. The vitrinite reflectivity was firstly used in anthracology field and later was regarded as important digital staff to analysis models of source rock. In this paper, the fundamental principle, experimental processes and data correction of the above method were described as well as its application in studying the thermal history of a basin. The method applied in thernal field began by simple engraving method and now develops into using complex computational simulation. Lerche calculation method and easy Ro numerical simulation are considered as classical model because they recover more exactly either in the thermal history of a basin or in the respect of the highest temperaturethat the basin has experienced. Besides, we can recover the eroded thickness by using discrete Ro values in stratum. By the same way, we can also analyze basin's types and displacement pressure. At present, the vitrinite reflectivity develops into the directions of micro and systematicness.At first, we learn the difference and different sources of aromatic benzene by bimodal or unimodal histogram, explaining its heating process by quantitative interpretation.Secondly, we use the advantage of irreversibility and stability of recording basin thermal events. Cooperated with other geologic thermometers,they restrict each other in order to reveal its heating process.【期刊名称】《当代化工》【年(卷),期】2015(000)005【总页数】9页(P1020-1028)【关键词】镜质体反射率;原理;实验方法;热史;剥蚀量【作者】孔伟思;方石;袁魏;王凯;聂镜奇;侯继盛【作者单位】吉林大学地球科学学院，吉林长春130021;吉林大学地球科学学院，吉林长春 130021;中国石化江苏油田物探技术研究院，江苏南京 210046;吉林大学地球科学学院，吉林长春 130021;吉林省第二地质调查所，吉林吉林 132013;吉林省第二地质调查所，吉林吉林 132013【正文语种】中文【中图分类】TE122随着石油、煤炭等能源工业的不断发展，盆地的热演化过程逐渐被越来越多的研究人员重视。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
P 0. 31 0. 35 0. 32 0. 35 0. 35 0. 35 0. 31 0. 35 0. 34 0. 31 0. 35 0. 35 0. 31 0. 35 0. 36 0. 31 0. 35 0. 36
注: A —— EA SY % R 0 法; B ——T TI 法; C ——LO M 法。
沉积有机质的热演化反应遵循 Ar rhenius 一级 反应动力学方程: dw / dt = - A ·ex p( - E / RT ) 。 式中 w —— 残留组分的浓度;
A ——频率因子; E ——反应活化能; R ——通用气体常数; T ——绝对温度, 为时间的函数。
图 1 沉积有机质平行化学反 应活化能与计算因子的分布
考虑到有机质化学组成的不均衡性, 可将复杂 的热降解反应视为一系列具有相同频率因子的平行 反应。对于组分 i 的化学反应, 在温度 T 时下式成 立:
dw / dt= - w i·A ·ex p[ - E i/ R T ( t ) ] ;
∑ 并且 dw / dt= idw i/ dt,
式中 w i —— 第 i 个平行反应中残留组分的浓度。 由此, 反应程度 F 为:
∑ F = 1- w / w 0 = 1- f i( w i / w 0i) ,
式中 w 0——全部反应物的初始浓度; w 0i —— 第 i 个平行反应反应物的初始浓度; f i ——第 i 个平行反应的化学计算因子。 反应程度 F 可以通过把时间温度历史分解为一
系列恒温或等受热速率阶段而求得, 详细求算步骤 见文献[ 3] 。根据反应程度 F , 建立起模拟镜质组反 射率及其演化的数学模式:
CR C Press , Inc, 1983: 173～180 10 Lyon W S. Trace el ement meas urment s at t he coal-f ired s t eam
plant . O xford : Bl ackw el l scient if ic, 1997: 144 ( 收稿日期 1997—08—21)
J2 1. 15 2. 19 0. 85 0. 96 1. 73 0. 85 0. 76 1. 60 0. 62 0. 73 1. 22 0. 71 0. 78 1. 39 0. 65 0. 74 1. 15 0. 60
T 1. 00 1. 14 0. 82 0. 85 0. 95 0. 78 0. 69 0. 95 0. 78 0. 69 0. 68 0. 61 0. 72 0. 28 0. 65 0. 68 0. 10 0. 60
T he contents and dist ributions o f t race elements in co als and the enrichment coefficient s o f t race elements relative to the crust
abundance ar e ex po unded, and t he hazar do us and useful elements ar e evaluated pr eliminarily . Keywords t race element ; coal; r eg ular ities o f dist ribution; Do ngsheng M i探
第 26 卷
虑反应活化能变化、小分子平行反应等长处, 创立起 一种镜质组反射率数值模拟的新方法—— EASY% Ro 数值模拟法, 并在北美、西欧等地煤化作用、盆地 古地温和油气资源评价的研究中加以成功应用, 引 起国际能源地质界的高度重视[ 2, 3, 4] 。为了深化我国 煤化作用以及盆地古地温研究, 作者引入 EASY% Ro 方法, 并编制计算机模拟软件对山西南部上古生 界煤层的煤化作用历史进行了数值模拟, 取得了良 好效果。
针对上述存在的主要问题, 美国学者 Burnham 等( 1990) 吸取先前平行反应化学动力学模式法中考
3 钱丽君, 白清昭 等 . 陕 西北部 侏罗纪 含煤地 层及聚 煤特 征 . 西 安: 西北大学出版社, 1987: 128～167
4 孙景信, J ervis R E. 煤中微量元素及其在燃烧过程中的分布特征. 中国科学 A 辑, 1986; ( 12) : 1287～1294
2 EA SY % Ro 法原理与模拟方法
与 T isso t 法、L OM 法、T T I 法等类似, EASY% Ro 法的理论依据也是 A rrhenius 一级反应动力学原 理。但与 T issot 法相比, EA SY % R o 法在考虑平行反 应的同时, 极大地简化了一系列反应活化能呈高斯 分布的平行反应的数目, 简明而又较为全面地描述 了有机质类型、平行反应与受热条件之间的函数关 系, 使模拟方法更为简便实用: 与 L OM 法和 T T I 法 相比, EA SY % Ro 法考虑到平行化学反应与有机质 演化进程之间的关系, 使得模拟结果更接近于真实。
宋党育, 秦勇: 镜质组反射率反演的 EASY% R o 数值模拟新方法
·1 7·
表 1 山西南部上古生界煤镜质组反射率演化历史数值模拟结果
R o/ %
地质 时代
阳城—翼城 ABC
临汾—洪洞 ABC
安泽 A BC
沁源—沁县 ABC
霍州 A BC
汾西 ABC
K1 3. 34 4. 95 3. 70 1. 78 2. 45 1. 70 2. 60 3. 02 2. 50 2. 87 4. 53 2. 50 1. 15 1. 50 1. 10 1. 46 2. 20 1. 30
PRELIMINARY STUDY ON THE TRACE ELEMENTS IN COALS OF
DONGSHENG MINING DISTRICT, SHENFU COALFIELD
Do u T inghuan Xiao Daxian Do ng Yaqin Zhang Suian Zhang Qing ling
本区煤层经历了 4 个主要阶段: 第一阶段, 二叠纪至 三叠纪末期, 正常古地热场, 古地温梯度在 3℃/ 100
高北低、东高西低”以及低煤化烟煤—无烟煤均有赋 存的煤级展布格局。
m 左右; 第二 阶段, 早—中侏 罗世, 古地 温波 动阶
应予强调, 本区上古生界煤的煤化作用在燕山
段; 第三阶段, 晚侏罗世至早白垩世, 煤层受热温度 早期曾出现间断, 由此将整个煤化作用过程划分为
170～280℃, 古地温梯度普遍超过 6℃/ 100 m, 大地 热流密度普遍大于 96 mW / m2, 具有异常古地热场 的典型特征; 第四阶段, 晚白垩世至第四纪, 古地温 梯度 2. 03～3℃/ 100 m , 古地热场恢复正常[ 5] 。
根 据 上 述 前 期 古 地 热 史 研 究 成 果, 采 用 EA SY % Ro 对山西南部上古生界煤层镜质组反射率 演化 历史进行 了数值模 拟, 并 与实测 反射率 以及 T T I 法、LOM 法的模拟结果进行了对比, 结果见表 1。可以看出: EA SY % Ro 法结果与镜质组最大反射 率的实测值最为接近, T T I 法结果误差较大, L OM 法也有一定误差。因此 EA SY % Ro 法数值模拟结果 较为真实地反映出山西南部上古生界煤的煤化作用
* 国家自然科学基金资助项目( 编号: 49472125) 部分研究内容。
其中, Kar weil 图解法、时间—温度指数( T T I ) 法、有 机成熟度水平( LOM ) 法和 T issot 法在一定的历史 阶段中得以广泛应用。然而这些方法均存在某些不 足, 或在理论上不够严密而只能得出经验性结果, 或 考虑因素过于简单而使结果存在较大偏差, 或模拟 方法及要求条件过于繁杂而导致应用不便。
第 26 卷 第 3 期 1998 年 6 月
煤田地质与勘探 CO A L G EOL O GY & EX PL O RA T IO N
·1 5·
镜质组反射率反演的 EA SY% Ro 数值 模拟新方法*
宋党育 秦 勇 ( 中国矿业大学资源与环境学院 徐州 221008)
摘要 在 以往众多 镜质组 反射率 数值 模拟方 法的基 础上, 美国学 者在 90 年代 初创立 了 EASY% Ro 方法, 并成功地运用于北美和西欧的研究实践。本文在对其原理和方法进行介绍的同 时, 报道了应用该法对我国山西南部上古生界煤层镜质组反射率演化历史或煤化作用过程的数值 模拟结果。
5 廖自基编 . 环境中微量元素的 污染危害与迁移转化 . 北京: 科学 出版社, 1989: 97～108
6 王运 泉, 任德贻 . 煤中 微量元素研 究的进展 . 煤田地质 与勘探, 1994; 22( 4) : 16～20
7 王运 泉等 . 煤中微量 元素分布 特征初步 研究 . 地质 科学, 1997;
% Ro = exp( - 1. 6+ 3. 7F) 。 在此模式中, F 的变化范围在 0～0. 85 之间, 所以% Ro 的最大值小于 4. 7% 。由于该模式的建立过程中 未考虑有机质化学结构的各向异性, 故所得结果应 为镜质组平均反射率。