CN103838984A - 一种自适应调整区域性解释参数的方法 - Google Patents
一种自适应调整区域性解释参数的方法 Download PDFInfo
- Publication number
- CN103838984A CN103838984A CN201410122200.4A CN201410122200A CN103838984A CN 103838984 A CN103838984 A CN 103838984A CN 201410122200 A CN201410122200 A CN 201410122200A CN 103838984 A CN103838984 A CN 103838984A
- Authority
- CN
- China
- Prior art keywords
- interpretation
- value
- parameters
- regional
- interpretation parameters
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
一种基于最优化测井解释技术的自适应调整区域性解释参数的方法。其特征包括：A、确定地层中矿物组分的名称和相应区域性解释参数值域范围；B、求得泥质区域性解释参数的调整顺序和调整方向；C、基于调整顺序和调整方向自动调整泥质区域性解释参数；D、调整的终止条件是理论计算曲线值与实际测井值的差值达到最小或者泥质解释参数值超出给定范围。该方法不仅能减少区域性解释参数的调试时间，而且能提高最优化测井解释成果的质量。
Description
技术领域
本发明涉及一种基于最优化测井解释技术的自适应调整区域性解释参数的方法，其特征是能降低最优化处理软件中区域性解释参数的调试难度和提高解释成果的质量。
背景技术
目前，最优化测井解释技术是根据地球物理学广义反演理论，以环境影响校正后反映地层特征的实际测井值为基础，根据适当的解释模型和测井响应方程，通过合理选择的区域性解释参数与储集层参数初始值，反算出相应的理论测井值，并与实际测井值作比较，按非线性加权最小二乘原理建立目标函数，用最优化技术不断调整未知储集层参数x，使目标函数达到极小值。一旦两者充分逼近，则此时计算理论测井值所采用的未知量x就是反映实际储集层参数值，即最优化测井解释结果x*。但上述的最优化测井解释技术存在以下两个方面的问题：（1）矿物解释模型不能完全描述地下地质情况，可能存在其他矿物对测井响应的影响，并且由于泥质的复杂性，选择固定的区域性解释参数会影响测井解释结果的可靠性；（2）由于将各测井值、储集层参数初始值、区域性解释参数值及误差等组成一个有机整体，它们之间的相互影响使解释参数的调试变得困难，资料处理所花费的精力和机时远多于常规分析程序。
发明内容
为了克服最优化测井解释技术中区域性解释参数调试时间长的不足，本发明提供一种自适应调整区域性解释参数的方法，该方法不仅能减少区域性解释参数的调试时间，而且能增加最优解的精确度。
本发明解决其技术问题所采用的技术方案是：A、提供关于地层中的多种矿物组分的信息，所述信息包括多种矿物组分的名称和相应区域性解释参数值域范围，建立测井解释模型。B、基于解释模型，在值域范围内任选一组区域性解释参数，利用最优化测井解释技术进行一次资料处理，得到目标函数的极小值并设定小于其平均值的预期值。C、根据各测井响应方程在全井段对目标函数极小值的总体贡献量按从大到小的顺序，确定各测井响应方程中泥质解释参数的调整顺序。D、对任意一个深度点，各测井响应方程中的泥质解释参数选择上一深度点的解释参数作为初始值；E、根据理论计算曲线值与实际测井值的大小确定泥质解释参数数值的调整方向；F、根据求得的测井响应方程中泥质解释参数的调整顺序和调整方向，对目标函数极小值大于预期值时的泥质区域解释参数在值域范围内从给定的初始值开始递增或者递减的调整；G、对调整过的区域性解释参数，利用最优化方法重新进行计算；H、自适应调整的终止条件是理论计算曲线值与实际测井值的差值达到最小或者泥质解释参数值超出值域范围。
其中，步骤A中区域性解释参数值域范围包括：判断测井响应方程的泥质解释参数是否为一定值，判断为否时，选择一定合理范围；其余矿物的区域性解释参数选择定值。
步骤E中泥质解释参数的调整方向为：理论计算曲线值大于实际测井值则减小解释参数数值，理论计算曲线值小于实际测井值则增加解释参数数值。
本发明的有益效果是，不再对区域性解释参数进行手动调整，可以减少程序调试时间的同时，还能提高最优化测井解释成果的质量。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1是本发明方法的流程图。
图2是常规人工调整区域性解释参数方法实施例的质检拟合图。
图3是自适应调整区域性解释参数方法实施例的质检拟合图。
图4是常规人工调整区域性解释参数方法实施例和自适应调整区域性解释参数方法实施例中的目标函数极小值变化图。
图5是密度测井响应方程的泥质区域性解释参数随深度变化图。
图6是声波时差测井响应方程的泥质区域性解释参数随深度变化图。
图7是中子测井响应方程的泥质区域性解释参数随深度变化图。
图8是自然伽马测井响应方程的泥质区域性解释参数随深度变化图。
具体实施方式
图1为本发明提供的一种自适应调整区域性解释参数方法的流程图，该方法步骤如下。
步骤（1）根据提供的地层信息选择出区域性解释参数，见下表，并建立测井解释模型。
| 石英 | 白云石 | 方解石 | 泥质 | 孔隙 | |
| DEN | 2.66 | 2.87 | 2.71 | 2.42（2.2-2.8） | 1 |
| DT | 55.5 | 43.5 | 47-49 | 125（70-150） | 189 |
| CN | -0.02 | 0.01 | -0.01 | 0.44（0.25-0.5） | 1 |
| GR | 4 | 2 | 7 | 200（150-200） | 15 |
| PE | 1.81 | 3.14 | 5.81 | 7 | 0.36 |
步骤（2）从表中选出一组解释参数。基于解释模型、响应方程、解释参数、储集层参数初始值和目标函数，对全井段进行一次最优化数据处理。对最优化处理后的理论计算曲线值与实际测井值进行拟合，见图2。
步骤（3）根据处理得到全井段的目标函数极小值的变化情况，设定图3实施例中目标函数极小值的预期值为20。
步骤（4）对任意一个深度点，各测井响应方程中的泥质解释参数选择上一深度点的解释参数作为初始值。
步骤（5）根据各测井响应方程在全井段对目标函数极小值的总体贡献量，按由大到小的顺序对各测井响应方程中泥质解释参数进行调整。图3实施例的调整顺序为CN、DT、GR、DEN。
步骤（6）根据理论计算曲线值与实际测井值的大小确定泥质解释参数数值的调整方向。
步骤（7）根据测井响应方程中泥质解释参数的调整顺序和调整方向，在全井段对目标函数极小值大于预期值时的泥质区域解释参数进行自适应调整。
步骤（8）采用调整过的解释参数，利用最优化方法重新计算。
步骤（9）将理论计算曲线值与实际测井值足够逼近或者泥质解释参数值超出给定范围作为自适应调整的终止条件。
步骤（10）优选出各测井响应方程的泥质解释参数和最优解。
图2实施例是采用人工调试的固定区域性解释参数进行最优化处理得到的测井曲线拟合图，实线曲线是实际测井曲线，虚线曲线是利用最优解计算的理论测井曲线。
图3实施例是采用本发明方法自适应调整的区域性解释参数进行最优化处理得到的测井曲线拟合图，实线曲线是实际测井曲线，虚线曲线是利用最优解计算的理论测井曲线。对比图3和图2可知，图3中理论测井曲线与实际测井曲线拟合的更好。
在图4中，虚线曲线是图2实施例的部分井段的目标函数极小值的变化曲线，实线曲线是图3实施例的部分井段的目标函数极小值的变化曲线。大部分深度点的实线曲线值均小于虚线曲线值，并且实线曲线值基本上小于预期设定值20。
图5、图6、图7和图8分别是图3实施例中密度、声波时差、中子和自然伽马测井响应方程的泥质区域性解释参数随深度变化图。其变化趋势符合实际地层的沉积模式。
综上所述，本发明的有益成果是：提供了一种自适应调整区域性解释参数的方法，通过创造性地在测井解释最优化处理中自动校正泥质解释参数，不仅能减少区域性解释参数的调试时间，而且能增加最优解的精确度。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想。同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
Claims (3)
1.一种自适应调整区域性解释参数的方法，包括：
A、提供关于地层中的多种矿物组分的信息，所述信息包括多种矿物组分的名称和相应区域性解释参数值域范围，建立测井解释模型；
B、基于解释模型，在值域范围内任选一组区域性解释参数，利用最优化测井解释技术进行一次资料处理，得到目标函数的极小值并设定小于其平均值的预期值；
C、基于各测井响应方程在全井段对目标函数极小值的总体贡献量按从大到小的顺序，确定各测井响应方程中泥质解释参数的调整顺序；
D、对任意一个深度点的各测井响应方程中的泥质解释参数，选择上一深度点的解释参数作为初始值；
E、基于理论计算曲线值与实际测井值的大小，确定泥质解释参数数值的调整方向；
F、基于求得的测井响应方程中泥质解释参数的调整顺序和调整方向，对目标函数极小值大于预期值时的泥质区域解释参数在值域范围内从给定的初始值开始递增或者递减的调整；
G、对调整过的区域性解释参数，利用最优化方法重新进行计算；
H、自适应调整的终止条件是理论计算曲线值与实际测井值的差值达到最小或者泥质解释参数值超出值域范围。
2.根据权利要求1所述的一种自适应调整区域性解释参数的方法，其特征是，步骤A中区域性解释参数值域范围包括：判断测井响应方程的泥质解释参数是否为一定值，判断为否时，选择一定合理范围；其余矿物的区域性解释参数选择定值。
3.根据权利要求1所述的一种自适应调整区域性解释参数的方法，其特征是，步骤E中泥质解释参数的调整方向为：理论计算曲线值大于实际测井值则减小解释参数数值，理论计算曲线值小于实际测井值则增加解释参数数值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410122200.4A CN103838984B (zh) | 2014-03-30 | 一种自适应调整区域性解释参数的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410122200.4A CN103838984B (zh) | 2014-03-30 | 一种自适应调整区域性解释参数的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103838984A true CN103838984A (zh) | 2014-06-04 |
| CN103838984B CN103838984B (zh) | 2016-11-30 |
Family
ID=
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1389298B1 (en) * | 1999-09-21 | 2011-03-02 | ExxonMobil Oil Corporation | Determining optimal well locations from a 3d reservoir model |
| CN103375166A (zh) * | 2012-04-26 | 2013-10-30 | 中国石油集团长城钻探工程有限公司 | 一种确定地层中多种矿物组分含量的方法 |
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1389298B1 (en) * | 1999-09-21 | 2011-03-02 | ExxonMobil Oil Corporation | Determining optimal well locations from a 3d reservoir model |
| CN103375166A (zh) * | 2012-04-26 | 2013-10-30 | 中国石油集团长城钻探工程有限公司 | 一种确定地层中多种矿物组分含量的方法 |
Non-Patent Citations (4)
| Title |
|---|
| 潘妍 等: "一种支持向量机的动态自适应加权算法", 《烟台大学学报》 * |
| 石林 等: "自适应岩石力学参数测井解释模型及其在水力压裂设计中的应用", 《钻采工艺》 * |
| 肖亮 等: "最优化方法在复杂岩性储集层测井评价中的应用", 《断块油气田》 * |
| 黄布宙 等: "复杂泥质砂岩储层测井解释模型研究", 《石油物探》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102360502B (zh) | 一种自动基线校正方法 | |
| DE112013007286T5 (de) | Algorithmus zur optimalen ZSV-Konfiguration mithilfe eines gekoppelten Bohrloch-Lagerstätten-Modells | |
| CN103106307B (zh) | 激光切割机的横梁设计方法 | |
| CN105631078B (zh) | 天然裂缝性油藏自适应介质的数值模拟方法 | |
| CN109061752B (zh) | 一种含灰质地层的电阻率曲线校正方法 | |
| CN105093331B (zh) | 获取岩石基质体积模量的方法 | |
| CN108444884A (zh) | 基于Harris模型的油相相对渗透率曲线表征方法 | |
| CN106326545A (zh) | 确定水平井页岩气藏地层参数的方法和系统 | |
| CN103616723A (zh) | 基于avo特征的crp道集真振幅恢复方法 | |
| CN105738952A (zh) | 一种水平井区储层岩石相建模方法 | |
| CN116523086A (zh) | 基于长短时记忆深度神经网络的单井生产动态预测方法 | |
| CN103838984A (zh) | 一种自适应调整区域性解释参数的方法 | |
| CN106481315A (zh) | 陆上砂岩油藏单井可采储量快速确定模型及建立方法 | |
| CN103838984B (zh) | 一种自适应调整区域性解释参数的方法 | |
| CN111666665A (zh) | 一种提高碳酸岩动静态杨氏模量转换模型精度的方法 | |
| CN107025330B (zh) | 一种单箱多室宽梁桥横向分布系数的计算方法 | |
| CN113283182B (zh) | 地层压力预测分析方法、装置、介质及设备 | |
| CN103105778B (zh) | 一种工业过程仿真数学模型参数的估计方法 | |
| CN109469465A (zh) | 注水井参数的处理方法、装置及存储介质 | |
| CN106351625B (zh) | 一个油田注水开发动态系统分析方法 | |
| CN110991084B (zh) | 一种基于流线数值试井的储层渗透率计算方法 | |
| CN112419493B (zh) | 页岩储层三维属性模型建立方法及装置 | |
| CN107704646A (zh) | 一种致密储层体积改造后的建模方法 | |
| CN112558179A (zh) | 远场应力状态的智能反演方法 | |
| CN102130730B (zh) | 一种链路仿真的方法及装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161130 Termination date: 20200330 |