CN117967288A - 一种用于油气田的井下压力监测系统和监测方法 - Google Patents
一种用于油气田的井下压力监测系统和监测方法 Download PDFInfo
- Publication number
- CN117967288A CN117967288A CN202410385363.5A CN202410385363A CN117967288A CN 117967288 A CN117967288 A CN 117967288A CN 202410385363 A CN202410385363 A CN 202410385363A CN 117967288 A CN117967288 A CN 117967288A
- Authority
- CN
- China
- Prior art keywords
- value
- oil
- conductivity
- gas field
- monitoring
- 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.)
- Pending
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 79
- 238000000034 method Methods 0.000 title claims abstract description 67
- 230000008859 change Effects 0.000 claims abstract description 26
- 239000000243 solution Substances 0.000 claims description 58
- 238000009826 distribution Methods 0.000 claims description 44
- 239000011435 rock Substances 0.000 claims description 30
- 238000001514 detection method Methods 0.000 claims description 28
- 230000015572 biosynthetic process Effects 0.000 claims description 27
- 238000005755 formation reaction Methods 0.000 claims description 27
- 238000002347 injection Methods 0.000 claims description 27
- 239000007924 injection Substances 0.000 claims description 27
- 239000003129 oil well Substances 0.000 claims description 20
- 239000004215 Carbon black (E152) Substances 0.000 claims description 17
- 229930195733 hydrocarbon Natural products 0.000 claims description 17
- 150000002430 hydrocarbons Chemical class 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 13
- 230000005540 biological transmission Effects 0.000 claims description 12
- 238000002922 simulated annealing Methods 0.000 claims description 12
- 238000005457 optimization Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 claims description 4
- 238000003860 storage Methods 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 238000005553 drilling Methods 0.000 description 8
- 230000007423 decrease Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 5
- 238000013016 damping Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 238000001914 filtration Methods 0.000 description 2
- 230000008054 signal transmission Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
- E21B47/07—Temperature
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/30—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with electromagnetic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- Mining & Mineral Resources (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Electromagnetism (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种用于油气田的井下压力监测系统和监测方法,涉及电磁波技术领域,该系统包括数据收发单元、藏储探测单元和监测修正单元,该方法包括在油井管道上实时采集井下压力数据,并在地表布设注入电极,在地表下方地层内布设电位电极,利用获取的电位差数据输出地电阻率数值,并用于构建电导率模型,计算出电磁波在目标地层中的衰减度值,最后计算出地面接收设备对电磁波的实际接收值。本发明对油气田井下的压力监测数据进行实时采集,通过记录地层的电位差变化值输出地电阻率数值,进而推导计算电磁波在地层中的衰减度值用于修正地面的电磁波接收信号参数,具有可以提升电磁波对井下油气压力监测精确度的优点和有益效果。
Description
技术领域
本发明涉及电磁波技术领域,具体涉及一种用于油气田的井下压力监测系统和监测方法。
背景技术
井下压力测量在油气开采行业中扮演着至关重要的角色。它不仅是评估油气井状态、监控生产效率和保障作业安全的关键指标,还对优化开采策略、评估储层性能和预防环境问题具有重要意义。目前在垂直油井下使用电磁波对传感器探测数据进行压力数据传输已成为油井压力数据监测的新型方式。由于油气田中同时存在有油藏和气藏的储备,为了提升数据探测的覆盖率,目前在符合条件的油气田内均有对油井底部进行横向掘进的形式,水平井段能够穿过地层中的水平延伸区域,与传统垂直井相比,水平井段可以减少钻井井数和钻井长度,从而降低钻井成本,这对于原本难以开采的油气储层特别有效。
在对油气田进行资源压力监测的过程中,油藏和气藏可能存在相邻、重叠甚至混合的分布情况,导致在水平方向不同位置的地层介质存在电导率差异,从而影响地层介质的整体电磁特性。当在油气田内将油井进行横向掘进,并在水平方向的不同位置进行电磁波压力监测时,会因为电导率差异而导致电磁波的传播数据出现多种不同程度变化,从而导致地下介质整体的不均匀性增加,进而使得电磁波在地层传播时出现不均匀的信号干扰和衰减,影响电磁波信号传输深度的精确监测,使得对油气田的资源藏储区域进行储量压力监测的测得数据与实际储量数据存在明显偏差。
发明内容
本发明提供一种用于油气田的井下压力监测系统和监测方法,解决在使用电磁波进行油气田井下压力监测时,因地层不均匀性影响电磁波传输深度使得压力监测精确度不足的问题。
本发明通过下述技术方案实现:
一种用于油气田的井下压力监测系统,该系统包括:
数据收发单元:在油井管道上设置压力传感器对井下压力数据进行采集,输出并记录在当前深度下的测定压力值,将测定压力值的电信号转换为数字信号,对数字信号进行编码调制以生成调制信号,并设置发射设备用以将调制信号发送至地面接收设备以输出第一监测值;
藏储探测单元:在油气田地表布设注入电极用以注入地下电流,在地表下方地层内布设电位电极用以记录目标地层的电位差变化值,基于电位差变化值使用反演方法输出地电阻率数值用于基于拟合曲线构建电导率模型,并将目标地层按厚度均分为数量若干的岩层;
监测修正单元:基于藏储探测单元的电导率模型获取油气藏储区域所占的岩层分布均数,用以结合油气藏储的电导率均值计算出电磁波在目标地层中的衰减度值,通过衰减度值与第一监测值计算出第二监测值用以表示地面接收设备对电磁波的实际接收值。
在对油气田进行资源压力监测的过程中,油藏和气藏可能存在相邻、重叠甚至混合的分布情况,当在油气田内将油井进行横向掘进,并在水平方向的不同位置进行电磁波压力监测时,会因为电导率差异而导致电磁波的传播数据出现多种不同程度变化,从而导致地下介质整体的不均匀性增加,进而使得电磁波在地层传播时出现信号干扰和衰减,影响电磁波信号传输深度的精确监测,使得对油气田的资源藏储区域进行储量压力监测的测得数据与实际储量数据存在明显偏差。基于此,本发明提供一种用于油气田的井下压力监测系统和监测方法,解决在使用电磁波进行油气田井下压力监测时,因地层不均匀性影响电磁波传输深度使得压力监测精确度不足的问题。
进一步地,所述注入电极和电位电极的布设形式包括:所述注入电极包括第一注入电极和第二注入电极,所述电位电极包括第一电位电极和第二电位电极;将地表区域划分为数量均为若干且分别独立监测的第一探测区和第二探测区,在第一探测区和第二探测区的地表分别布设第一注入电极和第二注入电极;所述第一探测区和第二探测区分别沿地表下方深度方向顺次设置第一电极层和第二电极层,在第一电极层和第二电极层分别布设第一电位电极和第二电位电极。
进一步地,使用拟合曲线构建电导率模型的过程还包括:对地电阻率数据进行拟合以得到拟合曲线,设置用于空间插值的空间变量用以描述地层电导率分布的自变量,使用空间插值方法对拟合曲线进行插值,生成电导率模型的空间分布图,并对插值结果进行拟合程度评估。
进一步地,所述衰减度值的计算过程包括:设地表到油井管道底部之间的油气藏储区域所占的岩层总数和序数分别为Z和t,岩层的平均厚度为h,油气藏储所占岩层中的平均电导率为L,并设置衰减度值为A,衰减比例常数为k,电导系数为μ,
则目标地层中总的衰减度值计算式为:。
进一步地,所述第二监测值的计算过程包括:设电磁波传输时的频率表示为f,设第一监测值的电磁波型号的复数函数表示为E0(f),第二监测值的电磁波型号的复数函数表示为Er(f),衰减度值表示为A(f),则第二监测值计算式表示为:,e为数学常数。
一种用于油气田的井下压力监测方法,该方法包括:
步骤S1:在油井管道上实时采集井下压力数据,将采集到的压力数据转换为数字信号并对数字信号进行编码调制以生成调制信号,将调制信号发送至地面输出第一监测值;
步骤S2:在地表布设注入电极,在地表下方地层内布设电位电极,使用反演优化方法处理电位电极获取的电位差数据以输出地电阻率数值,基于地电阻率数值构建电导率模型;
步骤S3:将目标地层按厚度均分为数量若干的岩层,电导率模型获取油气藏储区域所占的岩层分布均数,并结合油气藏储的电导率均值计算出电磁波在目标地层中的衰减度值;
步骤S4:基于衰减度值与第一监测值计算出第二监测值用以表示地面接收设备对电磁波的实际接收值,并回收井下的第一监测值发送记录用于与第二监测值比对以评估误差范围。
进一步地,将所述反演优化方法设为模拟退火优化方法用以获取地电阻率数值,该方法包括:
步骤A1:初始化设置模拟退火参数,包括初始温度、终止温度和温度下降率,并随机生成一个用于表示电导率分布解的初始解用作搜索起始点,并预设迭代次数;
步骤A2:将测得的电位差变化值设为观测值,基于每轮迭代时初始的电导率分布解正演计算出的电位差变化值设为预测值,并基于观测值和预测值的差异构建能量函数;
步骤A3:朝向能量函数减小的方向对电导率分布解进行邻域搜索以寻找解的更新,并在邻域的搜索空间内设置随机扰动参数以增加局部最优解的找寻范围;
步骤A4:每次迭代结束后更新模拟退火参数的温度下降值,当达到迭代次数或终止温度时得到电导率分布解的最终解,根据电导率分布解得出的电导率数值计算出地电阻率数值。
进一步地,所述能量函数的计算过程包括:
设能量函数表示为B,随机扰动参数表示为σ,观测值表示为Q,预测值表示为P,且设观测值的数量和序数分别为M和i,设随机扰动参数的数量和序数分别为N和j,设置正则化系数α,则能量函数的计算式设为:,
其中正则化系数α用于平衡能量函数的数据与随机扰动参数的拟合平滑度,且随机扰动参数的数量N与电导率分布解的数量一致。
本发明与现有技术相比,对油气田井下的压力监测数据进行实时采集,通过记录地层的电位差变化值输出地电阻率数值,进而推导计算电磁波在地层中的衰减度值用于修正地面的电磁波接收信号参数,具有可以提升电磁波对井下油气压力监测精确度的优点,对油气田开采工作具有降低风险和优化生产策略的有益效果。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不构成对本发明实施例的限定。在附图中:
图1为本发明结构示意图;
图2为本发明流程框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例1
如图1所示,本实施例为一种用于油气田的井下压力监测系统,该系统包括:
数据收发单元:在油井管道上设置压力传感器对井下压力数据进行采集,输出并记录在当前深度下的测定压力值,将测定压力值的电信号转换为数字信号,对数字信号进行编码调制以生成调制信号,并设置发射设备用以将调制信号发送至地面接收设备以输出第一监测值;
藏储探测单元:在油气田地表布设注入电极用以注入地下电流,在地表下方地层内布设电位电极用以记录目标地层的电位差变化值,基于电位差变化值使用反演方法输出地电阻率数值用于基于拟合曲线构建电导率模型,并将目标地层按厚度均分为数量若干的岩层;
监测修正单元:基于藏储探测单元的电导率模型获取油气藏储区域所占的岩层分布均数,用以结合油气藏储的电导率均值计算出电磁波在目标地层中的衰减度值,通过衰减度值与第一监测值计算出第二监测值用以表示地面接收设备对电磁波的实际接收值。
油气田井下的油井管道通常是用于开采和输送油气的管道系统。这些管道系统包括了从地表到油气层的井筒,其中包括了钻井设备以及用于注入水或其他辅助物质的管道。油井管道的主要组成部分包括钻井管、套管、井眼和井底装置等。钻井管是用于钻井的管道,通过钻进地下油气层;套管则是用于保护井壁和维持井孔稳定的管道;而井眼是指井筒的孔眼部分,用于油气的生产和采集;井底装置则是位于井底的设备,用于控制油气的生产流量和压力。所述压力传感器为常用的油气田井下测压设备,如电阻式压力传感器、压电式压力传感器和应变片式压力传感器等,其均可以适应高温、腐蚀和高压环境,并且具有较高的精度和稳定性,适用于极端的油井环境。所述测定压力值主要包括油井底部压力、油井管道内压力和地层压力。所述油井底部压力即油井井底的地下压力情况,这是评估油气储层压力状态的重要指标之一。所述油井管道内压力即为监测油井管道内的压力情况,以了解油气在管道中的流动情况和压力变化。
将测定压力值的电信号转换为数字信号可通过模数转换器完成,模数转换器将模拟电压信号转换为数字表示的压力值,通常使用脉冲密度调制(PDM)或者脉冲宽度调制(PWM)等技术进行转换。获取到数字信号后,可以进行数字信号处理,例如滤波、放大、校准等操作,以确保信号质量和稳定性。对数字信号进行编码调制以生成调制信号主要为了将其转换为适合传输的调制信号。常见的编码调制技术包括脉冲编码调制(PCM)、频移键控调制(FSK)、相移键控调制(PSK)等。编码调制的选择取决于传输介质、传输距离、抗干扰能力等因素。地面接收设备接收到调制信号后,需要进行解调处理,将其转换为可供分析和记录的数字信号。解调处理通常包括解调器解调和数字信号处理。所述第一监测值为在油气田地面接收生成的未经衰减修正处理的初始压力监测值。所述注入电极用于向地下注入电流的电极,通常与地面电源或电流发生器连接,并将电流引入地下。注入电极的位置通常与电位电极相对应,两者之间的距离可以根据需要调整。所述电位电极用于测量地下电位差,通常与地面测量仪器保持连接,并测量地下电位差的变化。
在具体实施中,可对采集到的电位差变化数据进行预处理,包括去除噪声、校正、滤波等操作,以提高数据质量和准确性。选择适合的反演方法对电位差变化数据进行处理。常见的反演方法包括有限元法、有限差分法、反演算法等。这些方法能够利用电位差变化值反推地下岩层的电导率分布,并利用电导率和电阻率之间的倒数计算关系推导出地电阻率数值。将地层按厚度均分为若干岩层可以使地下结构的模型更加精细化和准确化,通过对每个岩层进行定量分析,可以得到更多的地质参数,如电导率、电阻率、密度等,有助于更精确地描述地层的电性和物性特征,用于估算出基于藏储探测单元的电导率模型获取油气藏储区域所占的岩层分布均数。在建立电导率模型的基础上,需要识别出地下的油气藏储层,在现有技术中主要通过电导率模型中的异常区域来识别,因为油气藏通常具有较高的电导率,与周围的岩层有所不同。所述油气藏储的电导率均值可通过地球物理勘探邻域对油藏和气藏的电导率记录来计算获取。所述衰减度值即表示电磁波在地下传播过程中受到油气藏储阻碍的衰减程度,即受到常规岩层影响以外的信号强度降低程度。所述第二监测值即表示通过第一监测值与衰减度值计算后得到的测量修正值,即电磁波地面接收设备的实际接收值。
进一步地,作为一种可行的实施方式,所述注入电极和电位电极的布设形式包括:所述注入电极包括第一注入电极和第二注入电极,所述电位电极包括第一电位电极和第二电位电极;将地表区域划分为数量均为若干且分别独立监测的第一探测区和第二探测区,在第一探测区和第二探测区的地表分别布设第一注入电极和第二注入电极;所述第一探测区和第二探测区分别沿地表下方深度方向顺次设置第一电极层和第二电极层,在第一电极层和第二电极层分别布设第一电位电极和第二电位电极。通过将电位电极分设于不同深度,可以实现对地下不同深度区域的电位差变化值进行监测。通过在不同深度设置电位电极,可以获得不同深度处的电位差变化值,提高对地下结构和油气藏储层的垂直分辨率,使得监测结果更加精细和全面。不同深度处的电位差变化值反映了地下岩层的电性差异。通过分析不同深度处的电位差变化值,可以识别出更多不同地层的特征。这样可以获取更全面的地下电流分布信息,有助于准确评估油气田多个区域位置的井下不同深度的地层压力情况,提供了储层压力的多层次信息。电位电极记录地下电流经过不同深度时的电位差变化值,这些数据反映了地层的电性参数。通过分析不同深度处的电位差变化,可以推断地下岩层的电导率、介电常数等物理特性,为油气储层的识别和评估提供重要依据。
进一步地,作为一种可行的实施方式,使用拟合曲线构建电导率模型的过程还包括:对地电阻率数据进行拟合以得到拟合曲线,设置用于空间插值的空间变量用以描述地层电导率分布的自变量,使用空间插值方法对拟合曲线进行插值,生成电导率模型的空间分布图,并对插值结果进行拟合程度评估。拟合曲线通常是根据观测到的地电阻率数据,通过拟合方法找到最佳的拟合曲线,以最佳地反映地下岩层的电性参数。地电阻率数据通常是离散的观测值,而空间插值可以通过对这些观测值之间的空间关系进行推断,填补数据之间的空白,使得可以在更广泛的区域内生成连续的电导率分布模型。在插值过程中,需要使用空间变量来描述地下岩层的电导率分布情况,这些空间变量可以包括地表坐标、深度、岩性、地形等因素,以及拟合曲线的参数值。根据所选择的插值模型和空间变量,利用已有的观测数据和拟合曲线进行计算,以估计未观测位置的电导率数值。插值方法会根据空间相关性和观测数据的分布情况,对未观测位置的数值进行预测。所述空间分布图用于展示地下岩层的电导率分布情况,反映了地下结构的特征和复杂性。
进一步地,作为一种可行的实施方式,所述衰减度值的计算过程包括:设地表到油井管道底部之间的油气藏储区域所占的岩层总数和序数分别为Z和t,岩层的平均厚度为h,油气藏储所占岩层中的平均电导率为L,并设置衰减度值为A,衰减比例常数为k,电导系数为μ,
则目标地层中总的衰减度值计算式为:。
对每一层岩层的电导率与该层的平均厚度进行乘积运算。这个乘积代表了每一层岩层对电磁波传播的衰减程度的贡献,即每一层的衰减度。将所有岩层的电导率与厚度的乘积进行累加求和,得到总的电导率与厚度乘积的累加值,这个累加值反映了所有岩层对电磁波传播的综合衰减程度。传输衰减度的计算式中出现根号是缘于电磁波在地下传播时,其信号强度会随着传播距离的增加而减弱,因此需要使用非负的电导率来描述地下介质的传导特性,而衰减度的计算式中的根号通常是为了确保衰减度是一个正值,以便表示信号的衰减程度,因此设为根式开方的形式。通常采用地电阻率法或电磁法进行地下储藏层的探测,通过测量传输衰减度以及相应的电导率分布,可以利用反演算法来确定油气藏区域的地层数和电导率分布。所述衰减比例常数k用于调节电磁波在地下传播过程中的衰减程度,可以控制衰减度值A的尺度,使其与实际情况相符合。衰减比例常数的选择基于电磁波在地下传播时的衰减特性和地下岩层的电性参数,通常情况下,地下岩层的电导率和介电常数会影响电磁波的传播衰减程度,而衰减比例常数可以根据地下岩层的特征来确定。在具体运用中,衰减比例常数的选取可以是根据实验数据、模拟计算或者经验公式来确定,通过对不同地质条件下的数据分析和对比,可以确定适合的衰减比例常数,以使得计算得到的衰减度值能够反映地下结构的真实特征。所述电导系数是描述物质导电性能的指标,表示单位长度内物质对电流的导电能力。高电导系数意味着物质具有良好的导电性,电流能够在物质中快速传播;低电导系数则表示物质对电流的阻力较大,导电性能较差。电导系数μ的作用是用来描述地下岩层的导电性能,它影响了电磁波在地下岩层中传播时的衰减程度。
进一步地,作为一种可行的实施方式,所述第二监测值的计算过程包括:设电磁波传输时的频率表示为f,设第一监测值的电磁波型号的复数函数表示为E0(f),第二监测值的电磁波型号的复数函数表示为Er(f),衰减度值表示为A(f),则第二监测值计算式表示为:,e为数学常数。
第二监测值表示结合电磁波在地下传输过程中的衰减效应后,实际接收到的电磁波信号强度,通过将第一监测值乘以衰减度的负指数,可以更准确地模拟实际地下电磁波信号的接收情况。该计算式表示将第一监测值的电磁波型号的复数函数E0(f)乘上衰减度值的负指数即 e^(-A(f)),得到考虑衰减效应后的第二监测值的电磁波型号的复数函数Er(f)。其中所述第二监测值的电磁波型号Er(f)表示其振幅和相位。所述衰减度值A的计算通常需要考虑地下岩层的电导率分布,可以通过地电阻率测量数据或者电磁法测量数据来获取。第二监测值表示了考虑到电磁波在地下传播过程中的衰减程度后的实际接收值。通过将第一监测值乘以衰减度的负指数,可以得到考虑衰减效应后的电磁波信号强度,即第二监测值。
实施例2
如图2所示,一种用于油气田的井下压力监测方法,该方法包括:
步骤S1:在油井管道上实时采集井下压力数据,将采集到的压力数据转换为数字信号并对数字信号进行编码调制以生成调制信号,将调制信号发送至地面输出第一监测值;
步骤S2:在地表布设注入电极,在地表下方地层内布设电位电极,使用反演优化方法处理电位电极获取的电位差数据以输出地电阻率数值,基于地电阻率数值构建电导率模型;
步骤S3:将目标地层按厚度均分为数量若干的岩层,电导率模型获取油气藏储区域所占的岩层分布均数,并结合油气藏储的电导率均值计算出电磁波在目标地层中的衰减度值;
步骤S4:基于衰减度值与第一监测值计算出第二监测值用以表示地面接收设备对电磁波的实际接收值,并回收井下的第一监测值发送记录用于与第二监测值比对以评估误差范围。
更多地,作为一种可行的实施方式,将所述反演优化方法设为模拟退火优化方法用以获取地电阻率数值,该方法包括:
步骤A1:初始化设置模拟退火参数,包括初始温度、终止温度和温度下降率,并随机生成一个用于表示电导率分布解的初始解用作搜索起始点,并预设迭代次数;
步骤A2:将测得的电位差变化值设为观测值,基于每轮迭代时初始的电导率分布解正演计算出的电位差变化值设为预测值,并基于观测值和预测值的差异构建能量函数;
步骤A3:朝向能量函数减小的方向对电导率分布解进行邻域搜索以寻找解的更新,并在邻域的搜索空间内设置随机扰动参数以增加局部最优解的找寻范围;
步骤A4:每次迭代结束后更新模拟退火参数的温度下降值,当达到迭代次数或终止温度时得到电导率分布解的最终解,根据电导率分布解得出的电导率数值计算出地电阻率数值。
所述初始温度是模拟退火算法开始时的温度值,它决定了搜索过程中的起始状态的随机性。终止温度是模拟退火算法搜索过程中的停止条件。当温度降低到终止温度以下时,算法停止搜索并输出当前解作为最终解。终止温度通常是一个很小的值,当温度降低到这个阈值以下时,认为算法已经收敛到最优解或者近似最优解。温度下降率控制着模拟退火算法中温度的下降速度。对电导率分布解进行邻域搜索时,从初始解开始作为搜索的起始点,初始解可随机生成或根据油气藏的电导率参数知识评估设定。通过对当前解进行变化设定后在当前解的邻域内生成新的解,这些变化设定根据具体实施的油气田参数特征,可以是微小的,也可以是较大的。后对生成的邻域解进行评估,计算其对应的能量值。能量值可以是目标函数的值,也可以是一种衡量解优劣的指标。如果新解比当前解更优,那么就接受新解。如果新解比当前解差,那么以一定的概率接受新解,这个概率会随着温度的降低而逐渐减小。如果接受了邻域解而将其设为当前解,即更新搜索过程中的相关参数。否则保持当前解的参数不变。重复以上步骤直到达到最大迭代次数,或者达到设定的终止温度为止并输出最终解。引入所述随机扰动参数,可以在每次迭代中对当前解进行随机的微小变化,从而增加了搜索的多样性,有助于跳出局部最优解,更好地探索解空间。
进一步地,所述能量函数的计算过程包括:
设能量函数表示为B,随机扰动参数表示为σ,观测值表示为Q,预测值表示为P,且设观测值的数量和序数分别为M和i,设随机扰动参数的数量和序数分别为N和j,设置正则化系数α,则能量函数的计算式设为:,
其中正则化系数α用于平衡能量函数的数据与随机扰动参数的拟合平滑度,且随机扰动参数的数量N与电导率分布解的数量一致。
能量函数用于衡量观测值与预测值之间的差异以及随机扰动参数的拟合平滑度,用于衡量解的优劣程度。对于每个观测值和预测值的差的平方用于量化观测值与预测值之间的差异。所述随机扰动参数平方项用于衡量随机扰动参数的平滑度。所述正则化系数用于对随机扰动参数的平方项进行加权,以平衡能量函数的数据项与平滑度项之间的影响,确保能量函数在衡量数据拟合程度和平滑度时达到合适的平衡。将观测值与预测值之间的差异项和随机扰动参数的平滑度项加权相加,得到能量函数B的总值。所述能量函数最终的目标即为通过优化能量函数,使得其值最小化。通过优化能量函数,可以得到适合的模型参数,从而提高模型的拟合效果和稳定性。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种用于油气田的井下压力监测系统,其特征在于,该系统包括:
数据收发单元:在油井管道上设置压力传感器对井下压力数据进行采集,输出并记录在当前深度下的测定压力值,将测定压力值的电信号转换为数字信号,对数字信号进行编码调制以生成调制信号,并设置发射设备用以将调制信号发送至地面接收设备以输出第一监测值;
藏储探测单元:在油气田地表布设注入电极用以注入地下电流,在地表下方地层内布设电位电极用以记录目标地层的电位差变化值,基于电位差变化值使用反演方法输出地电阻率数值用于基于拟合曲线构建电导率模型,并将目标地层按厚度均分为数量若干的岩层;
监测修正单元:基于藏储探测单元的电导率模型获取油气藏储区域所占的岩层分布均数,用以结合油气藏储的电导率均值计算出电磁波在目标地层中的衰减度值,通过衰减度值与第一监测值计算出第二监测值用以表示地面接收设备对电磁波的实际接收值。
2.根据权利要求1所述的一种用于油气田的井下压力监测系统,其特征在于,所述注入电极和电位电极的布设形式包括:所述注入电极包括第一注入电极和第二注入电极,所述电位电极包括第一电位电极和第二电位电极;将地表区域划分为数量均为若干且分别独立监测的第一探测区和第二探测区,在第一探测区和第二探测区的地表分别布设第一注入电极和第二注入电极;所述第一探测区和第二探测区分别沿地表下方深度方向顺次设置第一电极层和第二电极层,在第一电极层和第二电极层分别布设第一电位电极和第二电位电极。
3.根据权利要求1所述的一种用于油气田的井下压力监测系统,其特征在于,使用拟合曲线构建电导率模型的过程还包括:对地电阻率数据进行拟合以得到拟合曲线,设置用于空间插值的空间变量用以描述地层电导率分布的自变量,使用空间插值方法对拟合曲线进行插值,生成电导率模型的空间分布图,并对插值结果进行拟合程度评估。
4.根据权利要求1所述的一种用于油气田的井下压力监测系统,其特征在于,所述衰减度值的计算过程包括:设地表到油井管道底部之间的油气藏储区域所占的岩层总数和序数分别为Z和t,岩层的平均厚度为h,油气藏储所占岩层中的平均电导率为L,并设置衰减度值为A,衰减比例常数为k,电导系数为μ,
则目标地层中总的衰减度值计算式为:。
5.根据权利要求4所述的一种用于油气田的井下压力监测系统,其特征在于,所述第二监测值的计算过程包括:设电磁波传输时的频率表示为f,设第一监测值的电磁波型号的复数函数表示为E0(f),第二监测值的电磁波型号的复数函数表示为Er(f),衰减度值表示为A(f),则第二监测值计算式表示为:,e为数学常数。
6.一种用于油气田的井下压力监测方法,基于权利要求1-5任意一项所述的一种用于油气田的井下压力监测系统,其特征在于,该方法包括:
步骤S1:在油井管道上实时采集井下压力数据,将采集到的压力数据转换为数字信号并对数字信号进行编码调制以生成调制信号,将调制信号发送至地面输出第一监测值;
步骤S2:在地表布设注入电极,在地表下方地层内布设电位电极,使用反演优化方法处理电位电极获取的电位差数据以输出地电阻率数值,基于地电阻率数值构建电导率模型;
步骤S3:将目标地层按厚度均分为数量若干的岩层,电导率模型获取油气藏储区域所占的岩层分布均数,并结合油气藏储的电导率均值计算出电磁波在目标地层中的衰减度值;
步骤S4:基于衰减度值与第一监测值计算出第二监测值用以表示地面接收设备对电磁波的实际接收值,并回收井下的第一监测值发送记录用于与第二监测值比对以评估误差范围。
7.根据权利要求6所述的一种用于油气田的井下压力监测方法,其特征在于,将所述反演优化方法设为模拟退火优化方法用以获取地电阻率数值,该方法包括:
步骤A1:初始化设置模拟退火参数,包括初始温度、终止温度和温度下降率,并随机生成一个用于表示电导率分布解的初始解用作搜索起始点,并预设迭代次数;
步骤A2:将测得的电位差变化值设为观测值,基于每轮迭代时初始的电导率分布解正演计算出的电位差变化值设为预测值,并基于观测值和预测值的差异构建能量函数;
步骤A3:朝向能量函数减小的方向对电导率分布解进行邻域搜索以寻找解的更新,并在邻域的搜索空间内设置随机扰动参数以增加局部最优解的找寻范围;
步骤A4:每次迭代结束后更新模拟退火参数的温度下降值,当达到迭代次数或终止温度时得到电导率分布解的最终解,根据电导率分布解得出的电导率数值计算出地电阻率数值。
8.根据权利要求7所述的一种用于油气田的井下压力监测方法,其特征在于,所述能量函数的计算过程包括:
设能量函数表示为B,随机扰动参数表示为σ,观测值表示为Q,预测值表示为P,且设观测值的数量和序数分别为M和i,设随机扰动参数的数量和序数分别为N和j,设置正则化系数α,则能量函数的计算式设为:,
其中正则化系数α用于平衡能量函数的数据与随机扰动参数的拟合平滑度,且随机扰动参数的数量N与电导率分布解的数量一致。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410385363.5A CN117967288A (zh) | 2024-04-01 | 2024-04-01 | 一种用于油气田的井下压力监测系统和监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410385363.5A CN117967288A (zh) | 2024-04-01 | 2024-04-01 | 一种用于油气田的井下压力监测系统和监测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117967288A true CN117967288A (zh) | 2024-05-03 |
Family
ID=90861342
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410385363.5A Pending CN117967288A (zh) | 2024-04-01 | 2024-04-01 | 一种用于油气田的井下压力监测系统和监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117967288A (zh) |
Citations (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4209747A (en) * | 1977-09-21 | 1980-06-24 | Schlumberger Technology Corporation | Apparatus and method for determination of subsurface permittivity and conductivity |
US5345179A (en) * | 1992-03-09 | 1994-09-06 | Schlumberger Technology Corporation | Logging earth formations with electromagnetic energy to determine conductivity and permittivity |
CA2375432A1 (en) * | 1999-06-21 | 2000-12-28 | Teruhiko Hagiwara | Method for radial profiling of resistivity at multiple depths of investigation |
CN101050700A (zh) * | 2007-05-09 | 2007-10-10 | 中国海洋石油总公司 | 一种用高导磁材料的电磁波测井仪校验方法 |
CN101609169A (zh) * | 2009-04-08 | 2009-12-23 | 中国石油集团钻井工程技术研究院 | 一种提高电磁波电阻率测量精度和扩展其测量范围的方法 |
CN102606138A (zh) * | 2012-03-31 | 2012-07-25 | 中国电子科技集团公司第二十二研究所 | 随钻电磁波电阻率测井仪相位幅度介电常数校正方法 |
CN103015970A (zh) * | 2012-11-30 | 2013-04-03 | 中国海洋石油总公司 | 一种随钻电阻率测井仪的模拟检测设备 |
CN103645512A (zh) * | 2013-12-10 | 2014-03-19 | 江苏大学 | 一种多层积水采空区探测方法 |
CN104912544A (zh) * | 2015-05-29 | 2015-09-16 | 中国航天科技集团公司烽火机械厂 | 一种油气井井下温度压力参数探测系统及方法 |
CN105700026A (zh) * | 2016-01-26 | 2016-06-22 | 中国科学院地质与地球物理研究所 | 一种地下目标体导电指数提取方法与装置 |
CN106054257A (zh) * | 2016-05-11 | 2016-10-26 | 中国科学院地质与地球物理研究所 | 一种含磁层状介质电磁响应的确定系统 |
CN107605475A (zh) * | 2017-10-27 | 2018-01-19 | 罗淮东 | 用于地层测试的设备、系统及方法 |
CN109061751A (zh) * | 2018-06-20 | 2018-12-21 | 西安石油大学 | 多层介质的视电导率计算方法 |
CN109252841A (zh) * | 2018-11-16 | 2019-01-22 | 中国海洋石油集团有限公司 | 基于电磁波无线监测提高注聚井注入能力的系统和方法 |
CN109752762A (zh) * | 2018-12-11 | 2019-05-14 | 中煤科工集团西安研究院有限公司 | 单发多收观测装置瞬变电场数据动校正方法及装置 |
CN109915120A (zh) * | 2019-04-16 | 2019-06-21 | 北京恒泰万博石油技术股份有限公司 | 一种基于环境因素的随钻电阻率测井系统校正方法 |
CN110967773A (zh) * | 2019-12-23 | 2020-04-07 | 中国煤炭地质总局地球物理勘探研究院 | 煤层内部富水性计算方法、装置及电子设备 |
CN112268913A (zh) * | 2020-09-18 | 2021-01-26 | 天津大学 | 消除水的矿化度影响的油气水三相流微波持水率测量方法 |
AU2021100785A4 (en) * | 2021-02-09 | 2021-05-13 | The Third Geological Exploration Institute Of Qinghai Province | Method and apparatus for determining the time window for a ground-airborne transient electromagnetic survey |
CN115097534A (zh) * | 2022-07-06 | 2022-09-23 | 中国地质调查局水文地质环境地质调查中心 | 干热岩压裂改造空间形态识别方法和系统 |
CN115980861A (zh) * | 2021-10-14 | 2023-04-18 | 中国石油化工股份有限公司 | 用于随钻方位电磁波探测界面的反演方法及其实现系统 |
US20230375739A1 (en) * | 2020-09-09 | 2023-11-23 | Schlumberger Technology Corporation | Methods and apparatus for determining dielectric constant and resistivity with electromagnetic propagation measurements |
WO2024030779A1 (en) * | 2022-08-01 | 2024-02-08 | Schlumberger Technology Corporation | Direct determination of formation apparent conductivity from em propagation measurements |
-
2024
- 2024-04-01 CN CN202410385363.5A patent/CN117967288A/zh active Pending
Patent Citations (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4209747A (en) * | 1977-09-21 | 1980-06-24 | Schlumberger Technology Corporation | Apparatus and method for determination of subsurface permittivity and conductivity |
US5345179A (en) * | 1992-03-09 | 1994-09-06 | Schlumberger Technology Corporation | Logging earth formations with electromagnetic energy to determine conductivity and permittivity |
CA2375432A1 (en) * | 1999-06-21 | 2000-12-28 | Teruhiko Hagiwara | Method for radial profiling of resistivity at multiple depths of investigation |
CN101050700A (zh) * | 2007-05-09 | 2007-10-10 | 中国海洋石油总公司 | 一种用高导磁材料的电磁波测井仪校验方法 |
CN101609169A (zh) * | 2009-04-08 | 2009-12-23 | 中国石油集团钻井工程技术研究院 | 一种提高电磁波电阻率测量精度和扩展其测量范围的方法 |
CN102606138A (zh) * | 2012-03-31 | 2012-07-25 | 中国电子科技集团公司第二十二研究所 | 随钻电磁波电阻率测井仪相位幅度介电常数校正方法 |
CN103015970A (zh) * | 2012-11-30 | 2013-04-03 | 中国海洋石油总公司 | 一种随钻电阻率测井仪的模拟检测设备 |
CN103645512A (zh) * | 2013-12-10 | 2014-03-19 | 江苏大学 | 一种多层积水采空区探测方法 |
CN104912544A (zh) * | 2015-05-29 | 2015-09-16 | 中国航天科技集团公司烽火机械厂 | 一种油气井井下温度压力参数探测系统及方法 |
CN105700026A (zh) * | 2016-01-26 | 2016-06-22 | 中国科学院地质与地球物理研究所 | 一种地下目标体导电指数提取方法与装置 |
CN106054257A (zh) * | 2016-05-11 | 2016-10-26 | 中国科学院地质与地球物理研究所 | 一种含磁层状介质电磁响应的确定系统 |
CN107605475A (zh) * | 2017-10-27 | 2018-01-19 | 罗淮东 | 用于地层测试的设备、系统及方法 |
CN109061751A (zh) * | 2018-06-20 | 2018-12-21 | 西安石油大学 | 多层介质的视电导率计算方法 |
CN109252841A (zh) * | 2018-11-16 | 2019-01-22 | 中国海洋石油集团有限公司 | 基于电磁波无线监测提高注聚井注入能力的系统和方法 |
CN109752762A (zh) * | 2018-12-11 | 2019-05-14 | 中煤科工集团西安研究院有限公司 | 单发多收观测装置瞬变电场数据动校正方法及装置 |
CN109915120A (zh) * | 2019-04-16 | 2019-06-21 | 北京恒泰万博石油技术股份有限公司 | 一种基于环境因素的随钻电阻率测井系统校正方法 |
CN110967773A (zh) * | 2019-12-23 | 2020-04-07 | 中国煤炭地质总局地球物理勘探研究院 | 煤层内部富水性计算方法、装置及电子设备 |
US20230375739A1 (en) * | 2020-09-09 | 2023-11-23 | Schlumberger Technology Corporation | Methods and apparatus for determining dielectric constant and resistivity with electromagnetic propagation measurements |
CN112268913A (zh) * | 2020-09-18 | 2021-01-26 | 天津大学 | 消除水的矿化度影响的油气水三相流微波持水率测量方法 |
AU2021100785A4 (en) * | 2021-02-09 | 2021-05-13 | The Third Geological Exploration Institute Of Qinghai Province | Method and apparatus for determining the time window for a ground-airborne transient electromagnetic survey |
CN115980861A (zh) * | 2021-10-14 | 2023-04-18 | 中国石油化工股份有限公司 | 用于随钻方位电磁波探测界面的反演方法及其实现系统 |
CN115097534A (zh) * | 2022-07-06 | 2022-09-23 | 中国地质调查局水文地质环境地质调查中心 | 干热岩压裂改造空间形态识别方法和系统 |
WO2024030779A1 (en) * | 2022-08-01 | 2024-02-08 | Schlumberger Technology Corporation | Direct determination of formation apparent conductivity from em propagation measurements |
Non-Patent Citations (5)
Title |
---|
刘洋: "三维井—地电位有限差分数值模拟正演研究", 《中国优秀硕士论文全文库工程科技I辑》, no. 4, 15 April 2018 (2018-04-15) * |
张长雷: "探地雷达技术在复杂地质条件下的应用研究", 《中国优秀硕士论文全文库工程科技II辑》, no. 4, 15 April 2013 (2013-04-15) * |
石亚兰: "浅海重磁电联合勘探方法试验及效果", 《工程地球物理学报》, vol. 13, no. 5, 1 July 2016 (2016-07-01) * |
赵云生: "井地差分电磁法的复电阻率研究及其在储层含油饱和度评价中的应用", 《中国优秀硕士论文全文库工程科技I辑》, no. 3, 15 March 2022 (2022-03-15) * |
陶果: "对感应测井曲线进行分层、环境校正的方法", 《华东石油学院学报(自然科学版)》, no. 2, 28 February 1987 (1987-02-28) * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Diaz et al. | Drilling data from an enhanced geothermal project and its pre-processing for ROP forecasting improvement | |
US6810331B2 (en) | Fixed-depth of investigation log for multi-spacing multi-frequency LWD resistivity tools | |
EP2810101B1 (en) | Improving efficiency of pixel-based inversion algorithms | |
US6182013B1 (en) | Methods and apparatus for dynamically estimating the location of an oil-water interface in a petroleum reservoir | |
BR112015009197B1 (pt) | Método e sistema para realizar uma operação de perfuração | |
CN107120111B (zh) | 一种基于多重分形的油藏井间连通程度的评估方法及系统 | |
US9341739B2 (en) | Apparatus and method for estimating geologic boundaries | |
BRPI1001536A2 (pt) | mÉtodo para determinar permeabilidades eficazes de formaÇÕes terrestres | |
MX2013012178A (es) | Metodos y sistemas para estimar la resistividad y porosidad de la formacion. | |
US20130176822A1 (en) | Method for estimating elastic parameters by inverting 4d seismic measurements | |
US10495524B2 (en) | Apparatus and method for monitoring production wells | |
WO2004095077A1 (en) | Method for predicting pore pressure | |
CN115469366B (zh) | 一种油基泥浆微电阻率扫描成像测井参数确定方法 | |
Garg et al. | Appropriate use of USGS volumetric “heat in place” method and Monte Carlo calculations | |
CN113176618A (zh) | 一种深部地温场预测方法、设备及存储介质 | |
EP2957716A2 (en) | Determining a quantitative bond using signal attenuation | |
WO2020209867A1 (en) | Machine learning for evaluating and categorizing inversion solutions | |
CN113486539B (zh) | 一种深水钻井全过程中地层压力不确定性分析方法 | |
CN117967288A (zh) | 一种用于油气田的井下压力监测系统和监测方法 | |
EP1558837A1 (en) | Method and device for determining the resistivity in a geological formation crossed by a cased well | |
CN112906465B (zh) | 一种基于地层因素的煤系地层声波曲线重构方法及系统 | |
CN112268923B (zh) | 一种基于测井曲线获取地层热导率的方法 | |
CN112241025B (zh) | 一种井震联合地层压力确定方法及系统 | |
NO20171368A1 (en) | Drill bit positioning system | |
CN103485772B (zh) | 测井设备、方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |