CN109711109A - 一种电阻率测井仪器电极系结构参数智能优化的方法及装置 - Google Patents
一种电阻率测井仪器电极系结构参数智能优化的方法及装置 Download PDFInfo
- Publication number
- CN109711109A CN109711109A CN201910132183.5A CN201910132183A CN109711109A CN 109711109 A CN109711109 A CN 109711109A CN 201910132183 A CN201910132183 A CN 201910132183A CN 109711109 A CN109711109 A CN 109711109A
- Authority
- CN
- China
- Prior art keywords
- objective function
- function
- overall goal
- optimization
- weight
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 125
- 238000000034 method Methods 0.000 title claims abstract description 28
- 230000000694 effects Effects 0.000 claims abstract description 76
- 238000010276 construction Methods 0.000 claims abstract description 23
- 238000001514 detection method Methods 0.000 claims description 46
- 238000011835 investigation Methods 0.000 claims description 28
- 238000012937 correction Methods 0.000 claims description 27
- 230000015572 biosynthetic process Effects 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 18
- 238000010606 normalization Methods 0.000 claims description 15
- 239000011435 rock Substances 0.000 claims description 14
- 238000012546 transfer Methods 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 8
- 230000005611 electricity Effects 0.000 claims description 7
- 238000007781 pre-processing Methods 0.000 claims description 6
- 238000003786 synthesis reaction Methods 0.000 claims description 6
- 230000009545 invasion Effects 0.000 claims description 5
- 238000009413 insulation Methods 0.000 claims description 4
- 238000005303 weighing Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 4
- 206010021703 Indifference Diseases 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 230000000873 masking effect Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000002922 simulated annealing Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
本发明涉及电阻率测井领域,公开了一种电阻率测井仪器电极系结构参数智能优化的方法及装置。方法包括获取输入的影响电阻率测井仪器电极系结构参数的多个因素;根据多个因素对应的多个目标函数以预先设定的权重求和获得可通过算法构造并计算多个由电极系结构参数组构成的仪器模型的总体目标函数;根据算法构造并计算总体目标函数,对多个目标函数的权重进行自适应调整以使多个目标函数的扰动对总体目标函数的影响相等;输出优化结果。本发明可以避免某些目标函数值过大掩盖了某些目标函数值很小但对电阻率测井仪器的电极系结构参数优化影响同样重要的因素导致优化的效果片面,从而获得对电阻率测井仪器电极系结构参数影响更全面的优化效果。
Description
技术领域
本发明涉及电阻率测井领域,具体的是涉及一种电阻率测井仪器电极系结构参数智能优化的方法及装置。
背景技术
在低阻泥浆条件下,电阻率测井是探测油气藏储层电阻率的重要方法,如双侧向电阻率测井仪器、阵列侧向电阻率测井仪器、阵列方位电阻率测井仪器等。电阻率测井仪器各电极结构及间隔尺寸等电极系结构参数决定了仪器的探测特性,好的电极系结构一方面能够使仪器测量值受井眼、围岩、侵入等因素影响小,测量值更接近地层电阻率真值;另一方面,仪器能提供多种探测模式,不同探测模式具有径向分布均匀的探测深度,能更清晰反映近井眼地层侵入情况;要得到较好的仪器特性,就必须研究电极系结构对仪器特性的影响,对电极系结构参数进行优化。
然而,在现有技术对电阻率测井仪器电极系结构参数的优化中,根据影响电阻率测井仪器特性的主要因素建立多个目标函数,利用算法进行优化时,未采取手段提高其计算的速度、扩大其计算规模,使得计算效率不高;也没有考虑到某些目标函数值过大,会掩盖某些目标函数值很小但是对电阻率测井仪器的电极系结构参数影响同样重要的因素,导致某些目标函数值很小的因素的影响无法与其重要性匹配地体现在优化结果中,使得优化结果片面化;并且现有技术中没有对仪器的主电极-屏蔽电极电流比进行控制,使得优化后的仪器由于主电极-屏蔽电极电流比过低,导致仪器无法在高阻地层条件下工作。
发明内容
(一)发明目的
为克服上述现有技术中存在的至少一种缺陷,本发明提供一种电阻率测井仪器电极系结构参数智能优化的方法及装置,该电阻率测井仪器电极系结构参数智能优化的方法可以在利用算法进行优化时,提高其计算速度并扩大计算规模,使得计算效率大大提高,还可以避免某些目标函数值过大掩盖了某些目标函数值很小但对电阻率测井仪器的电极系结构参数影响同样重要的因素导致优化的效果片面,从而获得对电阻率测井仪器电极系结构参数影响更全面的优化效果,并且对仪器的主电极-屏蔽电极电流比进行控制,使得仪器在高阻地层条件下也可以工作。
(二)技术方案
作为本发明的第一方面,本发明公开了一种电阻率测井仪器电极系结构参数智能优化的方法,包括:
在因素获取模块中,获取输入的影响电阻率测井仪器电极系结构参数的多个因素;
在预处理模块中,根据所述多个因素对应的多个目标函数以预先设定的权重求和获得可通过算法构造并计算多个由电极系结构参数组构成的仪器模型的总体目标函数;
在权重自适应调整模块中,根据算法构造并计算所述多个仪器模型的总体目标函数,并判断所述总体目标函数是否满足所述多个目标函数对其的影响相等的第二优化终止条件,如果满足,优化结束,输出优化结果,否则,基于所述多个目标函数的扰动对所述总体目标函数的影响相等的原则执行对所述多个目标函数的权重进行自适应调整,权重调整后,将执行按新的权重返回继续计算所述总体目标函数;根据算法构造并计算多个仪器模型的总体目标函数时采用迭代计算;每重复计算一次即为一次迭代,每次迭代得到的结果作为下一次迭代的初始值;根据算法生成多组用于构造仪器模型的电阻率测井仪器电极系结构参数后,可采用并行正演计算所述多组仪器模型对应的总体目标函数,对不收敛的总体目标函数计算其收敛速度,所述并行正演计算可以一次计算多个仪器模型的总体目标函数,达到提高计算速度并扩大计算规模,使得计算效率大大提高;
在优化结果输出模块中,输出满足所述优化终止条件的仪器电极系结构参数。
在一种可能的实施方式中,所述在权重自适应调整模块中的第二优化终止条件之前,还包括第一优化终止条件:
所述第一优化终止条件为当根据算法计算出的所述总体目标函数收敛时,优化结束,当其不收敛时,执行计算其收敛速度;
当计算出的所述总体目标函数的收敛速度与预先设定的总体目标函数收敛速度的最小值相比较数值较大时,则执行返回继续计算所述总体目标函数;
所述第二优化终止条件为当计算出的所述总体目标函数的收敛速度与预先设定的总体目标函数收敛速度的最小值相比较数值较小时,计算所述各个目标函数的扰动引起所述总体目标函数的变化贡献是否相等,若相等,优化结束;
当计算出的所述多个目标函数的扰动引起所述总体目标函数的变化贡献不相等时,则基于所述多个目标函数的扰动引起所述总体目标函数的变化贡献相等的原则执行对所述多个目标函数的权重进行自适应调整,权重调整后,将执行按新的权重返回继续计算所述总体目标函数。
在一种可能的实施方式中,所述计算多个目标函数的扰动引起所述总体目标函数的变化贡献,包括:
基于数据误差传递理论,所述多个目标函数的扰动体现为所述多个目标函数变量的扰动,以及所述多个目标函数变量的扰动相等,建立所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式;
所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式为:其中,Fi为根据所述多个目标函数中的变量u1,u2,...,un,由所述多个目标函数按F=F(u1,u2,...,un)合成的所述总体目标函数,n为变量的数量,i为整个优化过程的迭代次数;
根据所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式计算得出所述多个目标函数的扰动引起所述总体目标函数的变化贡献为 所述基于数据误差传递理论,即为基于测量结果的合成标准不确定度由各基本测量数据的误差传递而来。
在一种可能的实施方式中,所述基于所述多个目标函数的扰动引起所述总体目标函数的变化贡献相等的原则执行对所述多个目标函数的权重进行自适应调整,包括:
根据所述多个目标函数中每个目标函数的变量比例调整所述多个目标函数的权重比例,以使调整后所述多个目标函数中每个目标函数的变量比例满足将所述比例中的最小比例归一化,若最大比例大于Wmin,将所述多个目标函数的权重调整为最小比例归一化后的对应值,调整成功,否则调整不成功,其中Wmin为预先设定的所述多个目标函数权重调整的最小比例因子,其值大于1。
在一种可能的实施方式中,所述多个目标函数包括第一目标函数,所述第一目标函数包括:控制仪器探测深度的伪几何因子目标函数f1、控制仪器井眼影响的影响因子目标函数f2、控制仪器围岩影响的影响因子目标函数f3和控制仪器主电极-屏蔽电极电流比的电极电流目标函数f4中的任意至少一项;所述第一目标函数的计算方法为:
所述伪几何因子目标函数f1为不同探测模式的探测深度目标函数以预先设定的权重求和,其公式为:式中W1,i为第i个探测模式目标函数的权重,f1,i为第i个探测模式探测深度目标函数,目标函数f1,i的公式为式中ri为所优化的电阻率测井仪器第i种探测模式期望的探测深度,其中仪器探测深度定义为伪几何因子为0.5时的侵入半径r,式中Rt为伪几何因子地层模型原状地层电阻率,Ra(r)为伪几何因子地层模型在侵入半径r下的电极系对地层模型响应的视电阻率,Rxo为伪几何因子地层模型侵入带电阻率;
所述影响因子目标函数f2为不同探测模式的井眼影响目标函数以预先设定的权重求和,其公式为:式中W2,i,j,k为第i个探测模式、第j个控制点和第k个井眼半径目标函数的权重,f2,i,j,k为第i个探测模式、第j个控制点xi,j第k个井眼半径下的井眼影响目标函数,目标函数f2,i,j,k的公式为Rt=xi,j,wr=wrk,式中ti,j,k为不同控制点xi,j下不同井眼半径的井眼校正系数的控制目标,通常接近1,式中Rt/Ra为井眼影响地层模型井眼校正系数,Rt为井眼影响地层模型原状地层电阻率,Ra为井眼影响地层模型利用电阻率为Rt的均匀地层计算的视电阻率,wr为井眼影响地层模型井眼半径;其中仪器井眼影响通过不同Ra/Rm情况下的井眼校正系数Rt/Ra衡量,Rm为井眼影响地层模型井眼电阻率;理想情况下井眼校正系数为1;
所述影响因子目标函数f3为不同探测模式的围岩影响目标函数以预先设定的权重求和,其公式为:式中W3,i,j,k为第i个探测模式、第j个控制点和第k个围岩/目的层对比度下目标函数的权重,f3,i,j,k为第i个探测模式、第j个控制点hi,j和第k个围岩/目的层对比度下围岩影响目标函数,目标函数f3,i,j,k的公式为式中ti,j,k为不同控制点hi,j下不同围岩/目的层对比度下围岩校正系数的控制目标,通常接近1,式中Rt为围岩影响地层模型目的层电阻率,Rs为围岩影响地层模型围岩电阻率,Rt/Ra为围岩影响地层模型围岩校正系数,Rt/Rs为围岩影响地层模型地层对比度;其中围岩影响通过不同目的层厚度下的围岩校正系数Rt/Ra衡量,理想情况下围岩校正系数为1;
所述电极电流目标函数f4的公式为设主电极电流为单位电流,式中Ii为以仪器测量最高电阻率为均匀地层、最深探测模式第i个屏蔽电极的电流,threshold为设定的主电极-屏蔽电极电流比的最小值;其中仪器主电极-屏蔽电极电流比反映了仪器聚焦情况下所需屏蔽电流大小,其比值越小,电路硬件聚焦越容易实现;设地层为均匀地层,地层真电阻率为仪器工作的最高电阻率地层,在主电极发射1A电流情况下,计算此时仪器最远探测模式下各电极电流,进而根据电极电流目标函数f4的定义计算该目标函数;根据仪器探测的基本性质,仪器工作的目的层电阻率越高,探测深度越深,达到聚焦条件所需的屏蔽电流越大,因此在计算电极电流目标函数f4时,仅考虑仪器测量的最高电阻率和最深探测模式,使得优化后的仪器也可以在高阻地层条件下工作;
在一种可能的实施方式中,所述总体目标函数还包括控制电阻率测井仪器特征的第二目标函数,所述仪器特征包括仪器总长度、仪器主电极长度和仪器电极之间的最小绝缘长度中的任意至少一项。
作为本发明的第二方面,本发明公开了一种电阻率测井仪器电极系结构参数智能优化装置,包括:
因素获取模块,用于获取输入的影响电阻率测井仪器电极系结构参数的多个因素;
预处理模块,用于所述多个因素对应的多个目标函数以预先设定的权重求和获得可通过算法构造并计算多个由电极系结构参数组构成的仪器模型的总体目标函数;
权重自适应调整模块,用于根据算法构造并计算所述多个仪器模型的总体目标函数,并判断所述总体目标函数是否满足所述多个目标函数的扰动对其的影响相等的第二优化终止条件,如果满足,优化结束,输出优化结果,否则,基于所述多个目标函数的扰动对所述总体目标函数的影响相等的原则执行对所述多个目标函数的权重进行自适应调整,权重调整后,将执行按新的权重返回继续计算所述总体目标函数;根据算法构造并计算多个仪器模型的总体目标函数时采用迭代计算;每重复计算一次即为一次迭代,每次迭代得到的结果作为下一次迭代的初始值;根据算法生成多组用于构造仪器模型的电阻率测井仪器电极系结构参数后,可采用并行正演计算所述多组仪器模型对应的总体目标函数,对不收敛的总体目标函数计算其收敛速度,所述并行正演计算可以一次计算多个仪器模型的总体目标函数,达到提高计算速度并扩大计算规模,使得计算效率大大提高;
优化结果输出模块,用于输出满足所述优化终止条件的仪器电极系结构参数。
在一种可能的实施方式中,所述权重自适应调整模块包括含有所述第二优化终止条件的第二优化终止单元,还包括含有所述第一优化终止条件的第一优化终止单元:
所述第一优化终止单元,包括第一优化终止条件,用于当根据算法计算出的所述总体目标函数收敛时,优化结束,当其不收敛时,执行计算其收敛速度;
当计算出的所述总体目标函数的收敛速度与预先设定的总体目标函数收敛速度的最小值相比较数值较大时,则执行返回继续计算所述总体目标函数;
所述第二优化终止单元,包括第二优化终止条件,用于当计算出的所述总体目标函数的收敛速度与预先设定的总体目标函数收敛速度的最小值相比较数值较小时,计算所述多个目标函数的扰动引起所述总体目标函数的变化贡献是否相等,若相等,优化结束;
当计算出的所述多个目标函数的扰动引起所述总体目标函数的变化贡献不相等时,则基于所述多个目标函数的扰动引起所述总体目标函数的变化贡献相等的原则执行对所述多个目标函数的权重进行自适应调整,权重调整后,将执行按新的权重返回继续计算所述总体目标函数。
在一种可能的实施方式中,所述第二优化终止单元中的计算多个目标函数的扰动引起总体目标函数的变化贡献,包括:
第一计算单元,用于基于数据误差传递理论,所述多个目标函数的扰动体现为所述多个目标函数变量的扰动,以及所述多个目标函数变量的扰动相等,建立所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式;
所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式为:其中,Fi为根据所述多个目标函数中的变量u1,u2,...,un,由所述多个目标函数按F=F(u1,u2,...,un)合成的所述总体目标函数,n为变量的数量,i为整个优化过程的迭代次数;
根据所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式计算得出所述多个目标函数的扰动引起所述总体目标函数的变化贡献为 所述基于数据误差传递理论,即为基于测量结果的合成标准不确定度由各基本测量数据的误差传递而来。
在一种可能的实施方式中,所述第二优化终止单元中的基于所述多个目标函数的扰动引起所述总体目标函数的变化贡献相等的原则自适应调整所述多个目标函数的权重,包括:
第二计算单元,用于根据所述多个目标函数中每个目标函数的变量比例调整所述多个目标函数的权重比例,以使调整后所述多个目标函数中每个目标函数的变量比例满足将所述比例中的最小比例归一化,若最大比例大于Wmin,将所述多个目标函数的权重调整为最小比例归一化后的对应值,调整成功,否则调整不成功,其中Wmin为预先设定的所述多个目标函数权重调整的最小比例因子,其值大于1。
(三)有益效果
本发明提供的一种电阻率测井仪器电极系结构参数智能优化的方法及装置,具有如下有益效果:
通过获取输入的影响电阻率测井仪器电极系结构参数的多个因素,根据多个因素对应的多个目标函数以预先设定的权重求和获得可通过算法构造并计算多个由电极系结构参数组构成的仪器模型的总体目标函数;根据算法构造并迭代计算多个仪器模型的总体目标函数,并在根据算法生成多组用于构造仪器模型的电阻率测井仪器电极系结构参数后,采用并行正演计算多组仪器模型对应的总体目标函数,可以提高其计算速度并扩大计算规模,使得计算效率大大提高;而基于数据误差传递理论计算多个目标函数变量的扰动引起总体目标函数的变化贡献,自适应调整多个目标函数在总体目标函数中的权重以使多个目标函数的扰动对总体目标函数的影响相等,可以避免某些目标函数值过大掩盖了某些目标函数值很小但对电阻率测井仪器的电极系结构参数影响同样重要的因素导致优化的效果片面,从而获得对电阻率测井仪器电极系结构参数影响更全面的优化效果;对仪器的主电极-屏蔽电极电流比进行控制,使得优化后的仪器不会出现由于主电极-屏蔽电极电流比过低,而导致仪器无法在高阻地层条件下工作的情况。
附图说明
以下参考附图描述的实施例是示例性的,旨在用于解释和说明本发明,而不能理解为对本发明的保护范围的限制。
图1为本发明提供的第一实施例的电阻率测井仪器电极系结构参数智能优化的方法的流程图。
图2为本发明提供的第一实施例的较佳的实施流程图。
图3为本发明提供的伪几何因子地层模型示意图。
图4为本发明提供的井眼影响地层模型示意图。
图5为本发明提供的围岩影响地层模型示意图。
图6为本发明提供的第二实施例电阻率测井仪器电极系结构参数智能优化装置的结构示意图。
附图标记:11-伪几何因子地层模型原状地层电阻率Rt,12-伪几何因子地层模型侵入带电阻率Rxo,13-伪几何因子地层模型井眼电阻率Rm,21-井眼影响地层模型原状地层电阻率Rt,23-井眼影响地层模型井眼电阻率Rm,31-围岩影响地层模型目的层电阻率Rt,33-围岩影响地层模型井眼电阻率Rm,34-围岩影响地层模型围岩电阻率Rs。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚,下面结合附图和实施例对本发明进行更加详细的描述。
需要说明的是:所描述的实施例是本发明一部分实施例,而不是全部的实施例,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本文中,“第一”、“第二”、“f1”、“f2”“f3”、“f4”等仅用于彼此的区分,而非表示它们的重要程度及顺序等。
本文中的模块、单元的划分仅仅是一种逻辑功能的划分,在实际实现时可以有其他的划分方式,例如多个模块和/或单元可以结合或集成于另一个系统中。作为分离部件说明的模块、单元在物理上可以是分开的,也可以是不分开的。因此可以根据实际需要选择其中的部分或全部的单元来实现实施例的方案。
下面参考图1-图5详细描述本发明提供的一种电阻率测井仪器电极系结构参数智能优化的方法的第一实施例:本实施例主要应用于电阻率测井仪器。
如图所示,本实施例提供的电阻率测井仪器电极系结构参数智能优化的方法,包括:
101,在因素获取模块中,获取输入的影响电阻率测井仪器电极系结构参数的多个因素;所述因素包括仪器探测深度、仪器井眼影响、仪器围岩影响和仪器主电极-屏蔽电极电流比中的任意至少一个;所述因素还可以包括仪器总长度、仪器主电极长度和仪器电极之间的最小绝缘长度中的任意至少一项;
102,在预处理模块中,根据所述多个因素对应的多个目标函数以预先设定的权重求和获得可通过算法构造并计算多个由电极系结构参数组构成的仪器模型的总体目标函数;所述多个目标函数包括第一目标函数,还可以包括第二目标函数,所述多个目标函数以预先设定的权重求和获得总体目标函数,可以预先设定多个目标函数的初始权重为1;所述总体目标函数的公式为:式中Wi为第i个目标函数的权重,fi为第i个目标函数;
103,在权重自适应调整模块中,根据算法构造并计算所述多个仪器模型的总体目标函数,并判断所述总体目标函数是否满足所述多个目标函数的扰动对其的影响相等的第二优化终止条件,如果满足,优化结束,输出优化结果,否则,自适应调整所述多个目标函数在所述总体目标函数中的权重以使所述多个目标函数的扰动对所述总体目标函数的影响相等,避免某些目标函数值过大掩盖了某些目标函数值很小但对电阻率测井仪器的特性影响同样重要的因素导致优化的效果片面,从而获得对电阻率测井仪器特性影响更全面的优化效果;根据算法构造并计算多个仪器模型的总体目标函数时采用迭代计算;每重复计算一次即为一次迭代,每次迭代得到的结果作为下一次迭代的初始值;根据算法生成多组用于构造仪器模型的电阻率测井仪器电极系结构参数后,可采用并行正演计算所述多组仪器模型对应的总体目标函数,对不收敛的总体目标函数计算其收敛速度,所述并行正演计算可以一次计算多个仪器模型的总体目标函数,达到提高计算速度并扩大计算规模,使得计算效率大大提高;
所述算法包括遗传算法、人工神经网络、蚁群算法、粒子群算法、模拟退火算法、差分进化算法中的一种或者多种;
104,在优化结果输出模块中,输出满足所述优化终止条件的仪器电极系结构参数。
在一种可能的实施方式中,所述在权重自适应调整模块中的第二优化终止条件之前,还包括第一优化终止条件:
所述第一优化终止条件为当根据算法计算出的所述总体目标函数收敛时,优化结束,当其不收敛时,执行计算其收敛速度;
当计算出的所述总体目标函数的收敛速度与预先设定的总体目标函数收敛速度的最小值相比较数值较大时,则执行返回继续计算所述总体目标函数;
所述总体目标函数的收敛速度的公式为式中Fi表示第i次优化的总体目标函数;vmin为预先定义的总体目标函数收敛速度最小值,其值大于1;
所述第二优化终止条件为当计算出的所述总体目标函数的收敛速度与预先设定的总体目标函数收敛速度的最小值相比较数值较小时,计算所述多个目标函数的扰动引起所述总体目标函数的变化贡献是否相等,若相等,优化结束;
当计算出的所述多个目标函数的扰动引起所述总体目标函数的变化贡献不相等时,则基于所述多个目标函数的扰动引起所述总体目标函数的变化贡献相等的原则执行对所述多个目标函数的权重进行自适应调整,权重调整后,将执行按新的权重返回继续计算所述总体目标函数。
在一种可能的实施方式中,所述计算多个目标函数的扰动引起所述总体目标函数的变化贡献,包括:
基于数据误差传递理论,所述多个目标函数的扰动体现为所述多个目标函数变量的扰动,以及所述多个目标函数变量的扰动相等,建立所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式;
所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式为:其中,Fi为根据所述多个目标函数中的变量u1,u2,...,un,由所述多个目标函数按F=F(u1,u2,...,un)合成的所述总体目标函数,n为变量的数量,i为整个优化过程的迭代次数;所述每个目标函数包括至少一个变量,所述每个目标函数的权重为其所有变量对应值之和;
根据所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式计算得出所述多个目标函数的扰动引起所述总体目标函数的变化贡献为 所述基于数据误差传递理论,即为基于测量结果的合成标准不确定度由各基本测量数据的误差传递而来。
在一种可能的实施方式中,所述基于所述多个目标函数的扰动引起所述总体目标函数的变化贡献相等的原则执行对所述多个目标函数的权重进行自适应调整,包括:
根据所述多个目标函数中每个目标函数的变量比例调整所述多个目标函数的权重比例,以使调整后所述多个目标函数中每个目标函数的变量比例满足将所述比例中的最小比例归一化,若最大比例大于Wmin,将所述多个目标函数的权重调整为最小比例归一化后的对应值,调整成功,否则调整不成功,其中Wmin为预先设定的所述多个目标函数权重调整的最小比例因子,其值大于1。
表示所述总体目标函数对第j个目标函数变量的偏导在第i次优化计算的目标函数fj(i)的取值,这个值决定了当前状态下第j个目标函数对总体目标函数的贡献,即第j个目标函数变量的扰动对总体目标函数的影响,将目标函数权重比例调整为即实现了所述多个目标函数变量对所述总体目标函数的贡献相等,使得优化过程中各个变化因素对所述总体目标函数的影响无差别,避免了某些目标函数值过大掩盖了某些目标函数值很小但对电阻率测井仪器的特性影响同样重要的因素的情况,达到了更全面的优化效果。
在一种可能的实施方式中,所述多个目标函数包括第一目标函数,所述第一目标函数包括:控制仪器探测深度的伪几何因子目标函数f1、控制仪器井眼影响的影响因子目标函数f2、控制仪器围岩影响的影响因子目标函数f3和控制仪器主电极-屏蔽电极电流比的电极电流目标函数f4中的任意至少一项;所述第一目标函数的计算方法为:
所述伪几何因子目标函数f1为不同探测模式的探测深度目标函数以预先设定的权重求和,其公式为:式中W1,i为第i个探测模式目标函数的权重,f1,i为第i个探测模式探测深度目标函数,目标函数f1,i的公式为式中ri为所优化的电阻率测井仪器第i种探测模式期望的探测深度,其中仪器探测深度定义为伪几何因子为0.5时的侵入半径r,式中Rt为伪几何因子地层模型原状地层电阻率11,Ra(r)为伪几何因子地层模型在侵入半径r下的电极系对地层模型响应的视电阻率,Rxo为伪几何因子地层模型侵入带电阻率12;根据图3所示的伪几何因子地层模型,设侵入层半径为r,计算伪几何因子G,包括:
根据图3所示地层模型,计算侵入层半径为r时的视电阻率,记为Ra(r);
在图3所示地层模型中,Rm为伪几何因子地层模型井眼电阻率13,设地层为均匀地层,Rxo为伪几何因子地层模型侵入带电阻率12,计算此时的视电阻率替代原始给定的Rxo;
在图3所示地层模型中,设地层为均匀地层,Rt为伪几何因子地层模型原状地层电阻率11,计算此时的视电阻率替代原始给定的Rt;
根据公式计算侵入层半径为r时的伪几何因子;
根据仪器探测的基本性质,随侵入半径增大,仪器的伪几何因子递增,反过来说,如果某个探测模式在侵入半径ri下G(ri)<0.5,则仪器该模式的探测深度即大于ri,达到了优化目标。因此,建立伪几何因子目标函数仅考虑G(ri)>0.5;
所述影响因子目标函数f2为不同探测模式的井眼影响目标函数以预先设定的权重求和,其公式为:式中W2,i,j,k为第i个探测模式、第j个控制点和第k个井眼半径目标函数的权重,f2,i,j,k为第i个探测模式、第j个控制点xi,j第k个井眼半径下的井眼影响目标函数,目标函数f2,i,j,k的公式为Rt=xi,j,wr=wrk,式中ti,j,k为不同控制点xi,j下不同井眼半径的井眼校正系数的控制目标,通常接近1,式中Rt/Ra为井眼影响地层模型井眼校正系数,Rt为井眼影响地层模型原状地层电阻率21,Ra为井眼影响地层模型利用电阻率为Rt的均匀地层计算的视电阻率,wr为井眼影响地层模型井眼半径;其中仪器井眼影响通过不同Ra/Rm情况下的井眼校正系数Rt/Ra衡量,Rm为井眼影响地层模型井眼电阻率23;理想情况下井眼校正系数为1;根据图4所示的地层模型,计算井眼影响,包括:
根据图4所示地层模型及设定的井眼半径和地层真电阻率,计算视电阻率,记为Ra;
计算井眼校正系数Rt/Ra;
根据仪器探测的基本性质,仪器的井眼影响特性随井眼半径增大而单调变化,因此,在建立影响因子目标函数f2时,仅需要考虑控制仪器工作的最大井眼半径和最小井眼半径。此外,由于井眼校正系数越小越好,因此当满足时,使目标函数f2,i,j,k为0;
所述影响因子目标函数f3为不同探测模式的围岩影响目标函数以预先设定的权重求和,其公式为:式中W3,i,j,k为第i个探测模式、第j个控制点和第k个围岩/目的层对比度下目标函数的权重,f3,i,j,k为第i个探测模式、第j个控制点hi,j和第k个围岩/目的层对比度下围岩影响目标函数,目标函数f3,i,j,k的公式为式中ti,j,k为不同控制点hi,j下不同围岩/目的层对比度下围岩校正系数的控制目标,通常接近1,式中Rt为围岩影响地层模型目的层电阻率31,Rs为围岩影响地层模型围岩电阻率34,Rt/Ra为围岩影响地层模型围岩校正系数,Rt/Rs为围岩影响地层模型地层对比度;其中围岩影响通过不同目的层厚度下的围岩校正系数Rt/Ra衡量,理想情况下围岩校正系数为1;根据图5所示的地层模型,Rm为围岩影响地层模型井眼电阻率33,计算围岩影响,包括:
根据图5所示地层模型及设定的目的层厚度,设定不同的Rt,Rs,使得目的层与围岩对比度分别达到最大与最小,计算各自的视电阻率,记为Ra;
计算井眼围岩系数Rt/Ra;
根据仪器探测的基本性质,仪器的围岩影响特性随地层对比度Rt/Rs变化而单调变化,因此,在建立影响因子目标函数f3时,仅需要考虑控制仪器工作的最大地层对比度和最小地层对比度。此外,由于围岩校正系数越小越好,因此当满足时,使目标函数f3,i,j,k为0;
所述电极电流目标函数f4的公式为设主电极电流为单位电流,式中Ii为以仪器测量最高电阻率为均匀地层、最深探测模式第i个屏蔽电极的电流,threshold为设定的主电极-屏蔽电极电流比的最小值;其中仪器主电极-屏蔽电极电流比反映了仪器聚焦情况下所需屏蔽电流大小,其比值越小,电路硬件聚焦越容易实现;设置地层为均匀地层,地层真电阻率为仪器工作的最高电阻率地层,在主电极发射1A电流情况下,计算此时仪器最远探测模式下各电极电流,进而根据电极电流目标函数f4的定义计算该目标函数;根据仪器探测的基本性质,仪器工作的目的层电阻率越高,探测深度越深,达到聚焦条件所需的屏蔽电流越大,因此在计算电极电流目标函数f4时,仅考虑仪器测量的最高电阻率和最深探测模式,使得优化后的仪器不会出现由于主电极-屏蔽电极电流比过低,而导致仪器无法在高阻地层条件下工作的情况。
在一种可能的实施方式中,所述总体目标函数还包括控制电阻率测井仪器特征的第二目标函数,所述仪器特征包括仪器总长度、仪器主电极长度和仪器电极之间的最小绝缘长度中的任意至少一项,所述第二目标函数包括控制仪器总长度不超过Lmax的目标函数f5、控制仪器主电极A0长度LA0不小于LA0min的目标函数f6、控制仪器电极之间的绝缘长度L0不小于L0min的目标函数f7中的任意至少一项;
所述第二目标函数的公式为式中ti为第i个仪器特征的真实值,obji为第i个仪器特征的目标值。所述第二目标函数的公式分别为:根据算法计算的仪器尺寸结构计算f5~f7。
基于所述每个目标函数包括至少一个变量,所述各个目标函数的权重为其所有变量对应值之和;针对本实施例中的目标函数f1,其变量为各探测模式的伪几何因子Gk(r),其中rk为第k种探测模式的探测深度,为预先设定的优化期望值,M为探测模式数量,Gk(rk)为当前优化参数的测井仪器第k种探测模式在r=rk时的伪几何因子;针对本实施例中的目标函数f2,其变量为各探测模式的井眼影响,为仪器的M种探测模式,ti′,j′,k′,1为井眼影响控制目标,Rt1,Ra1分别为图4所示井眼影响地层模型下计算的地层真电阻率、视电阻率;针对本实施例中的目标函数f3,其变量为各探测模式的围岩影响,其对应的为仪器的M种探测模式,ti′,j′,k′,2为围岩影响控制目标,Rt2,Ra2分别为图5所示围岩影响地层模型下计算的地层真电阻率、视电阻率;针对本实施案例中的目标函数f4,其变量为各屏蔽电极电流,其对应的为其中N为仪器屏蔽电极数量;针对本实施例中的目标函数f5,其变量为仪器总长度,其对应的为2L;针对本实施例中的目标函数f6,其变量为仪器主电极A0长度,其对应的为2LA0;针对本实施例中的目标函数f7,其变量为仪器电极之间的最小绝缘长度,其对应的为2L0。将以上比例归一化后,最终多个目标函数变量权重确定为
将以上各比例因子的最小比例归一化,如果最大比例大于预先设定的目标函数权重调整最小比例因子Wmin(Wmin>1),将多个目标函数权重调整为最小比例归一化后的对应值,调整成功,否则调整不成功。
经多次迭代以及对所述多个目标函数的权重进行自适应调整后,优化结束时的状态包括:
所述总体目标函数收敛,优化结束;
所述总体目标函数不收敛,但所述多个目标函数的扰动引起所述总体目标函数的变化贡献相等,优化结束。
下面参考图2详细描述本第一实施例的较佳的实施流程图,包括:
S101,获取输入的多个因素;
S102,多个目标函数以预先设定的权重求和获得可通过算法构造并计算多个由电极系结构参数组构成的仪器模型的总体目标函数;
S103,根据算法构造并计算所述多个仪器模型的总体目标函数;
S104,判断总体目标函数是否收敛,是则执行S109,否则执行S105;
S105,计算总体目标函数的收敛速度并判断其是否小于预先设定的总体目标函数收敛速度的最小值,是则执行S106,否则执行S103;
S106,计算多个目标函数的扰动引起总体目标函数的变化贡献;
S107,比较多个目标函数的扰动引起总体目标函数的变化贡献是否相等,是则执行S109,否则执行S108;
S108,基于多个目标函数的扰动引起总体目标函数的变化贡献相等的原则自适应调整所述多个目标函数的权重,执行S103。
S109,优化结束。
下面参考图6详细描述本发明提供的一种电阻率测井仪器电极系结构参数智能优化装置的第二实施例:本实施例主要应用于电阻率测井仪器。
如图所示,本实施例提供的电阻率测井仪器电极系结构参数智能优化装置,包括:
因素获取模块201,用于获取输入的影响电阻率测井仪器电极系结构参数的多个因素;所述因素包括仪器探测深度、仪器井眼影响、仪器围岩影响和仪器主电极-屏蔽电极电流比中的任意至少一个;所述因素还可以包括仪器总长度、仪器主电极长度和仪器电极之间的最小绝缘长度中的任意至少一项;
预处理模块202,用于所述多个因素对应的多个目标函数以预先设定的权重求和获得可通过算法构造并计算多个由电极系结构参数组构成的仪器模型的总体目标函数;所述多个目标函数包括第一目标函数,还可以包括第二目标函数,所述多个目标函数以预先设定的权重求和获得总体目标函数,可以预先设定多个目标函数的初始权重为1;
权重自适应调整模块203,用于根据算法构造并计算所述多个仪器模型的总体目标函数,并判断所述总体目标函数是否满足所述多个目标函数的扰动对其的影响相等的第二优化终止条件,如果满足,优化结束,输出优化结果,否则,自适应调整所述多个目标函数在所述总体目标函数中的权重以使所述多个目标函数的扰动对所述总体目标函数的影响相等,避免某些目标函数值过大掩盖了某些目标函数值很小但对电阻率测井仪器的电极系结构参数影响同样重要的因素导致优化的效果片面,从而获得对电阻率测井仪器电极系结构参数影响更全面的优化效果;根据算法构造并计算多个仪器模型的总体目标函数时采用迭代计算;每重复计算一次即为一次迭代,每次迭代得到的结果作为下一次迭代的初始值;根据算法生成多组用于构造仪器模型的电阻率测井仪器电极系结构参数后,可采用并行正演计算所述多组仪器模型对应的总体目标函数,对不收敛的总体目标函数计算其收敛速度,所述并行正演计算可以一次计算多个仪器模型的总体目标函数,达到提高计算速度并扩大计算规模,使得计算效率大大提高;
所述算法包括遗传算法、人工神经网络、蚁群算法、粒子群算法、模拟退火算法、差分进化算法中的一种或者多种;
优化结果输出模块204,用于输出满足所述优化终止条件的仪器电极系结构参数。
在一种可能的实施方式中,所述权重自适应调整模块包括含有所述第二优化终止条件的第二优化终止单元,还包括含有所述第一优化终止条件的第一优化终止单元:
所述第一优化终止单元,包括第一优化终止条件,用于当根据算法计算出的所述总体目标函数收敛时,优化结束,当其不收敛时,执行计算其收敛速度;
当计算出的所述总体目标函数的收敛速度与预先设定的总体目标函数收敛速度的最小值相比较数值较大时,则执行返回继续计算所述总体目标函数;
所述总体目标函数的收敛速度的公式为式中Fi表示第i次优化的总体目标函数;vmin为预先定义的总体目标函数收敛速度最小值,其值大于1;
所述第二优化终止单元,包括第二优化终止条件,用于当计算出的所述总体目标函数的收敛速度与预先设定的总体目标函数收敛速度的最小值相比较数值较小时,计算所述多个目标函数的扰动引起所述总体目标函数的变化贡献是否相等,若相等,优化结束;
当计算出的所述多个目标函数的扰动引起所述总体目标函数的变化贡献不相等时,则基于所述多个目标函数的扰动引起所述总体目标函数的变化贡献相等的原则执行对所述多个目标函数的权重进行自适应调整,权重调整后,将执行按新的权重返回继续计算所述总体目标函数;
在一种可能的实施方式中,所述第二优化终止单元中的计算多个目标函数的扰动引起总体目标函数的变化贡献,包括:
第一计算单元,用于基于数据误差传递理论,所述多个目标函数的扰动体现为所述多个目标函数变量的扰动,以及所述多个目标函数变量的扰动相等,建立所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式;
所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式为:其中,Fi为根据所述多个目标函数中的变量u1,u2,...,un,由所述多个目标函数按F=F(u1,u2,...,un)合成的所述总体目标函数,n为变量的数量,i为整个优化过程的迭代次数;
根据所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式计算得出所述多个目标函数的扰动引起所述总体目标函数的变化贡献为 所述基于数据误差传递理论,即为基于测量结果的合成标准不确定度由各基本测量数据的误差传递而来。
在一种可能的实施方式中,所述第二优化终止单元中的基于所述多个目标函数的扰动引起所述总体目标函数的变化贡献相等的原则自适应调整所述多个目标函数的权重,包括:
第二计算单元,用于根据所述多个目标函数中每个目标函数的变量比例调整所述多个目标函数的权重比例,以使调整后所述多个目标函数中每个目标函数的变量比例满足将所述比例中的最小比例归一化,若最大比例大于Wmin,将所述多个目标函数的权重调整为最小比例归一化后的对应值,调整成功,否则调整不成功,其中Wmin为预先设定的所述多个目标函数权重调整的最小比例因子,其值大于1。
表示所述总体目标函数对第j个目标函数变量的偏导在第i次优化计算的目标函数fj(i)的取值,这个值决定了当前状态下第j个目标函数对总体目标函数的贡献,即第j个目标函数变量的扰动对总体目标函数的影响,将目标函数权重比例调整为即实现了所述多个目标函数变量对所述总体目标函数的贡献相等,使得优化过程中各个变化因素对所述总体目标函数的影响无差别,避免了某些目标函数值过大掩盖了某些目标函数值很小但对电阻率测井仪器的特性影响同样重要的因素的情况,达到了更全面的优化效果。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (10)
1.一种电阻率测井仪器电极系结构参数智能优化的方法,其特征在于,包括:
在因素获取模块中,获取输入的影响电阻率测井仪器电极系结构参数的多个因素;
在预处理模块中,根据所述多个因素对应的多个目标函数以预先设定的权重求和获得可通过算法构造并计算多个由电极系结构参数组构成的仪器模型的总体目标函数;
在权重自适应调整模块中,根据算法构造并计算所述多个仪器模型的总体目标函数,并判断所述总体目标函数是否满足所述多个目标函数的扰动对其的影响相等的第二优化终止条件,如果满足,优化结束,输出优化结果,否则,基于所述多个目标函数的扰动对所述总体目标函数的影响相等的原则执行对所述多个目标函数的权重进行自适应调整,权重调整后,将执行按新的权重返回继续计算所述总体目标函数;
在优化结果输出模块中,输出满足所述优化终止条件的仪器电极系结构参数。
2.根据权利要求1所述的电阻率测井仪器电极系结构参数智能优化的方法,其特征在于,所述在权重自适应调整模块中的第二优化终止条件之前,还包括第一优化终止条件:
所述第一优化终止条件为当根据算法计算出的所述总体目标函数收敛时,优化结束,当其不收敛时,执行计算其收敛速度;
当计算出的所述总体目标函数的收敛速度与预先设定的总体目标函数收敛速度的最小值相比较数值较大时,则执行返回继续计算所述总体目标函数;
所述第二优化终止条件为当计算出的所述总体目标函数的收敛速度与预先设定的总体目标函数收敛速度的最小值相比较数值较小时,计算所述多个目标函数的扰动引起所述总体目标函数的变化贡献是否相等,若相等,优化结束;
当计算出的所述多个目标函数的扰动引起所述总体目标函数的变化贡献不相等时,则基于所述多个目标函数的扰动引起所述总体目标函数的变化贡献相等的原则执行对所述多个目标函数的权重进行自适应调整,权重调整后,将执行按新的权重返回继续计算所述总体目标函数。
3.根据权利要求2所述的电阻率测井仪器电极系结构参数智能优化的方法,其特征在于,所述计算多个目标函数的扰动引起所述总体目标函数的变化贡献,包括:
基于数据误差传递理论,所述多个目标函数的扰动体现为所述多个目标函数变量的扰动,以及所述多个目标函数变量的扰动相等,建立所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式;
所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式为:其中,Fi为根据所述多个目标函数中的变量u1,u2,...,un,由所述多个目标函数按F=F(u1,u2,...,un)合成的所述总体目标函数,n为变量的数量,i为整个优化过程的迭代次数;
根据所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式计算得出所述多个目标函数的扰动引起所述总体目标函数的变化贡献为
4.根据权利要求3所述的电阻率测井仪器电极系结构参数智能优化的方法,其特征在于,所述基于所述多个目标函数的扰动引起所述总体目标函数的变化贡献相等的原则执行对所述多个目标函数的权重进行自适应调整,包括:
根据所述多个目标函数中每个目标函数的变量比例调整所述多个目标函数的权重比例,以使调整后所述多个目标函数中每个目标函数的变量比例满足将所述比例中的最小比例归一化,若最大比例大于Wmin,将所述多个目标函数的权重调整为最小比例归一化后的对应值,调整成功,否则调整不成功,其中Wmin为预先设定的所述多个目标函数权重调整的最小比例因子,其值大于1。
5.根据权利要求1-4中任一项权利要求所述的电阻率测井仪器电极系结构参数智能优化的方法,其特征在于,所述多个目标函数包括第一目标函数,所述第一目标函数包括:控制仪器探测深度的伪几何因子目标函数f1、控制仪器井眼影响的影响因子目标函数f2、控制仪器围岩影响的影响因子目标函数f3和控制仪器主电极-屏蔽电极电流比的电极电流目标函数f4中的任意至少一项;所述第一目标函数的计算方法为:
所述伪几何因子目标函数f1为不同探测模式的探测深度目标函数以预先设定的权重求和,其公式为:式中W1,i为第i个探测模式目标函数的权重,f1,i为第i个探测模式探测深度目标函数,目标函数f1,i的公式为式中ri为所优化的电阻率测井仪器第i种探测模式期望的探测深度,其中仪器探测深度定义为伪几何因子为0.5时的侵入半径r,式中Rt为伪几何因子地层模型原状地层电阻率,Ra(r)为伪几何因子地层模型在侵入半径r下的电极系对地层模型响应的视电阻率,Rxo为伪几何因子地层模型侵入带电阻率;
所述影响因子目标函数f2为不同探测模式的井眼影响目标函数以预先设定的权重求和,其公式为:式中W2,i,j,k为第i个探测模式、第j个控制点和第k个井眼半径目标函数的权重,f2,i,j,k为第i个探测模式、第j个控制点xi,j第k个井眼半径下的井眼影响目标函数,目标函数f2,i,j,k的公式为Rt=xi,j,wr=wrk,式中ti,j,k为不同控制点xi,j下不同井眼半径的井眼校正系数的控制目标,通常接近1,式中Rt/Ra为井眼影响地层模型井眼校正系数,Rt为井眼影响地层模型原状地层电阻率,Ra为井眼影响地层模型利用电阻率为Rt的均匀地层计算的视电阻率,wr为井眼影响地层模型井眼半径;其中仪器井眼影响通过不同Ra/Rm情况下的井眼校正系数Rt/Ra衡量,Rm为井眼影响地层模型井眼电阻率;理想情况下井眼校正系数为1;
所述影响因子目标函数f3为不同探测模式的围岩影响目标函数以预先设定的权重求和,其公式为:式中W3,i,j,k为第i个探测模式、第j个控制点和第k个围岩/目的层对比度下目标函数的权重,f3,i,j,k为第i个探测模式、第j个控制点hi,j和第k个围岩/目的层对比度下围岩影响目标函数,目标函数f3,i,j,k的公式为h=hi,j,式中ti,j,k为不同控制点hi,j下不同围岩/目的层对比度下围岩校正系数的控制目标,通常接近1,式中Rt为围岩影响地层模型目的层电阻率,Rs为围岩影响地层模型围岩电阻率,Rt/Ra为围岩影响地层模型围岩校正系数,Rt/Rs为围岩影响地层模型地层对比度;其中围岩影响通过不同目的层厚度下的围岩校正系数Rt/Ra衡量,理想情况下围岩校正系数为1;
所述电极电流目标函数f4的公式为设主电极电流为单位电流,式中Ii为以仪器测量最高电阻率为均匀地层、最深探测模式第i个屏蔽电极的电流,threshold为设定的主电极-屏蔽电极电流比的最小值;其中仪器主电极-屏蔽电极电流比反映了仪器聚焦情况下所需屏蔽电流大小,其比值越小,电路硬件聚焦越容易实现;设地层为均匀地层,地层真电阻率为仪器工作的最高电阻率地层,在主电极发射1A电流情况下,计算此时仪器最远探测模式下各电极电流,进而根据电极电流目标函数f4的定义计算该目标函数;根据仪器探测的基本性质,仪器工作的目的层电阻率越高,探测深度越深,达到聚焦条件所需的屏蔽电流越大,因此在计算电极电流目标函数f4时,仅考虑仪器测量的最高电阻率和最深探测模式,使得优化后的仪器也可以在高阻地层条件下工作。
6.根据权利要求5所述的电阻率测井仪器电极系结构参数智能优化的方法,其特征在于,所述总体目标函数还包括控制电阻率测井仪器特征的第二目标函数,所述仪器特征包括仪器总长度、仪器主电极长度和仪器电极之间的最小绝缘长度中的任意至少一项。
7.一种电阻率测井仪器电极系结构参数智能优化装置,其特征在于,包括;
因素获取模块,用于获取输入的影响电阻率测井仪器电极系结构参数的多个因素;
预处理模块,用于所述多个因素对应的多个目标函数以预先设定的权重求和获得可通过算法构造并计算多个由电极系结构参数组构成的仪器模型的总体目标函数;
权重自适应调整模块,用于根据算法构造并计算所述多个仪器模型的总体目标函数,并判断所述总体目标函数是否满足所述多个目标函数的扰动对其的影响相等的第二优化终止条件,如果满足,优化结束,输出优化结果,否则,基于所述多个目标函数的扰动对所述总体目标函数的影响相等的原则执行对所述多个目标函数的权重进行自适应调整,权重调整后,将执行按新的权重返回继续计算所述总体目标函数;
优化结果输出模块,用于输出满足所述优化终止条件的仪器电极系结构参数。
8.根据权利要求7所述的电阻率测井仪器电极系结构参数智能优化装置,其特征在于,所述权重自适应调整模块包括含有所述第二优化终止条件的第二优化终止单元,还包括含有所述第一优化终止条件的第一优化终止单元:
所述第一优化终止单元,包括第一优化终止条件,用于当根据算法计算出的所述总体目标函数收敛时,优化结束,当其不收敛时,执行计算其收敛速度;
当计算出的所述总体目标函数的收敛速度与预先设定的总体目标函数收敛速度的最小值相比较数值较大时,则执行返回继续计算所述总体目标函数;
所述第二优化终止单元,包括第二优化终止条件,用于当计算出的所述总体目标函数的收敛速度与预先设定的总体目标函数收敛速度的最小值相比较数值较小时,计算所述多个目标函数的扰动引起所述总体目标函数的变化贡献是否相等,若相等,优化结束;
当计算出的所述多个目标函数的扰动引起所述总体目标函数的变化贡献不相等时,则基于所述多个目标函数的扰动引起所述总体目标函数的变化贡献相等的原则执行对所述多个目标函数的权重进行自适应调整,权重调整后,将执行按新的权重返回继续计算所述总体目标函数。
9.根据权利要求8所述的电阻率测井仪器电极系结构参数智能优化装置,其特征在于,所述第二优化终止单元中的计算多个目标函数的扰动引起总体目标函数的变化贡献,包括:
第一计算单元,用于基于数据误差传递理论,所述多个目标函数的扰动体现为所述多个目标函数变量的扰动,以及所述多个目标函数变量的扰动相等,建立所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式;
所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式为:其中,Fi为根据所述多个目标函数中的变量u1,u2,...,un,由所述多个目标函数按F=F(u1,u2,...,un)合成的所述总体目标函数,n为变量的数量,i为整个优化过程的迭代次数;
根据所述多个目标函数变量的扰动引起所述总体目标函数的变化贡献的公式计算得出所述多个目标函数的扰动引起所述总体目标函数的变化贡献为
10.根据权利要求9所述的电阻率测井仪器电极系结构参数智能优化装置,其特征在于,所述第二优化终止单元中的基于所述多个目标函数的扰动引起所述总体目标函数的变化贡献相等的原则自适应调整所述多个目标函数的权重,包括:
第二计算单元,用于根据所述多个目标函数中每个目标函数的变量比例调整所述多个目标函数的权重比例,以使调整后所述多个目标函数中每个目标函数的变量比例满足将所述比例中的最小比例归一化,若最大比例大于Wmin,将所述多个目标函数的权重调整为最小比例归一化后的对应值,调整成功,否则调整不成功,其中Wmin为预先设定的所述多个目标函数权重调整的最小比例因子,其值大于1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910132183.5A CN109711109B (zh) | 2019-02-22 | 2019-02-22 | 一种电阻率测井仪器电极系结构参数智能优化的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910132183.5A CN109711109B (zh) | 2019-02-22 | 2019-02-22 | 一种电阻率测井仪器电极系结构参数智能优化的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109711109A true CN109711109A (zh) | 2019-05-03 |
CN109711109B CN109711109B (zh) | 2023-12-22 |
Family
ID=66264941
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910132183.5A Active CN109711109B (zh) | 2019-02-22 | 2019-02-22 | 一种电阻率测井仪器电极系结构参数智能优化的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109711109B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114526063A (zh) * | 2022-02-14 | 2022-05-24 | 北京工业大学 | 获取探边电磁波测井仪器结构参数的方法及装置 |
CN117930787A (zh) * | 2024-03-21 | 2024-04-26 | 南京航空航天大学 | 一种数控机床加工的工艺参数优化方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102678106A (zh) * | 2012-05-02 | 2012-09-19 | 中国电子科技集团公司第二十二研究所 | 随钻电磁波电阻率测井仪器的数据处理方法 |
CN104481525A (zh) * | 2014-11-28 | 2015-04-01 | 中国石油天然气集团公司 | 一种快测阵列感应测井系统及其测井方法 |
CN106837299A (zh) * | 2017-01-19 | 2017-06-13 | 中国科学院声学研究所 | 一种井眼校正的系统及方法 |
CN107630697A (zh) * | 2017-09-26 | 2018-01-26 | 长江大学 | 基于随钻电磁波电阻率测井的地层电阻率联合反演方法 |
CN107947196A (zh) * | 2017-11-16 | 2018-04-20 | 国网四川省电力公司 | 一种基于改进粒子群算法的超低频振荡抑制方法 |
CN108073765A (zh) * | 2017-12-11 | 2018-05-25 | 中海石油(中国)有限公司湛江分公司 | 一种水平井常规随钻测井地层界面识别与边界距反演方法 |
-
2019
- 2019-02-22 CN CN201910132183.5A patent/CN109711109B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102678106A (zh) * | 2012-05-02 | 2012-09-19 | 中国电子科技集团公司第二十二研究所 | 随钻电磁波电阻率测井仪器的数据处理方法 |
CN104481525A (zh) * | 2014-11-28 | 2015-04-01 | 中国石油天然气集团公司 | 一种快测阵列感应测井系统及其测井方法 |
CN106837299A (zh) * | 2017-01-19 | 2017-06-13 | 中国科学院声学研究所 | 一种井眼校正的系统及方法 |
CN107630697A (zh) * | 2017-09-26 | 2018-01-26 | 长江大学 | 基于随钻电磁波电阻率测井的地层电阻率联合反演方法 |
CN107947196A (zh) * | 2017-11-16 | 2018-04-20 | 国网四川省电力公司 | 一种基于改进粒子群算法的超低频振荡抑制方法 |
CN108073765A (zh) * | 2017-12-11 | 2018-05-25 | 中海石油(中国)有限公司湛江分公司 | 一种水平井常规随钻测井地层界面识别与边界距反演方法 |
Non-Patent Citations (6)
Title |
---|
VAHIDREZA NASIRIAN: "Distributed Cooperative Control of DC Microgrids", IEEE TRANSACTIONS ON POWER ELECTRONICS * |
刘文颖: "一种计及交叉权重的有功校正控制算法", 电工技术学报 * |
康正明: ""钻头电阻率测井仪器探测特性研究"", pages 459 - 460 * |
苏航: ""基于时变条件的钢筋混凝土桥梁耐久性AHP-模糊综合预测研究"", pages 9 - 32 * |
郭涛: ""阵列侧向测井井眼影响研究"", pages 24 - 53 * |
顾明亮: "基于动态调整惯性权重的混合粒子群算法", 计算机与现代化 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114526063A (zh) * | 2022-02-14 | 2022-05-24 | 北京工业大学 | 获取探边电磁波测井仪器结构参数的方法及装置 |
CN114526063B (zh) * | 2022-02-14 | 2024-07-12 | 北京工业大学 | 获取探边电磁波测井仪器结构参数的方法及装置 |
CN117930787A (zh) * | 2024-03-21 | 2024-04-26 | 南京航空航天大学 | 一种数控机床加工的工艺参数优化方法 |
CN117930787B (zh) * | 2024-03-21 | 2024-06-11 | 南京航空航天大学 | 一种数控机床加工的工艺参数优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109711109B (zh) | 2023-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107239859B (zh) | 基于串联长短期记忆循环神经网络的供热负荷预测方法 | |
CN109376850A (zh) | 一种基于改进bp神经网络状态估计中不良数据的检测方法 | |
CN108920812B (zh) | 一种机械加工表面粗糙度预测方法 | |
CN111191835A (zh) | 基于c-gan迁移学习的ies不完备数据负荷预测方法及系统 | |
CN110650153B (zh) | 一种基于聚焦损失深度神经网络的工控网络入侵检测方法 | |
CN110262233A (zh) | 一种基于遗传算法和bp神经网络的磁控镀膜仪工艺参数的优化方法 | |
CN109711109A (zh) | 一种电阻率测井仪器电极系结构参数智能优化的方法及装置 | |
CN110348046B (zh) | 一种基于模糊综合评价的钻进轨迹多目标优化方法 | |
CN110472840A (zh) | 一种基于神经网络技术的农业水利调度方法及系统 | |
CN113300380B (zh) | 一种基于负荷曲线分段的配电网无功优化补偿方法 | |
CN105956715A (zh) | 一种土壤墒情预测方法及装置 | |
CN108564592A (zh) | 基于动态多种群集成差分进化算法的图像分割方法 | |
CN107868979A (zh) | 一种基于恒拉速控制结构的硅单晶直径控制方法 | |
CN108204944A (zh) | 基于apso优化的lssvm的埋地管道腐蚀速率预测方法 | |
CN103092074A (zh) | 半导体先进过程控制的参数优化控制方法 | |
CN115048804A (zh) | 一种城市管网多传感器优化部署方法 | |
CN109919313A (zh) | 一种梯度传输的方法及分布式训练系统 | |
CN112149883A (zh) | 基于fwa-bp神经网络的光伏功率预测方法 | |
CN115132302B (zh) | 一种pvc改性聚合物的研发配方验证方法及系统 | |
CN109472105A (zh) | 半导体产品良率上限分析方法 | |
CN110147890A (zh) | 一种基于狮群算法优化极限学习机集成学习的方法及系统 | |
CN111914487B (zh) | 一种基于对抗神经网络的无资料地区水文参数率定方法 | |
CN110174841A (zh) | 基于参数优化的工业过程大滞后惯性系统平行控制方法 | |
CN107392315B (zh) | 一种优化大脑情感学习模型的乳腺癌数据分类方法 | |
CN111242463B (zh) | 一种基于bp神经网络的at单线供电系统的故障测距方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |