CN112597108A - 一种直流电测深反演数据处理系统和方法 - Google Patents
一种直流电测深反演数据处理系统和方法 Download PDFInfo
- Publication number
- CN112597108A CN112597108A CN202011378823.XA CN202011378823A CN112597108A CN 112597108 A CN112597108 A CN 112597108A CN 202011378823 A CN202011378823 A CN 202011378823A CN 112597108 A CN112597108 A CN 112597108A
- Authority
- CN
- China
- Prior art keywords
- inversion
- direct current
- data processing
- data
- storing
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/10—File systems; File servers
- G06F16/16—File or folder operations, e.g. details of user interfaces specifically adapted to file systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Human Computer Interaction (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明具体涉及一种直流电测深反演数据处理系统,包括存储直流电测深原始数据的EXCEL工作簿和基于VB语言开发的数据处理程序,所述数据处理程序包括写入指令、运算指令和存储指令,所述写入指令用于调用所述EXCEL工作簿存储的直流电测深原始数据,所述运算指令用于将所述直流电测深原始数据根据公式运算,所述存储指令用于将运算后的数据存储在所述EXCEL工作簿中。利用本发明提供的系统进行直流电测深反演数据处理,能够自动给出编制地电断面图所需要的深度‑电阻率转换数据。
Description
技术领域
本发明涉及地球物理勘探领域数据处理技术领域,特别是涉及一种直流电测深反演数据处理系统和方法。
背景技术
直流电测深方法原理是以地质体的电学特性为前提,探测不同地质体的电性差异,从而达到解决地质问题的一种方法。与地震探测法相比,它具有简单易学、耗资少,资料处理相对简单,短时间能取得地质成果的优点。适用于解决地震方法探测盲区段(0~500米)的断裂、地层中砂岩和泥岩展布状态等地质问题。
电测深资料反演法有许多种,Dar-Zarrouk反演法是其中的一种方法。Dar-Zarrouk反演法以层状模型的总横向电阻和总向电导为基础,根据Dar-Zarrouk曲线在电阻率随深度变化梯度不大情况下与视电阻率曲线基本重合这一结论推导出的一套直接反演公式,经后人改良,反演公式演变为以下形式:
(ρs)i为i个电极距所测的视电阻率;
grd(i)为第i个电极距处的梯度;
ΔH(i)为第i个极距所反应的地层的厚度;
H(i)为第i个极距所反应的地层的层底深度。
这六组公式计算简单,用Excel电子表格就能快速地将探测数据进行反演计算。但实际应用中发现存在如下问题:
1.如图1所示,由于公式(4)、(5)对探测数据有适用范围的要求,计算结果中出现大量的错误值(如#NUM、#DIV/0!等)。这些错误值夹在数据中,需要手工剔除。
2.由于反演后的数据在同一测线的不同测点之间其地层底界埋深(H)均不同,因此在垂面上不能组成规范的网格,无法用地质绘图软件绘制地层电性等值线剖面图,需手工用线性内插法将反演的真实电阻率换算为在对数坐标上大致等间距深度的真实电阻率。
发明内容
基于此,有必要针对上述问题,提供一种直流电测深反演数据处理系统和方法,能够方便、快捷的将野外探测的直流电测深原始数据转化为实际应用的地质成果。
为了实现上述目的,本发明提供如下技术方案:
一种直流电测深反演数据处理系统,包括存储直流电测深原始数据的EXCEL工作簿和基于VB语言开发的数据处理程序,所述数据处理程序包括写入指令、运算指令和存储指令,所述写入指令用于调用所述EXCEL工作簿存储的直流电测深原始数据,所述运算指令用于将所述直流电测深原始数据根据公式运算,所述存储指令用于将运算后的数据存储在所述EXCEL工作簿中。
进一步地,所述存储直流电测深原始数据的EXCEL工作簿包括存储直流电测深原始数据的工作表。
进一步地,所述存储直流电测深原始数据的工作表包含以下信息:线号、点号、点距、电极距和视电阻率,所述线号、点号、点距、电极距和视电阻率的名称分别位于所述工作表的第一行第一列、第一行第二列、第一行第三列、第一行第四列和第一行第五列,所述行号、点号、点距、电极距和视电阻率的数据值分别位于所述工作表的第N行第一列、第N行第二列、第N行第三列、第N行第四列和第N行第五列,所述N=2,3……t,t为自然数,其余列为空白列。
进一步地,所述写入指令包含以下信息:Excel工作簿的文件名,Excel工作簿的存储位置,以及以输入对话框的形式输入的所述存储直流电测深原始数据的工作表的信息。
进一步地,所述存储直流电测深原始数据的EXCEL工作簿的存储位置与所述写入指令中的Excel工作簿的存储位置一致,所述存储直流电测深原始数据的EXCEL工作簿的文件名需与写入指令的Excel工作簿的文件名一致,所述以输入对话框的形式输入的所述存储直流电测深原始数据的工作表的信息包括所述存储直流电测深原始数据的工作表的第一行行号、电极距所在列列号和电极距的个数。
进一步地,所述公式包括Dar-Zarrouk反演公式和内插公式。
进一步地,所述运算指令根据Dar-Zarrouk反演公式运算时,Dar-Zarrouk反演公式中的数据采用写入指令调用的所述EXCEL工作簿存储的直流电测深原始数据,运算后得到反演地层厚度、反演地层底界埋深和反演地层真实电阻率。
进一步地,所述写入指令调用的所述EXCEL工作簿的工作表中存储的直流电测深原始数据采用的是电极距≤3000米的直流电测深原始数据。
进一步地,所述运算指令根据内插公式运算时,内插公式中的数据采用是反演地层真实电阻率,运算后得到反演地层埋深和反演地层内插真实电阻率。
进一步地,所述运算指令能够自动按预先规范好的深度段格式将反演地层真实电阻率线性內插到各个深度段,实现反演地层真实电阻率的深度-电阻率变换,得到反演地层埋深和反演地层内插真实电阻率。
进一步地,所述深度段是在对数坐标上大致等间距的,包括:10、15、20、30、40、50、75、100、150、200、300、400、500、700、1000、1500、2000和2400米。
进一步地,所述运算指令进行反演地层真实电阻率的深度-电阻率变换时,会自动将每个缺失电阻率值的反演地层埋深做补零处理。
进一步地,所述存储指令中包含反演数据存储顺序,所述存储指令按照所述反演数据存储顺序将所述反演地层厚度、反演地层底界埋深、反演地层真实电阻率、反演地层埋深和反演地层内插真实电阻率存储在所述存储直流电测深原始数据的工作表中。
进一步地,所述存储指令将所述反演地层厚度、反演地层底界埋深、反演地层真实电阻率、反演地层埋深和反演地层内插真实电阻率分别存储在所述存储直流电测深原始数据的工作表的第六列、第七列、第八列、第十列和第十一列中。
一种基于上述直流电测深反演数据处理系统的直流电测深反演数据处理方法,包括如下步骤:
建立存储直流电测深原始数据的Excel工作簿,并将其存储在数据处理程序要求的存储位置上;
将数据处理程序的Excel工作簿的文件名修改为与所述Excel工作簿的文件名一致,然后将修改后的数据处理程序安装到VB的可执行文件上;
用VB的可执行文件启动数据处理程序,数据处理程序在执行过程中显示输入对话框,向输入对话框输入相应信息;
数据处理程序运算,并将运算完成后的数据存储在所述Excel工作簿中。
本发明的有益技术效果:
本发明提供的直流电测深反演数据处理系统和方法,能够自动给出编制地层电性等值线剖面图所需要的深度-电阻率转换数据,省略了人工大量内插计算的繁琐工作,节约了人力和时间,在最短的时间里实现了物探成果向地质成果的转化。
附图说明
图1为经excel处理后的数据样式图;
图2为经本发明的系统经Dar-Zarrouk反演公式运算后的数据样式图;
图3为经本发明的系统经Dar-Zarrouk反演公式和内插公式运算后的数据样式图;
图4为本发明的系统的Excel工作簿的工作表行、列示意图。
具体实施方式
以二连盆地赛汗高毕地区电法测量成果项目为例,详细说明本发明的直流电测深反演数据处理系统和方法。
一种直流电测深反演数据处理系统,包括存储直流电测深原始数据的EXCEL工作簿和基于VB语言开发的数据处理程序,所述数据处理程序包括写入指令、运算指令和存储指令,所述写入指令用于调用所述EXCEL工作簿存储的直流电测深原始数据,所述运算指令用于将所述直流电测深原始数据根据公式运算,所述存储指令用于将运算后的反演数据存储在所述EXCEL工作簿中。
进一步地,所述存储直流电测深原始数据的EXCEL工作簿包括存储直流电测深原始数据的工作表。
进一步地,所述存储直流电测深原始数据的工作表包含以下信息:线号(Lo)、点号(No)、点距(D)、电极距和视电阻率(ρs),所述线号(Lo)、点号(No)、点距(D)、电极距和视电阻率(ρs)的名称分别位于所述工作表的第一行第一列、第一行第二列、第一行第三列、第一行第四列和第一行第五列,所述行号(Lo)、点号(No)、点距(D)、电极距和视电阻率(ρs)的数据值分别位于所述工作表的第N行第一列、第N行第二列、第N行第三列、第N行第四列和第N行第五列,所述N=2,3……t,t为自然数,其余列为空白列。
进一步地,所述写入指令包含以下信息:Excel工作簿的文件名,Excel工作簿的存储位置,以及以输入对话框的形式输入的所述存储直流电测深原始数据的工作表的信息。
进一步地,所述存储直流电测深原始数据的EXCEL工作簿的存储位置与所述写入指令的Excel工作簿的存储位置一致,所述存储直流电测深原始数据的EXCEL工作簿的文件名与写入指令的Excel工作簿的文件名一致,所述以输入对话框的形式输入的所述存储直流电测深原始数据的工作表的信息包括所述存储直流电测深原始数据的工作表的第一行行号、电极距所在列列号和电极距的个数。
电极距的个数就是直流电测深原始数据中使用的电极距数量,本实施例使用了18个电极距,分别为:3、4.5、6、9、15、25、40、65、100、150、250、350、500、750、1000、1500、2000和3000米。
进一步地,所述公式包括Dar-Zarrouk反演公式和内插公式。
进一步地,所述运算指令根据Dar-Zarrouk反演公式运算时,Dar-Zarrouk反演公式中的数据采用写入指令调用的所述EXCEL工作簿存储的直流电测深原始数据,运算后得到反演地层厚度(△H)、反演地层底界埋深(H)和反演地层真实电阻率(ρt)。
进一步地,所述写入指令调用的所述EXCEL工作簿的工作表中存储的直流电测深原始数据采用的是电极距≤3000米的直流电测深原始数据。
野外生产一般都使用对称四极测量装置进行直流电测深测量,获得直流电测深原始数据。对称四极测量装置的标准电极距包括:3、4.5、6、9、15、25、40、65、100、150、250、350、500、750、1000、1500、2000、3000米等。
进一步地,所述运算指令根据内插公式运算时,内插公式中的数据采用是反演地层真实电阻率(ρt),运算后得到反演地层埋深(Hi)和反演地层内插真实电阻率(ρtr)。
进一步地,所述运算指令能够自动按预先规范好的深度段格式将反演地层真实电阻率(ρt)线性內插到各个深度段,实现反演地层真实电阻率(ρt)的深度-电阻率变换,得到反演地层埋深(Hi)和反演地层内插真实电阻率(ρtr)。
进一步地,所述深度段是在对数坐标上大致等间距的,包括:10、15、20、30、40、50、75、100、150、200、300、400、500、700、1000、1500、2000和2400米。
进一步地,所述运算指令进行反演地层真实电阻率(ρt)的深度-电阻率变换时,会自动将每个缺失电阻率值的反演地层埋深(Hi)做补零处理(当反演底界地层埋深(H)远远小于反演地层埋深(Hi)时,会出现底部数个反演地层埋深(Hi)缺失电阻率值现象)。
进一步地,所述存储指令中包含反演数据存储顺序,所述存储指令按照所述反演数据存储顺序将所述反演地层厚度(△H)、反演地层底界埋深(H)、反演地层真实电阻率(ρt)、反演地层埋深(Hi)和反演地层内插真实电阻率(ρtr)存储在所述存储直流电测深原始数据的工作表中。
进一步地,所述存储指令将所述反演地层厚度(△H)、反演地层底界埋深(H)、反演地层真实电阻率(ρt)、反演地层埋深(Hi)、反演地层内插真实电阻率(ρtr)分别存储在所述存储直流电测深原始数据的工作表的第六列、第七列、第八列、第十列和第十一列中。
本发明还提供基于上述直流电测深反演数据处理系统的直流电测深反演数据处理方法,包括如下步骤:
1、打开Microsoft Excel软件,新建一个Excel工作簿,命名所述Excel工作簿的文件名;
2、在所述Excel工作簿上设置一个Sheet1工作表,将直流电测深原始数据按照规定格式存储在所述Sheet1工作表中;
3、将所述Excel工作簿存储到数据处理程序要求的存储位置(计算机驱动器C盘)上;
4、将数据处理程序的Excel工作簿的文件名修改为步骤1的Excel工作簿的文件名,然后通过装入事件(Load)或点击事件(Click)将修改后的数据处理程序安装到VB的可执行文件(.exe)上;
6、数据处理程序运算;
7、数据处理程序将运算完成后的数据存储在Excel工作簿中。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种直流电测深反演数据处理系统,其特征在于,包括存储直流电测深原始数据的EXCEL工作簿和基于VB语言开发的数据处理程序,所述数据处理程序包括写入指令、运算指令和存储指令,所述写入指令用于调用所述EXCEL工作簿存储的直流电测深原始数据,所述运算指令用于将所述直流电测深原始数据根据公式运算,所述存储指令用于将运算后的数据存储在所述EXCEL工作簿中。
2.根据权利要求1所述的一种直流电测深反演数据处理系统,其特征在于,所述存储直流电测深原始数据的EXCEL工作簿包括存储直流电测深原始数据的工作表。
3.根据权利要求1所述的一种直流电测深反演数据处理系统,其特征在于,所述写入指令包含以下信息:Excel工作簿的文件名,Excel工作簿的存储位置,以及以输入对话框的形式输入的所述存储直流电测深原始数据的工作表的信息,所述存储直流电测深原始数据的EXCEL工作簿的存储位置与所述写入指令中的Excel工作簿的存储位置一致,所述存储直流电测深原始数据的EXCEL工作簿的文件名需与写入指令的Excel工作簿的文件名一致,所述以输入对话框的形式输入的所述存储直流电测深原始数据的工作表的信息包括所述存储直流电测深原始数据的工作表的第一行行号、电极距所在列列号和电极距的个数。
4.根据权利要求1所述的一种直流电测深反演数据处理系统,其特征在于,所述公式包括Dar-Zarrouk反演公式和内插公式。
5.根据权利要求5所述的一种直流电测深反演数据处理系统,其特征在于,所述运算指令根据Dar-Zarrouk反演公式运算时,Dar-Zarrouk反演公式中的数据采用写入指令调用的所述EXCEL工作簿存储的直流电测深原始数据,所述直流电测深原始数据电极距≤3000米,运算后得到反演地层厚度、反演地层底界埋深和反演地层真实电阻率。
6.根据权利要求5所述的一种直流电测深反演数据处理系统,其特征在于,所述运算指令根据内插公式运算时,内插公式中的数据采用是反演地层真实电阻率,运算后得到反演地层埋深和反演地层内插真实电阻率。
7.根据权利要求6所述的一种直流电测深反演数据处理系统,其特征在于,所述运算指令能够自动按预先规范好的深度段格式将反演地层真实电阻率线性內插到各个深度段,实现反演地层真实电阻率的深度-电阻率变换,得到反演地层埋深和反演地层内插真实电阻率。
8.根据权利要求7所述的一种直流电测深反演数据处理系统,其特征在于,所述运算指令进行反演地层真实电阻率的深度-电阻率变换时,会自动将每个缺失电阻率值的反演地层埋深做补零处理。
9.根据权利要求2所述的一种直流电测深反演数据处理系统,其特征在于,所述存储指令中包含反演数据存储顺序,所述存储指令按照所述反演数据存储顺序将运算后的数据存储在所述存储直流电测深原始数据的工作表中。
10.一种直流电测深反演数据处理方法,其特征在于,使用根据权利要求1-9任意一项所述的一种直流电测深反演数据处理系统,包括如下步骤:
建立存储直流电测深原始数据的Excel工作簿,并将其存储在数据处理程序要求的存储位置上;
将数据处理程序的Excel工作簿的文件名修改为与所述Excel工作簿的文件名一致,然后将修改后的数据处理程序安装到VB的可执行文件上;
用VB的可执行文件启动数据处理程序,数据处理程序在执行过程中显示输入对话框,向输入对话框输入相应信息;
数据处理程序运算,并将运算完成后的数据存储在所述Excel工作簿中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011378823.XA CN112597108A (zh) | 2020-11-30 | 2020-11-30 | 一种直流电测深反演数据处理系统和方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011378823.XA CN112597108A (zh) | 2020-11-30 | 2020-11-30 | 一种直流电测深反演数据处理系统和方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112597108A true CN112597108A (zh) | 2021-04-02 |
Family
ID=75187336
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011378823.XA Pending CN112597108A (zh) | 2020-11-30 | 2020-11-30 | 一种直流电测深反演数据处理系统和方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112597108A (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108107478A (zh) * | 2017-12-25 | 2018-06-01 | 湖南科技大学 | 大地电磁同步探测与实时反演方法及系统 |
CN108169801A (zh) * | 2018-01-16 | 2018-06-15 | 陕西铁道工程勘察有限公司 | 高分辨地电阻率快速成像方法 |
CN108535783A (zh) * | 2018-04-16 | 2018-09-14 | 中煤科工集团西安研究院有限公司 | 一种对电阻率断面进行层位校正的方法及装置 |
CN110333543A (zh) * | 2019-07-03 | 2019-10-15 | 山东大学 | 基于反射系数分析的低阻体解释及成像方法与系统 |
CN111305834A (zh) * | 2020-02-25 | 2020-06-19 | 北京工业大学 | 基于多探测模式电阻率测井的三维反演初始模型构建方法 |
-
2020
- 2020-11-30 CN CN202011378823.XA patent/CN112597108A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108107478A (zh) * | 2017-12-25 | 2018-06-01 | 湖南科技大学 | 大地电磁同步探测与实时反演方法及系统 |
CN108169801A (zh) * | 2018-01-16 | 2018-06-15 | 陕西铁道工程勘察有限公司 | 高分辨地电阻率快速成像方法 |
CN108535783A (zh) * | 2018-04-16 | 2018-09-14 | 中煤科工集团西安研究院有限公司 | 一种对电阻率断面进行层位校正的方法及装置 |
CN110333543A (zh) * | 2019-07-03 | 2019-10-15 | 山东大学 | 基于反射系数分析的低阻体解释及成像方法与系统 |
CN111305834A (zh) * | 2020-02-25 | 2020-06-19 | 北京工业大学 | 基于多探测模式电阻率测井的三维反演初始模型构建方法 |
Non-Patent Citations (4)
Title |
---|
喻忠鸿 等: "高密度电法二维反演程序反演常规电法电测深资料探讨", 《工程地球物理学报》 * |
喻忠鸿 等: "高密度电法二维反演程序反演常规电法电测深资料探讨", 《工程地球物理学报》, vol. 6, no. 06, 30 December 2009 (2009-12-30), pages 724 - 727 * |
熊彬 等: "MATLAB地球物理科学计算实战", 中国地质大学出版社, pages: 341 - 346 * |
陆柏树等: "基于数据库的直流电测深数据处理系统及其应用", 《工程地球物理学报》, no. 01, 15 February 2008 (2008-02-15), pages 21 - 26 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103975341B (zh) | 基于4d饱和度模型和仿真模型的储层建模 | |
EP3298437B1 (en) | Inversion for tectonic stress | |
Benesh et al. | Geometry, kinematics, and displacement characteristics of tear-fault systems: An example from the deep-water Niger Delta | |
CN105223617A (zh) | 利用地震属性控制初始波阻抗的建模方法及系统 | |
CN109799540B (zh) | 基于地质信息约束下的火山岩型铀矿床磁化率反演方法 | |
CN112150582B (zh) | 一种面向多模态数据的地质剖面图近似表达方法 | |
Maesano et al. | Vel-IO 3D: a tool for 3D velocity model construction, optimization and time-depth conversion in 3D geological modeling workflow | |
US20140163943A1 (en) | System for modeling geologic structures | |
CN113552621B (zh) | 页岩气地应力确定方法和装置 | |
NO20101734A1 (no) | Hastighetsmodeller for en enkeltbronner og for et sett av bronner | |
CN114542056A (zh) | 水平井地层解释模型快速构建方法和装置 | |
CN108957554B (zh) | 一种地球物理勘探中的地震反演方法 | |
Putz-Perrier et al. | Spatial distribution of brittle strain in layered sequences | |
EP3929630A1 (en) | Computer implemented method for correcting a reservoir model of a reservoir geological formation based on seismic images | |
CN112597108A (zh) | 一种直流电测深反演数据处理系统和方法 | |
CN108074282B (zh) | 基于三维地质建模软件的物性参数的厚度校正方法 | |
Miroslav et al. | Tutorial: The gravity-stripping process as applied to gravity interpretation in the eastern Mediterranean | |
CN112394392B (zh) | 对烃源岩分布情况进行评价的方法和装置 | |
CN111260791A (zh) | 一种更新地质导向模型的方法 | |
Wen et al. | Seismic structure-constrained inversion of CSAMT data for detecting karst caves | |
CN113496070B (zh) | 地层俘获截面曲线的处理方法、装置、设备及介质 | |
CN110927789B (zh) | 一种基于欠失数据预测页岩平面分布的方法及装置 | |
CN117238405B (zh) | 一种基于深度学习的地球化学数据分析方法和装置 | |
AU2015259067A1 (en) | Geomechanical modeling using dynamic boundary conditions from time-lapse data | |
Cabarcas | Pitfalls locating microseismic events from borehole measurements—Practical observations from field applications |
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 |