CN104407390A - 一种确定海洋可控电磁源勘探激发频率的方法及系统 - Google Patents
一种确定海洋可控电磁源勘探激发频率的方法及系统 Download PDFInfo
- Publication number
- CN104407390A CN104407390A CN201410697345.7A CN201410697345A CN104407390A CN 104407390 A CN104407390 A CN 104407390A CN 201410697345 A CN201410697345 A CN 201410697345A CN 104407390 A CN104407390 A CN 104407390A
- Authority
- CN
- China
- Prior art keywords
- stimulating frequency
- oily
- earth
- curve
- module
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请提供一种确定海洋可控电磁源勘探激发频率的方法及系统。所述方法包括:建立一维地电模型;利用快速正演算法,计算含油气以及不含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线;根据含油气以及不含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线,计算不同激发频率点对应的归一化异常曲线;对所述不同激发频率点对应的归一化异常曲线进行编制,得到归一化异常幅值随激发频率和偏移距变化的平面分布图;根据所述平面分布图确定最佳的发射源激发频率。本申请实施例通过快速正演算法,得到归一化异常幅值随偏移距和激发频率点变化的平面分布图,从而可以精确地确定最佳的发射源激发频率点。
Description
技术领域
本申请涉及海洋可控电磁源勘探技术领域,特别涉及一种确定海洋可控电磁源勘探激发频率的方法及系统。
背景技术
近年来,随着陆上资源的逐步枯竭和科学技术的高速发展,海洋资源因其经济上的巨大潜力和战略上的重要地位越来越被人们所重视。海洋资源中的海洋油气资源更是各国的战略储备资源之首。在对海洋油气资源的勘探过程中,海洋可控电磁源勘探技术(CSEM,Controlled-Source Electromagnetic Method)作为一种新型的海洋油气资源勘探技术,日益受到各国研究人员的关注。海洋可控电磁源勘探技术可以通过对海底地表电磁场的观测,有效地追踪到与含油气层有关的异常信息,从而提高海洋勘探的成功率。
在实施海洋可控电磁源勘探工作时,确定合适的发射源激发频率对获取海底含油气层的异常信息具有重要作用。现有技术对发射源激发频率进行选择的一般做法为:选取不同的发射源激发频率点,针对每个发射源激发频率点计算出随偏移距变化的归一化异常曲线。通过比较不同发射源激发频率点对应的异常幅值的大小,从而确定出最佳的发射源激发频率点。
在实现本申请的过程中,发明人发现现有技术中至少存在以下问题:
现有技术中通过发射源激发频率点对应的异常幅值进行对比,由于计算的激发频率点相对离散并且数量有限,不利于对最佳发射源激发频率点的确定。
发明内容
本申请实施例的目的在于提供一种确定海洋可控电磁源勘探激发频率的方法及系统,能够精确地确定最佳的发射源激发频率点。
本申请实施例提供的一种确定海洋可控电磁源勘探激发频率的方法及系统是这样实现的:
一种确定海洋可控电磁源勘探激发频率的方法,包括:
建立一维地电模型,所述一维地电模型包括:含油气的地电模型以及不含油气的地电模型;
利用快速正演算法,计算含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线;
利用快速正演算法,计算不含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线;
根据含油气以及不含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线,计算不同激发频率点对应的归一化异常曲线;
对所述不同激发频率点对应的归一化异常曲线进行编制,得到归一化异常幅值随激发频率和偏移距变化的平面分布图;
根据所述平面分布图确定最佳的发射源激发频率。
一种确定海洋可控电磁源勘探激发频率的系统,包括:一维地电模型建立模块,快速正演计算模块,归一化异常曲线计算模块,平面分布图编制模块,最佳激发频率确定模块,其中:
所述一维地电模型建立模块,用来建立一维地电模型,所述一维地电模型包括:含油气的地电模型以及不含油气的地电模型;
所述快速正演计算模块,用来利用快速正演算法,计算含油气以及不含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线;
所述归一化异常曲线计算模块,用来根据含油气以及不含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线,计算不同激发频率点对应的归一化异常曲线;
所述平面分布图编制模块,用来对所述不同激发频率点对应的归一化异常曲线进行编制,得到归一化异常幅值随激发频率和偏移距变化的平面分布图;
所述最佳激发频率确定模块,用来根据所述平面分布图确定最佳的发射源激发频率。
本申请实施例通过快速正演算法,计算出大量不同激发频率点对应的随偏移距变化的归一化异常曲线,并且将得到的归一化的异常曲线编制为归一化异常幅值随偏移距和激发频率点变化的平面分布图,从而可以精确地确定最佳的发射源激发频率点。
附图说明
图1是本申请一实施例提供的一种确定海洋可控电磁源勘探激发频率的方法流程图;
图2为本申请一实施例中建立的一维地电模型示意图;
图3为激发频率为5Hz时含油气以及不含油气的地电模型中的MVO曲线示意图;
图4为激发频率为5Hz时归一化后的异常曲线;
图5为本申请一实施例中编制的归一化异常幅值随激发频率和偏移距变化的平面分布图;
图6a为储层含油气饱和度30%以下的电磁异常平面分布图;
图6b为储层含油气饱和度40%-60%的电磁异常平面分布图;
图6c为储层含油气饱和度60%以上的电磁异常平面分布图;
图7为本申请一实施例提供的一种确定海洋可控电磁源勘探激发频率的系统的功能模块图;
图8为本申请一实施例提供的一种确定海洋可控电磁源勘探激发频率的系统中快速正演计算模块的功能模块图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
图1是本申请一实施例提供的一种确定海洋可控电磁源勘探激发频率的方法流程图。如图1所示,所述方法包括:
S1:建立一维地电模型。
图2为本申请一实施例中建立的一维地电模型示意图。如图2所示,一维地电模型包括空气,海水以及海底地层。靠近海底的发射源发射的低频电磁信号,不仅向水中传播,同时还向水面以及向海底地层传播。因此,位于海底的电磁接收站可以接收到三条路径的信号,一是直达波,而是来自海底地层的反射和折射信号,三是来自水面的反射和折射信号。
S2:利用快速正演算法,计算含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线。
在海洋可控电磁源勘探技术中,通过可以用电磁场振幅随偏移距变化的MVO(MagnitudeValue Offset)曲线来显示含油气地层异常。在计算发射源产生的电磁场响应时,通常可以采用有限元正演方法来实现。传统的有限元正演方法中,一般是采用矩形剖分有限元的方法,将海底地层结构利用矩形进行划分,进而对发射源产生的电磁场响应进行求解。在本申请的一实施例中,可以采用不规则三角剖分有限元的快速正演算法,对发射源产生的电磁场响应进行求解。具体方法如下所述:
S201:采用非规则的三角单元对地电模型结构进行网格剖分。
采用非规则三角单元进行网格剖分的优势在于:
1、对比传统的矩形单元或规则三角剖分,非规则三角形剖分能更容易模拟复杂的、任意形状的地电模型结构。
2、非规则三角形剖分可以很好对不同尺度的模型结构进行剖分,例如:在模型变化剧烈的地方采用加密网格剖分,而在模型变化平缓的地方采用稀疏网格剖分,从而可以减少网格节点个数,进而减少对内存大小的要求。在保证精度的条件下,提高计算速度。
3、采用非规则三角形剖分可以不必对模型进行扩边处理,减小对内存的要求。
S202:利用自适应的有限元计算方法对网格剖分后的海底地层结构进行计算。
对地电模型进行非规则的三角单元剖分后,可以从麦克斯韦方程入手,进行自适应的有限元计算方法。所述麦克斯韦方程具体表示为:
▽×E=iμ0ωH
▽×H-σE=Js
其中,E代表电场响应,H代表磁场响应,Js为电流密度,σ为电导率,μ0为真空磁导率。
利用加权残值法(Zienkiewicz,1977)以及齐次Dirichlet边界条件,可以将总的电磁场分为若干个非规则三角元的和。其结果具体可以表示为:
Ku=p
其中K为常数,u为列向量,包括各个非规则三角单元的电磁场响应,p为已知分量。
通过上式可以构建出一系列线性方程组,从而求解出每个非规则三角单元的第磁场响应。
S203:利用后验误差估计方法,对网格产生的数值进行误差估计。
利用有限元法进行计算的结果与剖分网格的大小有关,误差会随着网格剖分的细化程度的加深而减小。为了确保得到精确的值,就需要在地电模型结构中的每一个地方都进行足够小的网格剖分。但这样会导致大数据量的产生,造成计算效率下降的问题。本申请实施例中可以采用后验误差估计的方法,将误差超过预设阈值的网格再一次进行更加精细的网格剖分,经过重复迭代的过程,直到得到符合误差允许条件为止。
本申请一实施例中计算网格误差具体可以表示为:
其中,e代表剖分后的网格,ηe为网格e的后验误差,R为恢复算子,uh为求解出的网格e的电磁场响应,为求解出的网格e的电磁场响应的梯度,L2表示欧几里德空间中矢量的范数内积算子。
进一步地,可以将恢复算子R用SmQh来表示,其中,Qh是L2的投影算子,S为平滑算子,m是平滑处理的迭代次数。这样,计算网格误差具体可以改进为:
S204:对误差超过预设阈值的网格重新进行网格剖分和计算,直到误差估计不大于预设阈值为止。
S3:利用快速正演算法,计算不含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线。
与步骤S2采用的方法类似,可以得到不含油气的地电模型中的不同激发频率点对应的MVO曲线。现在可以将同一激发频率点对应的含油气以及不含油气的地电模型的MVO曲线绘制到同一坐标系下。例如,图3为激发频率为5Hz时含油气以及不含油气的地电模型中的MVO曲线示意图。如图3所示,实线代表含油气的地电模型的MVO曲线,虚线代表不含油气的地电模型的MVO曲线示意图。从图3中可以看出,在偏移距较小时,海底的电磁接收站仅仅能接收到直达波的信号,因此两条MVO曲线都重合在一起。当偏移距与海底高阻储层埋深相当或稍大时,在含油气的地电模型中,向海底地层传播的电磁波在遇到高阻油气层时,可以沿着岩层传播并且信号能量衰减得较慢。最终沿着岩层传播的电磁波折射回海底电磁接收站后,能量比直达波以及向上经海水面传播的电磁波的能量均较大,在三种信号中占据主导地位,反映到MVO曲线中就会呈现比较缓慢的下降趋势。而在不含油气的地电模型中,由于三种信号在传播过程中均衰减较快,因此反映到MVO曲线中就会呈现比较快速的下降趋势。这样两条MVO曲线便会出现分离。当偏移距进一步增大时,向下经海底传播的电磁波经过长距离较快的衰减,比向上经海水面传播的电磁波要弱得多。此时向上经海水面传播的电磁波就占据主导地位,两条MVO曲面又重新重合在一起。
S4:根据含油气以及不含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线,计算不同激发频率点对应的归一化异常曲线。
为了确定最佳的发射源激发频率点,通常要对含油气与不含油气的MVO曲线进行归一化,得到不同激发频率点对应的归一化异常曲线。具体方法如下所述:
在同一激发频率点,将含油气的MVO曲线与不含油气的MVO曲线中相同偏移距处对应的电磁场响应相除,便可以得到随偏移距变化的归一化后的异常幅值。图4为与图3对应的归一化后的异常曲线。对比图3和图4可以看出,当偏移距较小时,由于两条MVO曲线重合,因此得到的归一化异常幅值为1;当偏移距与海底高阻储层埋深相当或稍大时,由于两条MVO曲线出现分离,因此进行相除运算得到的归一化异常幅值会突起;当偏移距进一步增大时,由于两条MVO曲线重新重合,因此得到的归一化异常幅值又成为1。
S5:对所述不同激发频率点对应的归一化异常曲线进行编制,得到归一化异常幅值随激发频率和偏移距变化的平面分布图,并且根据所述平面分布图确定最佳的发射源激发频率。
以发射源的激发频率作为横坐标,以偏移距作为纵坐标,以颜色的深浅代表归一化异常幅值的大小,编制归一化异常幅值随激发频率和偏移距变化的平面分布图。图5为本申请一实施例中编制的归一化异常幅值随激发频率和偏移距变化的平面分布图。从图5中可以非常直观地看出归一化异常幅值最大的点,从而可以方便地从平面分布图中确定归一化异常幅值最大的点对应的发射源激发频率以及对应的偏移距。所述归一化异常幅值最大的点对应的发射源激发频率就可以为最佳的发射源激发频率。
在本申请的另一实施例中,可以利用本申请提供的一种确定海洋可控电磁源勘探激发频率的方法来预测海底储层含油气饱和度的定量变化。利用上述实施例得到归一化异常幅值随激发频率和偏移距变化的平面分布图的方法,可以得到如图6a,图6b,图6c所示的含油气储层电磁异常随激发频率和偏移距变化的平面分布图。其中,图6a为储层含油气饱和度30%以下的电磁异常平面分布图,图6b为储层含油气饱和度40%-60%的电磁异常平面分布图,图6c为储层含油气饱和度60%以上的电磁异常平面分布图。从图中可以直观地看出,颜色越深的地方含油气越饱和,并且可以方便地找出含油气最饱和的点对应的发射源激发频率以及偏移距。
本申请实施例还提供一种确定海洋可控电磁源勘探激发频率的系统。图7为本申请一实施例提供的一种确定海洋可控电磁源勘探激发频率的系统的功能模块图。如图7所示,所述系统包括:
一维地电模型建立模块1,用来建立一维地电模型,所述一维地电模型包括:含油气的地电模型以及不含油气的地电模型;
快速正演计算模块2,用来利用快速正演算法,计算含油气以及不含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线;
归一化异常曲线计算模块3,用来根据含油气以及不含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线,计算不同激发频率点对应的归一化异常曲线;
平面分布图编制模块4,用来对所述不同激发频率点对应的归一化异常曲线进行编制,得到归一化异常幅值随激发频率和偏移距变化的平面分布图;
最佳激发频率确定模块5,用来根据所述平面分布图确定最佳的发射源激发频率。
图8为本申请一实施例提供的一种确定海洋可控电磁源勘探激发频率的系统中快速正演计算模块2的功能模块图.如图8所示,所述快速正演计算模块2包括:
网格剖分模块201,用来采用非规则的三角单元对地电模型结构进行网格剖分;
自适应有限元计算模块202,用来利用自适应的有限元计算方法对网格剖分后的海底地层结构进行计算;
后验误差估算模块203,用来利用后验误差估计方法,对网格产生的数值进行误差估计;
误差循环计算模块204,用来对误差超过预设阈值的网格重新进行网格剖分和计算,直到误差估计不大于预设阈值为止。
进一步地,在本申请的一优选实施例中,所述后验误差估计方法具体可以由下式表示:
其中,e代表剖分后的网格,ηe为网格e的后验误差,R为恢复算子,uh为求解出的网格e的电磁场响应,为求解出的网格e的电磁场响应的梯度,L2表示欧几里德空间中矢量的范数内积算子。
进一步地,可以将恢复算子R用SmQh来表示,其中,Qh是L2的投影算子,S为平滑算子,m是平滑处理的迭代次数。这样,计算网格误差具体可以改进为:
通过上述实施例可以发现,本申请实施例通过快速正演算法,计算出大量不同激发频率点对应的随偏移距变化的归一化异常曲线,并且将得到的归一化的异常曲线编制为归一化异常幅值随偏移距和激发频率点变化的平面分布图,从而可以精确地确定最佳的发射源激发频率点。
在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进(例如,对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable Gate Array,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片2。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现,它与程序开发撰写时所用的软件编译器相类似,而要编译之前的原始代码也得用特定的编程语言来撰写,此称之为硬件描述语言(HardwareDescription Language,HDL),而HDL也并非仅有一种,而是有许多种,如ABEL(AdvancedBoolean Expression Language)、AHDL(Altera Hardware Description Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware DescriptionLanguage)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby Hardware Description Language)等,目前最普遍使用的是VHDL(Very-High-Speed Integrated Circuit Hardware DescriptionLanguage)与Verilog2。本领域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。
控制器可以按任何适当的方式实现,例如,控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit,ASIC)、可编程逻辑控制器和嵌入微控制器的形式,控制器的例子包括但不限于以下微控制器:ARC 625D、AtmelAT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320,存储器控制器还可以被实现为存储器的控制逻辑的一部分。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本申请可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (8)
1.一种确定海洋可控电磁源勘探激发频率的方法,其特征在于,包括:
建立一维地电模型,所述一维地电模型包括:含油气的地电模型以及不含油气的地电模型;
利用快速正演算法,计算含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线;
利用快速正演算法,计算不含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线;
根据含油气以及不含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线,计算不同激发频率点对应的归一化异常曲线;
对所述不同激发频率点对应的归一化异常曲线进行编制,得到归一化异常幅值随激发频率和偏移距变化的平面分布图;
根据所述平面分布图确定最佳的发射源激发频率。
2.如权利要求1所述的一种确定海洋可控电磁源勘探激发频率的方法,其特征在于,所述快速正演算法具体为:
采用非规则的三角单元对地电模型结构进行网格剖分;
利用自适应的有限元计算方法对网格剖分后的海底地层结构进行计算;
利用后验误差估计方法,对网格产生的数值进行误差估计;
对误差超过预设阈值的网格重新进行网格剖分和计算,直到误差估计不大于预设阈值为止。
3.如权利要求2所述的一种确定海洋可控电磁源勘探激发频率的方法,其特征在于,所述后验误差估计方法具体为:
其中,e代表剖分后的网格,ηe为网格e的后验误差,R为恢复算子,uh为求解出的网格e的电磁场响应,为求解出的网格e的电磁场响应的梯度,L2表示欧几里德空间中矢量的范数内积算子。
4.如权利要求3所述的一种确定海洋可控电磁源勘探激发频率的方法,其特征在于,所述恢复算子R具体为:
R=SmQh
其中,Qh是L2的投影算子,S为平滑算子,m是平滑处理的迭代次数。
5.一种确定海洋可控电磁源勘探激发频率的系统,其特征在于,包括:一维地电模型建立模块,快速正演计算模块,归一化异常曲线计算模块,平面分布图编制模块,最佳激发频率确定模块,其中:
所述一维地电模型建立模块,用来建立一维地电模型,所述一维地电模型包括:含油气的地电模型以及不含油气的地电模型;
所述快速正演计算模块,用来利用快速正演算法,计算含油气以及不含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线;
所述归一化异常曲线计算模块,用来根据含油气以及不含油气的地电模型中不同激发频率点处电磁场振幅随偏移距变化的曲线,计算不同激发频率点对应的归一化异常曲线;
所述平面分布图编制模块,用来对所述不同激发频率点对应的归一化异常曲线进行编制,得到归一化异常幅值随激发频率和偏移距变化的平面分布图;
所述最佳激发频率确定模块,用来根据所述平面分布图确定最佳的发射源激发频率。
6.如权利要求5所述的一种确定海洋可控电磁源勘探激发频率的系统,其特征在于,所述快速正演计算模块包括:网格剖分模块,自适应有限元计算模块,后验误差估算模块,误差循环计算模块,其中:
所述网格剖分模块,用来采用非规则的三角单元对地电模型结构进行网格剖分;
所述自适应有限元计算模块,用来利用自适应的有限元计算方法对网格剖分后的海底地层结构进行计算;
所述后验误差估算模块,用来利用后验误差估计方法,对网格产生的数值进行误差估计;
所述误差循环计算模块,用来对误差超过预设阈值的网格重新进行网格剖分和计算,直到误差估计不大于预设阈值为止。
7.如权利要求6所述的一种确定海洋可控电磁源勘探激发频率的系统,其特征在于,所述后验误差估计方法具体为:
其中,e代表剖分后的网格,ηe为网格e的后验误差,R为恢复算子,uh为求解出的网格e的电磁场响应,为求解出的网格e的电磁场响应的梯度,L2表示。
8.如权利要求7所述的一种确定海洋可控电磁源勘探激发频率的系统,其特征在于,所述恢复算子R具体为:
R=SmQh
其中,Qh是L2的投影算子,S为平滑算子,m是平滑处理的迭代次数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410697345.7A CN104407390A (zh) | 2014-11-27 | 2014-11-27 | 一种确定海洋可控电磁源勘探激发频率的方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410697345.7A CN104407390A (zh) | 2014-11-27 | 2014-11-27 | 一种确定海洋可控电磁源勘探激发频率的方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104407390A true CN104407390A (zh) | 2015-03-11 |
Family
ID=52645034
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410697345.7A Pending CN104407390A (zh) | 2014-11-27 | 2014-11-27 | 一种确定海洋可控电磁源勘探激发频率的方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104407390A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107065017A (zh) * | 2016-12-02 | 2017-08-18 | 北京工业大学 | 一种激励源编码选择方法及系统 |
CN109977585A (zh) * | 2019-04-04 | 2019-07-05 | 中南大学 | 一种高精度大地电磁正演模拟方法 |
CN111948722A (zh) * | 2019-05-15 | 2020-11-17 | 中国石油天然气集团有限公司 | 海底含油气异常识别方法及装置 |
-
2014
- 2014-11-27 CN CN201410697345.7A patent/CN104407390A/zh active Pending
Non-Patent Citations (4)
Title |
---|
YUGUO LI ET AL.: "《Special Section — Marine Controlled-Source Electromagnetic Methods 2D marine controlled-source electromagnetic modeling:Part 1—An adaptive finite-element algorithm》", 《GEOPHYSICS》 * |
何展翔等: "《基于三维模拟的海洋CSEM资料处理》", 《地球物理学报》 * |
史增园: "《海洋可控源电磁法正演研究与实现》", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
张秀丽等: "《海洋CSEM法探测高阻储层的勘探参数设计》", 《海洋地质与第四纪地质》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107065017A (zh) * | 2016-12-02 | 2017-08-18 | 北京工业大学 | 一种激励源编码选择方法及系统 |
CN107065017B (zh) * | 2016-12-02 | 2019-02-01 | 北京工业大学 | 一种激励源编码选择方法及系统 |
CN109977585A (zh) * | 2019-04-04 | 2019-07-05 | 中南大学 | 一种高精度大地电磁正演模拟方法 |
CN111948722A (zh) * | 2019-05-15 | 2020-11-17 | 中国石油天然气集团有限公司 | 海底含油气异常识别方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3018502A2 (en) | Modeling fluid-conducting fractures in reservoir simulation grids | |
Trani et al. | Seismic history matching of fluid fronts using the ensemble Kalman filter | |
Costabile et al. | Hydraulic characterization of river networks based on flow patterns simulated by 2‐D shallow water modeling: Scaling properties, multifractal interpretation, and perspectives for channel heads detection | |
EP3096252A2 (en) | Adaptive multiscale multi-fidelity reservoir simulation | |
Lindsey et al. | Improved workflow for 3D inverse modeling of magnetotelluric data: Examples from five geothermal systems | |
Sharifi et al. | Novel permeability upscaling method using fast marching method | |
CN103278848A (zh) | 基于mpi并行预条件迭代的地震成像正演方法 | |
CN104407390A (zh) | 一种确定海洋可控电磁源勘探激发频率的方法及系统 | |
Kânoğlu et al. | Focusing of long waves with finite crest over constant depth | |
Paulatto et al. | Heterogeneous and asymmetric crustal accretion: New constraints from multibeam bathymetry and potential field data from the Rainbow area of the M id‐A tlantic R idge (36° 15'N) | |
CN108614293A (zh) | 砂体预测方法及装置 | |
Pareja‐Roman et al. | Impact of historical channel deepening on tidal hydraulics in the Delaware Estuary | |
CN104391319A (zh) | 一种地震资料采集系统的确定方法及装置 | |
CN110334365A (zh) | 一种非均质压裂后储层流动数值模拟方法及系统 | |
US10393923B2 (en) | Method of determining a map of height of liquid hydrocarbon in a reservoir | |
Karlsson et al. | Comparison of methods for predicting regolith thickness in previously glaciated terrain, Stockholm, Sweden | |
Goloviznin et al. | Balance-characteristic scheme as applied to the shallow water equations over a rough bottom | |
Bazaikin et al. | Numerical analysis of topological characteristics of three-dimensional geological models of oil and gas fields | |
CN104375186A (zh) | 一种去除异常初至的方法及系统 | |
RU2757386C1 (ru) | Способ проведения электромагнитного мониторинга ГРП | |
CN109709602A (zh) | 一种远探测声波偏移成像方法、装置及系统 | |
Brown et al. | The evaluation of uncertainty in surfactant eor performance prediction | |
Oliveira et al. | On the generation of regular long waves in numerical wave flumes based on the particle finite element method | |
Mulder | Geophysical modelling of 3D electromagnetic diffusion with multigrid | |
Hashemi et al. | Two-point versus multiple-point geostatistics: the ability of geostatistical methods to capture complex geobodies and their facies associations—an application to a channelized carbonate reservoir, southwest Iran |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150311 |
|
RJ01 | Rejection of invention patent application after publication |