水文测井资料的精细解释方法研究.pdfVIP

·792· 中国地球物理2011 方法，而且由于这种方法的电极系没有聚流能力，很容易受井眼的影响。为此采用有限元方法来实现其迭 代反演，不仅解决了仪器响应函数变化的问题，而且还将井眼的影响考虑在内。 本区普通电阻率测井采用的是电位电极系。首先将地下岩层考虑为轴对称，并对其进行网格剖分，考 虑到岩相有时出现厘米级变化，对网格在纵向上按l厘米进行剖分：径向上按指数变化进行剖分，离井眼 越近剖分越密，并且考虑井眼的变化，使得剖分线与井眼重合。然后根据之前得到的岩层界面信息对每个 岩层在网格中进行分配，并且赋予每个岩层一个初始的电阻率值。随后采用有限元技术对初始模型进行计 算，得到一条初始电阻率测井曲线，并将所得到的曲线与实际曲线进行对比，根据曲线之间的差异调整每 个岩层的电阻率值，调整之后再进行有限元计算，重复这一过程，使得计算曲线与实际曲线之间的均方差 小于预设值，则停止迭代。这时所得到的每个岩层电阻率值就是每个岩层实际的电阻率值。当然根据OPTV 所得到的初步岩相也作为约束条件参与到实际处理过程中。 对于声波测井和自然伽马测井，首先根据实际测量仪器可以得到仪器响应函数，然后对每个岩层赋予 一个初始测井值，根据仪器响应函数与岩层测井值的褶积便可以得到一条初始的测井曲线，并将所得到的 曲线与实际曲线进行对比，根据曲线之间的差异调整每个岩层的测井值，调整之后再进行褶积计算，重复 这一过程，使得计算曲线与实际曲线之间的均方差小于预设值，则停止迭代。这时所得到的每个岩层测井 值就是每个岩层实际的测井值。当然根据OPTV所得到的初步岩相同样也作为约束条件参与到实际处理过 程中。 3．测井资料的最优化处理 由于该区测井方法相对比较少，为了能够充分地利用这些信息，采用最优化方法进行处理。首先建立 各种测井方法的响应方程，考虑岩石由基质有效孔隙度、裂缝孔隙度、泥质以及岩石骨架组成，这样便可 以得到各种测井方法的响应方程。 在上述测井响应方程的基础上，便可以建立最优化目标函数，通过最优化技术求解该目标函数便可以 得到包括基质有效孔隙度、泥质含量、裂缝孔隙度、泥质电导率和渗透率在内的各种物性参数。 4．测并资料的解释 在得到每个岩层的相应物性参数之后，结合之前从OPTV资料大致得到的岩相信息，便可以对岩层进 行综合评价。如果对于某一岩相，根据迭代反演方法得到的测井曲线与实际测井曲线差异在迭代数较大时 仍然较大，这时候就需要考虑初始岩相的假设有问题，对初始岩相进行调整，相应岩相的约束条件也随之 改变，根据新的岩相假设重新进行迭代反演，直至迭代反演方法得到的测井曲线与实际测井曲线之间差异 小于设定值。 EXCEL界面下进行人机交互解释。通过这一交互平台，解释人员可以根据自己的经验对岩层模型进行调 整，并根据调整后的模型重新计算，从而使得解释结论更合理。 5．结论 与石油领域的测井处理与解释相比，水文领域的测井处理与解释的井段要小的多，在这种情况下，采 用有限元方法对电阻率测井进行迭代反演变得可以实现，一般情况下，迭代次数达到4~5次便可以满足停 止迭代的条件。在得到每个岩层精确的岩相和物性资料之后，便可以根据相应关系式来确定水文参数、土 层参数和污染物与土层发生反应的界面反应参数，或者对这些参数的取值进行一定的约束限制，从而为下 一步对污染物在含水层中运移的定量模拟中参数的选取提供了重要的依据。 本研究由英国自然环境研究委员会(NERC)资助。 参考文献 eta1．ThePermo-TriassicSandstonesoftheCheshireBasin．Baseline 【1】Griffiths，K～J QualityReportsSeries[R]．British Commissioned GeologicalSurvey Report,2002． a1．Fluid inSandstones：TheOnshoreUKPermo-TriassicRedBed 【2】 R．D．Barker，etFlowandSoluteMovement Scqu明ce- Special SocietyofLondon，2006． Publication[M]