牛哥工程数学报告.docVIP

CHINA UNIVERSITY OF PETROLEUM 工程数学大作业 班 级： 地质工程14-6班 姓 名： 韩磊 学 号： 2014211518 日期 2014.12.18 正交多项式最小二乘拟合：镜质体反射率和埋深的关系 1.问题描述 沉积岩中有机质的丰度和类型是生成油气的物质基础，但是有机质只有达到一定的热演化程度才能开始大量生烃。勘探实践证明，只有在成熟烃源岩分布的地区才有较高的油气勘探成功率。所以，烃源岩的成熟度是决定油气勘探成的关键。 成熟度是表示沉积有机质向油气转化的热演化程度。由于在烃源岩的演化过程中，烃源岩中的有机质的许多物理性质、化学性质都发生相应的变化，并且这一变化过程是不可逆的，因而可以应用有机质的某些物理性质和化学组成的变化特点来判断有机质热演化程度，划分有机质演化阶段。在常用的判断指标中，运用最多的是精致提反射率（Ro）。 镜质体（Vitrinite）是一组富氧的显微组分，由同泥炭成因相关的腐殖质组成，具镜煤（Vitrain）的特征。镜质体反射率（Ro）是镜质体反射光的能力，被认为是目前研究干酪根热演化和成熟度的最佳参数之一。干酪根的光学研究结果表明，其基本成分为镜质体碎片和非晶质有机物，主要来源于陆生高等植物碎片。干酪根的热解过程与镜质体的演化过程相符，镜质体以芳香环为核，带有不同的支链烷基。在热演化过程中，链烷热解析出，芳环稠和，出现微片状结构，芳香片间距逐渐缩小，致使反射率增大、透射率减小、颜色变暗，这是一种不可逆反应。 镜质体反射率与成岩作用关系密切，热变质作用越深，镜质体反射率越大。在生物化学生气阶段镜质体反射率为低值，即低于0.5%。随着埋藏深度而逐渐变化，在热催化生油气阶段和热裂解生湿气阶段，反射率增加较快，约从0.5%上升至2%；至深部高温生气阶段，反射率继续增大。因此，测定烃源岩中有机质或煤夹层的镜质体反射率，可以预测油气的分布。 本次实验目的是要想要了解镜质体反射率和埋深的关系，进行离散数据的最小二乘拟合，运用相关数据建立二者的函数表达式。 2.数值算法 (1).专业方面 本次实验使用琼东南盆地Ya21-1-1井的深度和镜质体反射率数据（表1）。根据资料，该口井钻遇了渐新世、中新世、上新世的地层，是古近纪和新近纪的沉积，距今约数十个百万年（Ma）。 表1 琼东南盆地Ya-21-1-1井深度与Ro数据 深度(m) Ro（%） 深度(m) Ro（%） 深度(m) Ro（%） 深度(m) Ro（%） 1829.00 0.36 2825.00 0.54 3727.00 0.91 4247.00 0.97 2030.00 0.39 3018.00 0.58 3786.00 0.80 4630.00 1.23 2231.00 0.42 3228.00 0.53 3901.00 0.84 2423.00 0.46 3435.00 0.69 4075.00 0.90 2624.00 0.41 3627.00 0.75 4193.00 1.04 ——引自琼东南盆地三口井样品的干酪根ESR和Ro数据 (2).数值分析方面 本次实验主要运用最小二乘法对离散数据进行二次拟合，建立二者的关系。其中，在拟合过程中主要运用由正交函数组构造最佳平方逼近元的方法。 在内积空间中，设，但，在中寻找一个元素，使得即，若存在，则称其为在上的最佳平方毕竟函数(最小二乘拟合曲线)。 由正交函数组构造最佳平方逼近多项式 原理：若函数组是关于点集和带权的正交函数组，即，则法方程 其中 ， 最佳平方逼近多项式为 构造关于点集的正交函数组，由构造关于点集的首1正交多项式 其中： 式中内积定义为： 离散数据最佳平方逼近，若使用正交函数组构造最佳平方逼近，则： 因此有 （3）计算过程 按照教材中P268算法7-3，在Matlab中编写子程序DiscreteOrthoFit（x,y,m,w,xx）。以表1中深度数据为，Ro值为，对应输入向量，代入相应程序中，即得到关于和的拟合函数关系，并绘制出对应的函数曲线。(如图1) 所求的2次拟合多项式为： 0.000000091635395*x^2 - 0.00028776445*x + 0平方误差为：0.035594 图1 Ro-Depth数值及其拟合曲线 从图中可以看出，离散的*点为实际的数据点，连续的曲线为拟合的结果，根据Matlab得出的结论可知其 =0.035594，函数趋势也与地质实际相符合，即随着深度的增加，镜质体反射率（Ro）增大。 实际用途：根据由该井数据得出