RU2399931C2 - Способ количественного разделения эффектов электромагнитной индукции и вызванной поляризации - Google Patents
Способ количественного разделения эффектов электромагнитной индукции и вызванной поляризации Download PDFInfo
- Publication number
- RU2399931C2 RU2399931C2 RU2008141889/28A RU2008141889A RU2399931C2 RU 2399931 C2 RU2399931 C2 RU 2399931C2 RU 2008141889/28 A RU2008141889/28 A RU 2008141889/28A RU 2008141889 A RU2008141889 A RU 2008141889A RU 2399931 C2 RU2399931 C2 RU 2399931C2
- Authority
- RU
- Russia
- Prior art keywords
- field
- electromagnetic induction
- model
- functions
- medium
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 230000005674 electromagnetic induction Effects 0.000 title claims abstract description 32
- 230000000694 effects Effects 0.000 title abstract description 9
- 238000009434 installation Methods 0.000 claims abstract description 5
- 230000010287 polarization Effects 0.000 claims description 49
- 230000015572 biosynthetic process Effects 0.000 claims description 31
- 238000005259 measurement Methods 0.000 claims description 19
- 238000000926 separation method Methods 0.000 claims description 10
- 238000012905 input function Methods 0.000 claims description 7
- 230000005672 electromagnetic field Effects 0.000 claims description 6
- 230000002123 temporal effect Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 19
- 239000000126 substance Substances 0.000 abstract 1
- 238000005755 formation reaction Methods 0.000 description 29
- 230000005520 electrodynamics Effects 0.000 description 9
- 230000006698 induction Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 230000001052 transient effect Effects 0.000 description 6
- 238000009826 distribution Methods 0.000 description 4
- 230000007704 transition Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 238000009827 uniform distribution Methods 0.000 description 2
- RLLPVAHGXHCWKJ-IEBWSBKVSA-N (3-phenoxyphenyl)methyl (1s,3s)-3-(2,2-dichloroethenyl)-2,2-dimethylcyclopropane-1-carboxylate Chemical compound CC1(C)[C@H](C=C(Cl)Cl)[C@@H]1C(=O)OCC1=CC=CC(OC=2C=CC=CC=2)=C1 RLLPVAHGXHCWKJ-IEBWSBKVSA-N 0.000 description 1
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 241000566515 Nedra Species 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000002500 effect on skin Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Images
Classifications
-
- 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
-
- 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/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
-
- 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/12—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation 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/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
- G01V2003/086—Processing
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Изобретение относится к области электроразведочных исследований. Сущность: измеряют процесс становления поля над поляризующейся средой дипольно-осевой установкой при пропускании импульсов тока. Формируют несколько функций так, что они по-разному зависят от полей электромагнитной индукции и вызванной поляризации. Одну из этих функций формируют так, чтобы повысить соотношение электромагнитной индукции и вызванной поляризации по сравнению с DU(t)=ΔU(t)/ΔU0, где ΔU0 - разность потенциалов ΔU, измеренная во время пропускания тока. Вторую из функций формируют так, чтобы понизить указанное соотношение по сравнению с DU(t). Третью из функций формируют как комбинацию временных и пространственных производных поля становления. Осуществляют инверсию одновременно для всех функций, включая DU(t), полученных в одной точке записи. Получают геоэлектрическую модель разреза среды. В полученной модели обнуляют поляризуемость для всех слоев и путем решения прямой задачи рассчитывают поле электромагнитной индукции. В той же модели обнуляют волновые числа и путем решения прямой задачи рассчитывают поле гальванической составляющей вызванной поляризации IP. Оценивают изменение гальванической составляющей по площади и осуществляют ее геологическую интерпретацию. Технический результат: возможность осуществить количественное определение интенсивности полей электромагнитной индукции и вызванной поляризации с погрешность не более 0,5% и обеспечить большой временной интервал разделения полей от 1 мс до нескольких секунд, в том числе и над разрезами с высокой проводимостью. 3 з.п. ф-лы, 13 ил.
Description
Изобретение относится к области электроразведочных исследований. Настоящее изобретение относится к способам определения свойств подземных формаций, при этом осуществляют измерения параметров электрических свойств формации с разделением измеренных параметров.
Изобретение используется как в наземной, так и в морской электроразведке с контролируемыми источниками электромагнитного поля и применяется в комплексе нефтегазопоисковых работ для поисков и оконтуривания нефтегазовых залежей на основе выделения отклика от вторично измененных под воздействием миграции углеводородов пород для разделения измеренного сигнала на какие-либо составляющие. Изобретение позволяет обеспечить совокупность технологических приемов, обеспечивающих послойное определение значения геоэлектрических параметров, а также характеристик процессов вызванной поляризации среды и электромагнитной индукции.
Известно изобретение «Способ геоэлектроразведки (варианты)». Патент RU 2235347, опубл. 2004.08.27, МПК G01V 3/06, в соответствии с которым возбуждают электромагнитное поле в толще исследуемой среды, пропуская через нее периодическую последовательность прямоугольных импульсов тока с паузами между каждым из них. В конце каждого импульса тока и между импульсами тока в паузах на всем протяжении измеряют мгновенные значения осевых и ортогональных разностей потенциалов, определяя разности потенциалов. Рассчитывают n нормированных параметров по формуле, решая обратную задачу на основе дифференциального волнового уравнения математической физики для напряженности дипольного источника в поляризующейся проводящей среде, и находят модель среды, наиболее близкую по геометрическому строению и электрическим параметрам к исследуемой. Строят временные разрезы этой модели по входящим в данное уравнение электрофизическим параметрам. Данный способ позволяет обнаружить, осуществить оконтуривание нефтегазовых залежей и оценку качества их насыщения, а также обеспечить возможность разделения параметров электропроводности и вызванной поляризации. Однако не позволяет решить главную проблему количественного раздельного изучения процессов электромагнитной индукции и вызванной поляризации, не учитывает эквивалентностей электрофизических параметров найденной модели, не определяет их информативности.
Известно изобретение «Способ морской геоэлектроразведки с фокусировкой электрического тока (варианты)», патент RU 2284555, опубл. 2006.09.27, 7 МПК G01V 3/06, в соответствии с которым возбуждают электромагнитное поле в толще исследуемой среды, измеряют мгновенные значения первых и вторых разностей электрических потенциалов, обеспечивают условие равенства нулю результирующей осевой разности электрических потенциалов вдоль профиля и определяют коэффициент из уравнения, рассчитывают три множества независимых нормированных электрических параметров, решают математическую обратную задачу, определяя три электрофизических параметра: удельную электропроводность, вызванную поляризацию и постоянную времени спада разности потенциалов вызванной поляризации и строят три временных разреза по этим параметрам. Изобретение заключается в обеспечении возможности разделения параметров электропроводности и вызванной поляризации, а также дополнительно дает возможность определения постоянной времени спада разности потенциалов вызванной поляризации - важного третьего наряду с двумя первыми параметра. Однако получают лишь суммарные сведения о всех элементах строения исследуемой среды, в которой развивается поле, так как в ней распределение в пространстве измеряемого тока I источника ничем не контролируется, и информации об указанном распределении в реально существующих трехмернонеоднородных средах нет. Изобретение не позволяет решить главную проблему количественного раздельного изучения процессов электромагнитной индукции и вызванной поляризации, не учитывает эквивалентностей электрофизических параметров найденной модели, не определяет их информативности.
Известно изобретение «Способ геоэлектроразведки (варианты)», патент RU 2231089 от 2003.07.08 г., опубл. 2004.06.20, МПК G01V 3/06, в соответствии с которым возбуждают электромагнитное поле в толще исследуемой среды, пропуская через нее периодическую последовательность прямоугольных импульсов тока с паузами после каждого из них, в точках наблюдения измеряют разности электрических потенциалов, используя значения нормированных параметров и дифференциальное уравнение математической физики для напряженности электрического поля дипольного источника в электрохимически поляризующейся проводящей среде, решают обратную задачу и определяют присущие каждому элементу среды электрофизические параметры: удельную электропроводность, вызванную поляризацию и постоянную времени спада разности потенциалов вызванной поляризации и строят временные разреза по этим параметрам. Изобретение позволяет обнаружить и оконтурить нефтегазовые залежи, а также произвести оценку качества их насыщения. Обеспечивается разделение параметров электропроводности и вызванной поляризации. Однако не удается до конца разделить присущую элементам геологической среды, вызванную поляризацию от переходных электродинамических процессов, связанных с электропроводностью указанных элементов толщи слагающих разрез геологических пород. Не позволяет решить главную проблему количественного раздельного изучения полей электромагнитной индукции и вызванной поляризации, не учитывает эквивалентностей электрофизических параметров найденной модели, не определяет их информативности.
Известно изобретение «Способ анализа геоэлектроразведочного сигнала», патент RU 2253137, опубл. 2005.01.10, 7 МПК G01V 3/38, включающий преобразование каждого измеряемого сигнала одновременно по крайней мере двумя модуляторами, после чего находят интервалы с наименьшим шагом. Изобретение позволяет повысить точность и информативность за счет одновременного анализа сигнала (регулярной компоненты) и помех (случайной компоненты). Однако целью изобретения является использование данного метода, сравнивают выходные сигналы интеграторов дельта-модуляторов с входным сигналом соответствующего канала на стробируемом компараторе. Однако предложенный способ анализа можно осуществить только при применении, по крайней мере, двух дельта-модуляторов с существенно различающимися шагами квантования, системы приемных электродов, датчиков магнитного поля либо любых иных датчиков, обладающих достаточной широкополостностью и пространственной угловой селективностью (направленностью). Данный способ невозможно применить для анализа параметров, полученных при использовании процессов вызванной поляризации для поисков и разведки месторождений нефти и газа.
Известно изобретение «Способ определения электрического удельного сопротивления геологической формации, окружающей ствол скважины, заполненный скважинным флюидом», патент RU 2209451, опубл. 2003.07.27, 7 МПК G01V 3/38, включающий получение некоторой совокупности диаграмм, выбора моделированного профиля, введения дополнительно моделированных профилей, повторения операций корректировки моделированного профиля до тех пор, пока разница между диаграммой и соответствующей моделированной диаграммой не будет ниже выбранного порогового значения, а также включает в себя корректировку каждого моделированного профиля как функцию отношения величин. Изобретение позволяет обеспечить усовершенствования способа определения электрического удельного сопротивления геологической формации, в котором учитывается вторжение скважинного флюида в формацию, окружающую ствол скважины. Однако изобретение не позволяет обнаружить, осуществить оконтуривание нефтегазовых залежей и оценку качества их насыщения, а также обеспечить возможность разделения параметров электропроводности и вызванной поляризации. Кроме того, не позволяет решить главную проблему количественного раздельного изучения процессов электромагнитной индукции и вызванной поляризации, не учитывает эквивалентностей электрофизических параметров найденной модели, не определяет их информативности.
Наиболее близким по способу обработки геоэлектрических характеристик является изобретение «Способ определения свойств подземных формаций», патент RU 2294547, опубл. 2007.02.27, 7 МПК G01V 3/38, включающий разработку модели параметров формации, моделирование измерений электрических свойств, сравнение моделированных измерений с измеренными электрическими свойствами, определение физических/нефтефизических свойств формации на основании различных измерений для обеспечения комбинированной инверсии, оценки физических/нефтефизических свойств формации при помощи модели и различных измерений электрических свойств, а также использование комплексной модели. Изобретение позволяет обеспечить реализацию усовершенствованных способов обработки данных при осуществлении разреза буровой скважины с применением электричества, а также дать более точное определение величины Rt с целью оценки характеристик формации с использованием определения залегания формации с введением в схему инверсии первичной информации - геометрическая информация для стабилизации неточно поставленных задач «расширения возможностей». Однако изобретение не позволяет обнаружить, осуществить оконтуривание нефтегазовых залежей и оценку качества их насыщения, а также обеспечить возможность разделения параметров электропроводности и вызванной поляризации. Кроме того, не позволяет решить главную проблему количественного раздельного изучения процессов электромагнитной индукции и вызванной поляризации, не учитывает эквивалентностей электрофизических параметров найденной модели, не определяет их информативности. Изобретение имеет другую цель, осуществить метод, когда модель среды просто уточняется посредством итерационного использования инверсии и не происходит разделения измеренного сигнала на какие-либо составляющие. В нашем же случае найденная в результате решения обратной задачи модель среды служит для разделения измеренного сигнала ΔU(t) на составляющие EM(t) и Ipg(t).
При пропускании электрического тока в среде проявляются эффекты, которые определены различными физическими явлениями. При этом выделяется электромагнитная индукция, которая характеризуется индукционными полями, что определено волновой теорией, и
вызванная поляризация, которая определена электрохимическими процессами в среде.
До настоящего времени задача разделения этих эффектов решалась путем увеличения времени регистрации переходных процессов от импульса возбуждения в среде, вызывающего электромагнитную индукцию и поляризацию. При этом с увеличением времени до момента измерения поле электромагнитной индукции успевало уменьшаться до пренебрежительно малых величин. Однако с уменьшением времени до момента измерения, в начальные моменты времени после импульса возбуждения в среде, эффект поляризации пренебрежимо мал не будет. При измерении проводимости платформ, имеющих большую мощность осадочного чехла, электромагнитные колебания затухают медленно, вследствие большой (сотни и тысячи Сименс) проводимости осадочного чехла. Наличие данной проблемы и отсутствие ее решения указано в книге Уэйта (Уэйт Дж. Геоэлектромагнетизм. М., Недра, 1987. - 235 с.).
Техническим результатом предложенного способа является то, что он позволяет решить главную проблему количественного раздельного изучения полей электромагнитной индукции и вызванной поляризации и учесть эквивалентность электрофизических параметров найденной модели, а также определить их информативность.
Данный способ позволяет
- осуществить количественное определение интенсивности полей электромагнитной индукции и вызванной поляризации постоянного тока с погрешность не более 0,5%,
- обеспечить большой временной интервал разделения полей от 1 мс до нескольких секунд, в том числе и над разрезами с высокой проводимостью.
Данный способ можно применить, поскольку функция ΔU является всегда аналитична, т.е. всегда непрерывна, а также непрерывны и ее производные 2-го порядка. Кроме того, поля электромагнитной индукции и вызванной поляризации - не аддитивны, т.к. часть поля вызванной поляризации (ВП) вызвана токами индукции, и наоборот - часть полей индукции обусловлена полями ВП. Мы будем говорить только о ВП, созданной постоянным током за время его пропускания в среду.
Заявленный технический результат обеспечен тем, что способ количественного разделения эффектов электромагнитной индукции и вызванной поляризации включает:
- измерение процесса становления над поляризующейся средой;
- моделирование электромагнитных измерений в рамках горизонтально-слоистой поляризующейся среды;
- выделение электродинамической составляющей из процесса становления;
- выделение гальванической составляющей вызванной поляризации из процесса становления;
- оценки изменения гальванической составляющей вызванной поляризации по площади и ее геологическая интерпретация.
Способ отличается тем, что регистрируют несколько входных функций, зависящих от параметров среды, и измеренных на поверхности среды, формируют эти функции так, что они по-разному зависят от полей электромагнитной индукции и вызванной поляризации, - при этом одна из этих функций, например Ps, рассчитывается так, чтобы повысить соотношение электромагнитной индукции и вызванной поляризации по сравнению с DU(t), при этом одна из этих функций, например P1, рассчитывается так, чтобы понизить соотношение электромагнитной индукции и вызванной поляризации по сравнению с DU(t), и одна из этих функций, например Dϕ, является комбинацией временных и пространственных производных поля становления. Далее осуществляют инверсию одновременно для всех входных функций и получают геоэлектрическую модель разреза, при этом число параметров среды, характеризующих вызванную поляризацию, является заведомо избыточным, так что найденная модель среды не является единственной, а лежит в некоторой области эквивалентности. В полученной модели обнуляют поляризуемость для всех слоев и путем решения прямой задачи рассчитывают поле электромагнитной индукции ЕМ, в той же модели обнуляют волновые числа и путем решения прямой задачи рассчитывают поле гальванической составляющей вызванной поляризации IP. Кроме того, например, формируют параметр P1, принимая его за фильтр, повышающий, в сравнении с ΔU(t), соотношение ВП/Электродинамика, и определяют P1 как отношение второй осевой разности потенциалов поля становления Δ2U(t) к первой осевой разности потенциалов поля становления ΔU(t) в режиме выключения тока по формуле P1(t)=Δ2U(t)/ΔU(t)/time off; или, например, формируют параметр Ps, принимая его за фильтр, понижающий, в сравнении с ΔU(t), соотношение ВП/Электродинамика, и определяют Ps как отношение 2-ой осевой разности потенциалов поля становления Δ2U(t) к первой осевой разности потенциалов ΔU(t) в режиме включения тока по формуле Ps(t)=Δ2U(t)/ΔU(t)/time on; в частности, формируют параметр Dφ(t) как разность параметра It (t) и параметра P1(t) в соответствии с формулой Dф(t)=It(t)-P1(t), при этом It (t) рассчитывают по формуле It(t)=Δt(Δ2U(t))/Δt(ΔU(t)), как отношение производных по времени Δ2U(t) и ΔU(t).
После формирования входных параметров, полученных в одной точке записи, осуществляют их минимизацию, для чего подают на вход программы минимизации параметры
DU(t)=ΔU(t)/ΔU0, где ΔU0 - разность потенциалов ΔU, измеренная во время пропускания тока,
- P1(t)=Δ2U(t)/ΔU(t), который рассчитан в режиме выключения тока,
- Ps(t)=Δ2U(t)/ΔU(t), который рассчитан в режиме включения тока,
- Dф(t)=It(t)-P1(t), где It(t)=Δt(Δ2U(t))/Δt(ΔU(t))
и осуществляют их инверсию в рамках горизонтально-слоистой поляризующейся среды при заведомо избыточном числе геоэлектрических параметров модели, характеризующих ее поляризационные свойства. При этом, как частный случай, определяют поле IP как результат решения прямой задачи при обнуленных волновых числах и определяют поле ЕМ как результат решения прямой задачи при обнуленной поляризуемости.
После совершения всех операций осуществляют проверку воспроизводимости решения.
Для решения поставленной задачи разделения эффектов следует выполнить несколько этапов.
I этап
Регистрируют несколько входных функций, которые можно измерить и построить на поверхности среды.
Функции конструируют так, что они по-разному зависят от электромагнитной индукции и вызванной поляризации. Одна из них повышает соотношение электромагнитной индукции к вызванной поляризации по сравнению с DU - кривой спада напряжения, другая - понижает. Данные функции являются комбинацией поля переходных процессов и его пространственных производных. Кроме того, одна из функций является комбинацией временных и пространственных производных.
II этап
1. Осуществляют инверсию (решают обратную задачу) одновременно для всех входных функций и получают некоторую модель разреза. При этом послойно определяются параметры среды: удельное электрическое сопротивление ρ, коэффициент поляризуемости η, время релаксации τ, показатель степени с. Коэффициент поляризуемости, время релаксации и показатель степени определяются в соответствии с моделью Cole-Cole. Параметризация модели по этим параметрам задается заведомо избыточной, т.е. мы получаем некоторое решение из класса эквивалентных.
2. В полученной модели зануляют поляризуемость для всех слоев и путем решения прямой задачи определяют поле электромагнитной индукции ЕМ.
3. В той же модели зануляют волновые числа и путем решения прямой задачи определяют поле вызванной поляризации постоянного тока IP.
Повторение последовательности действий (2-го этапа) показывает совпадение рассчитанных полей ЕМ и IP с точностью решения обратной задачи. Полученные модели будут различаться, но поля совпадать.
Причем поле IP - это поле, полученное при решении прямой задачи над данной моделью при задании волновых чисел в любой точке разреза, равными нулю, а поле ЕМ - это поле, полученное при решении прямой задачи над данной моделью при задании поляризуемости в любой точке разреза, равной нулю.
Способ иллюстрируется и подтверждается схемами и диаграммами, показанными на чертежах.
На Фиг.1 показаны графики ΔU1 и ΔU2 - разности потенциалов на приемных линиях ОМ и ON для дипольно-осевой установки.
Она иллюстрирует то, что существует бесконечное множество различных сочетаний полей IP и ЕМ. Для синтетической модели была просчитана кривая спада DU и дважды осуществлена ее инверсия с разных нулевых приближений с точностью до 0.1%. Получены два различных эквивалентных (по параметру DU) решения, каждому из которых соответствуют существенно различные значения полей IP и ЕМ.
На Фиг.2 показан интервал распределения токов ЕМ и IP в различных стадиях (фазах) переходного процесса: ранней - А), промежуточной (переходной) - В) и поздней - С).
При этом поле ЕМ, подчиняясь скин-эффекту, диффундирует в нижнее полупространство во время переходного процесса, а вихревые токи, вызывающие его, стремятся к равномерному распределению в нижнем полупространстве, поле IP, напротив, всегда остается пространственно неоднородным и уменьшается с удалением от источника тока - излучателя.
Исходя из смысла поставленной задачи, требуется ввести в рассмотрение такие трансформанты поля становления, которые были бы записаны в одной точке и имели различную функциональную зависимость от полей ВП и индукционных. Наиболее естественное решение - привлечение, в дополнение к самим полям становления, их нормированных пространственных производных (точнее - конечных разностей), поскольку пространственно-временная структура индукционных полей и полей ВП различна.
На Фиг.3 показаны - а) схема измерительной установки для регистрации параметров DU(t), P1(t), Ps(t), Dφ(t) во время движения судна; б) форма тока в генераторной линии АВ, для возбуждения переходного процесса; в) измеряемые сигналы Δ2U(t) и ΔU(t) на приемной линии.
На Фиг.4 показано затухание кривой DU(t) над земной поверхностью во времени при наличии поляризации (В) и при ее отсутствии (А), на приемной линии МО (на расстоянии 1000 м (1)) и на приемной линии ON (на расстоянии 2400 м (2)).
На Фиг.5 показана кривая параметра P1(t) над земной поверхностью при наличии поляризации (2) и при ее отсутствии (1). При η=1% и ρ=50 Om·m и η=0%.
На Фиг.6 показана кривая параметра Ps над земной поверхностью относительно, при η=4% и η=0, где η - коэффициент поляризуемости в терминах модели Cole-Cole.
На Фиг.7 показана кривая Dφ, рассчитанная при η=0 и при η=1%, где η - коэффициент поляризуемости в терминах модели Cole-Cole. Кривые строят относительно нулевой линии Dφ. Использование временных производных, наряду с пространственными характеристиками, позволяет получить более детальную характеристику переходных процессов.
На Фиг.8 показана трехмерная референтная модель, от которой были рассчитаны синтетические кривые параметров, соответствующие получаемым при полевых измерениях (DU(t), P1(t), Ps(t), Dφ(t)), и значения полей ЕМ и IP. Также приведена эквивалентная ей одномерная модель.
на Фиг.9 показаны поля ЕМ и IP, которые для приведенных на фиг.8 моделей совпали с точностью 0.15%.
на Фиг.10 показаны кривые, наблюденные в полевых измерениях на Северо-Гуляевской площади (шельф Баренцева моря) и соответствующие им модельные. Приведена одна из эквивалентных геоэлектрическая моделей.
на Фиг.11 показаны кривые, наблюденные в полевых измерениях на Северо-Гуляевской площади (шельф Баренцева моря) и соответствующие им модельные. Приведена одна из эквивалентных геоэлектрическая моделей, отличная от показанной на фиг.10.
на Фиг.12 - кривые спада полей ЕМ и IP, полученные от эквивалентных моделей отличаются друг от друга на величину, меньшую погрешности измерений, что составляет 0.5%. При этом поле IP выделяется на уровне поля ЕМ, превосходящего его на 2 порядка, в широком временном диапазоне.
на Фиг.13 - изолинии поля IР по профилю над известным нефтегазовым месторождением и выделение нефтегазовой залежи по полю IP.
Способ осуществляется следующим образом.
Имеется набор входных данных, полученных в результате измерений на земной поверхности.
Данные сформированы таким образом, что функциональная зависимость вызванной поляризации и электромагнитной индукции у всех у них разная. На основе этих данных определяется путем инверсии геоэлектрическая модель при условии заведомой избыточности числа параметров разреза, характеризующих процессы вызванной поляризации.
Если для этой модели определить поле IР как результат решения прямой задачи при зануленных волновых числах, а поле ЕМ как результат решения прямой задачи при зануленной поляризуемости, то при этом соблюдается условие воспроизводимости полей IP и ЕМ. Это означает, что расхождения, полученные при расчетах полей IР и ЕМ над различными эквивалентными моделями, не превышают реальной погрешности наблюдений (около 0.5%).
Для реализации способа
1. формируют, при измерении в режиме выключения токового импульса, параметр Р1 (фиг.1) - принимая его за фильтр, повышающий в сравнении с ΔU(t) соотношение ВП/Электродинамика
- для этого определяют Р1 как отношение второй осевой разности потенциалов поля становления Δ2U к первой осевой разности потенциалов поля становления ΔU в режиме выключения тока,
т.е. P1(t)=Δ2U(t)/ΔU(t).
На фиг.1 видны ΔU1 и ΔU2 - разности потенциалов на приемных линиях ОМ и ON для дипольно-осевой установки. Параметр Р1 реагирует на пространственную неоднородность поля переходных процессов. Поле электромагнитных токов описывается уравнением диффузии, поэтому при увеличении времени спада оно стремится к равномерному распределению в среде. Поле же ВП пропорционально плотности поляризующего тока, поэтому оно всегда пространственно неоднородно (Фиг.2) [Легейдо П.Ю. Теория и технологии дифференциально-нормированной геоэлектроразведки для изучения поляризующихся разрезов в нефтегазовой геофизике. Диссертация на соискание ученой степени доктора геолого-минералогических наук, 1998]. В силу этого параметр P1 можно рассматривать как фильтр, повышающий, в сравнении с полем становления ΔU, соотношение ВП/Электродинамика. (см. Фиг.4 и 5).
2. Формируют, при измерении в режиме включения токового импульса, параметр Ps - принимая его за фильтр, понижающий в сравнении с ΔU(t) соотношение ВП/Электродинамика
- для этого определяют отношение 2-ой осевой разности потенциалов Δ2U поля становления к 1-ой осевой разности потенциалов ΔU в режиме включения тока.
В ранней стадии становления соотношение ВП/Электродинамика мало, в поздней же стадии параметр Ps в основном определяется распределением сопротивлений, поляризуемость на него влияет мало. Поэтому параметр Ps можно рассматривать как фильтр, понижающий, в сравнении с полем становления ΔU, соотношение ВП/Электродинамика (Фиг.6).
3. Формируют параметр Dφ - для более детальной характеристики переходного процесса, т.к. он рассчитывается на основе временных производных Δ2U(t) и ΔU(t)
- получают как разность параметра It и параметра Р1 относительно условных единиц, являющихся логарифмом времени переходного процесса, т.е. Dϕ=It-P1
- при этом It рассчитывают как lt=Δt(Δ2U)/ Δt(ΔU), те. отношение производных по времени Δ2U(t) и ΔU(t).
Скорость спада поля ВП и электромагнитной индукции обычно различается, поэтому целесообразно использовать временные производные для более детальной характеристики переходных процессов (Фиг.7).
4. Осуществляют минимизацию параметров, полученных в одной точке записи
- для чего подают на вход программы минимизации параметры (Фиг.3)
- - DU(t)=ΔU(t)/ ΔU0, где ΔU0 - разность потенциалов ΔU, измеренная во время пропускания тока,
- - P1(t)=Δ2U/ΔU(t) рассчитан при режиме выключенного тока,
- - Ps(t)=Δ2U/ΔU(t) рассчитан при режиме включенного тока,
- - Dϕ(t)=It-P1, где It=Δt(Δ2U)/Δt(ΔU),
- осуществляют их инверсию в рамках горизонтально-слоистой поляризующейся среды при заведомо избыточном числе геоэлектрических параметров модели, чтобы как можно точнее аппроксимировать полевые кривые, записанные над неоднородной анизотропной и градиентной средой.
5. Осуществляют расчет полей IP и ЕМ по полученной модели среды
- для полученной модели зануляют поляризуемость во всех областях среды и путем решения прямой задачи получают поле электромагнитной индукции ЕМ,
- для полученной модели зануляют волновые числа во всех областях среды и путем решения прямой задачи получают поле вызванной поляризации IP.
6. Осуществляют проверку воспроизводимости решения
- проверка воспроизводимости полей ЕМ и IP показывает, что для практически всех типов разрезов они воспроизводятся с погрешностью, не превышающей погрешность измерений (т.е. около 0.5%) (см. Фиг 10 и 13).
Иллюстрацию способа разделения индукционных и поляризационных полей можно продемонстрировать на следующем примере.
Для трехмерной модели были получены синтетические кривые параметров, соответствующие получаемым при полевых измерениях, и значения полей ЕМ и IP (Фиг.8).
Для синтетических кривых была найдена эквивалентная модель в рамках горизонтально-слоистой поляризующейся среды, на основе которой были рассчитаны поля ЕМ и IP (Фиг.8)
Было произведено сравнение полей, полученных от 3-мерной и горизонтально-слоистой моделей (Фиг.9). Погрешность между ними (0.15%) не превышает погрешность расхождения между синтетическими от 3-мерной модели и модельными кривыми от горизонтально-слоистой поляризующейся среды. В силу малости расхождения различие графиков визуально незаметно.
На практике способ количественного разделения эффектов электромагнитной индукции и вызванной поляризации опробовался в самых различных геолого-геофизических условиях.
Над известным нефтегазовым месторождением на Северо-Гуляевской площади (шельф Баренцева моря) были выполнены электроразведочные работы с измерением параметров DU(t), P1(t), Ps(t), Dφ(t). Для полевых кривых были найдены горизонтально-слоистые поляризующиеся модели, с погрешностью аппроксимации менее 0.5% (Фиг.11 и 12). При этом для одного и того же пикета измерений было найдено несколько эквивалентных моделей и рассчитанные графики полей ЕМ и IP имели расхождение менее 0.5% (Фиг.12). При этом удалось выделить поле IP, которое на два порядка меньше поля ЕМ. В результате построения изолиний поля IP четко выделяется нефтегазовая залежь (Фиг.13),
причем Ro - сопротивление, Ом·м
С - показатель степени в модели Cole-Cole
Н - мощность слоя, м
η - коэффициент поляризуемости модели Cole-Cole, %
t - время релаксации модели Cole-Cole, с
Например, для модели с параметрами:
Ro=10 Омм, Eta=5%, Tau=0.5 с, С=0.5, Н=100 m
Ro=50 Омм, Eta=2%, Tau=0.5 с, С=0.5, Н=100 m
Ro=20 Омм, Eta=2%, Tau=1 с, С=0.5, Н=200 m
Ro=5 Омм, Eta=5%, Tau=1 с, С=0.5, Н=300 m
Ro=200 Омм, Eta=7%, Tau=1 с, С=0.5
Ro=500 Омм, Eta=0, H=450 m
Ro=5 Омм, Eta=0, H=800 m
Ro=1000 Омм, Eta=0, H=∞
Для этой модели были рассчитаны значения входных данных, соответствующих получаемым при полевых измерениях, а также полей IР и ЕМ. И получена для шельф Баренцева моря. Северо-Гуляевское месторождение диаграмма - см. Фиг.13.
Иными словами, исходя из смысла поставленной задачи, осуществили введение в рассмотрение такие трансформанты поля становления, которые были бы записаны в одной точке и имели различную функциональную зависимость от полей ВП и индукционных полей. При этом осуществили привлечение, в дополнение к самим полям становления, их нормированных пространственных производных (точнее - конечных разностей), поскольку пространственно-временная структура индукционных полей и полей ВП различна. Тем самым был получен заявленный технический результат.
Claims (4)
1. Способ количественного разделения электромагнитной индукции и вызванной поляризации, включающий измерение процесса становления поля над поляризующейся средой дипольно-осевой установкой при пропускании импульсов тока, получение модели разреза решением обратной задачи путем инверсии функций, полученных на основании измерения и зависящих от параметров среды, в рамках горизонтально-слоистой поляризующейся среды, при этом формируют несколько функций так, что они по-разному зависят от полей электромагнитной индукции и вызванной поляризации, причем одну из этих функций формируют так, чтобы повысить соотношение электромагнитной индукции и вызванной поляризации по сравнению с DU(t)=ΔU(t)/ΔU0, где ΔU0 - разность потенциалов ΔU, измеренная во время пропускания тока, вторую из этих функций формируют так, чтобы понизить соотношение электромагнитной индукции и вызванной поляризации по сравнению с DU(t), и третью из этих функций формируют как комбинацию временных и пространственных производных поля становления, после чего осуществляют инверсию одновременно для всех функций, включая DU(t), полученных в одной точке записи, и получают геоэлектрическую модель разреза среды, при этом число параметров среды, характеризующих вызванную поляризацию, является заведомо избыточным, так что найденная модель не является единственной, а лежит в некоторой области эквивалентности, в полученной модели обнуляют поляризуемость для всех слоев и путем решения прямой задачи рассчитывают поле электромагнитной индукции ЕМ, в той же модели обнуляют волновые числа и путем решения прямой задачи рассчитывают поле гальванической составляющей вызванной поляризации IP, оценивают изменение гальванической составляющей вызванной поляризации по площади и осуществляют ее геологическую интерпретацию.
2. Способ по п.1, отличающийся тем, что формируют входную функцию P1(t) как отношение второй осевой разности потенциалов Δ2U(t) к первой осевой разности потенциалов ΔU(t) в режиме выключения тока по формуле P1(t)=Δ2U(t)/ΔU(t)time off.
3. Способ по п.1, отличающийся тем, что формируют входную функцию Ps(t) как отношение второй осевой разности потенциалов Δ2U(t) к первой осевой разности потенциалов ΔU(t) в режиме включения тока по формуле Ps(t)=Δ2U(t)/ΔU(t)Itime on.
4. Способ по п.1, отличающийся тем, что формируют входную функцию Dφ(t) как разность It(t) и P1(t) в соответствии с формулой Dφ(t)= It (t)-P1(t), при этом It (t) рассчитывают по формуле It(t)=Δt(Δ2U(t))/Δt(ΔU(t)) как отношение производных по времени Δ2U(t) и ΔU(t).
Priority Applications (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2008141889/28A RU2399931C2 (ru) | 2008-10-23 | 2008-10-23 | Способ количественного разделения эффектов электромагнитной индукции и вызванной поляризации |
| MX2011002163A MX2011002163A (es) | 2008-10-23 | 2008-11-17 | Metodo para la separacion cuantitativa de los efectos de induccion electromagnetica y polarizacion inducida. |
| EP08877586.1A EP2352041A4 (en) | 2008-10-23 | 2008-11-17 | Method for quantitatively separating the effects of electromagnetic induction and induced polarization |
| PCT/RU2008/000711 WO2010047613A1 (ru) | 2008-10-23 | 2008-11-17 | Способ количественного разделения эффектов электромагнитной индукции и вызванной поляризации |
| CN200880132434.XA CN102265188B (zh) | 2008-10-23 | 2008-11-17 | 电磁感应和诱导极化效应的量化分解方法 |
| US12/588,507 US8138763B2 (en) | 2008-10-23 | 2009-10-19 | Method for quantitative separation of electromagnetic induction and induced polarization effects |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2008141889/28A RU2399931C2 (ru) | 2008-10-23 | 2008-10-23 | Способ количественного разделения эффектов электромагнитной индукции и вызванной поляризации |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| RU2008141889A RU2008141889A (ru) | 2010-04-27 |
| RU2399931C2 true RU2399931C2 (ru) | 2010-09-20 |
Family
ID=42119492
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2008141889/28A RU2399931C2 (ru) | 2008-10-23 | 2008-10-23 | Способ количественного разделения эффектов электромагнитной индукции и вызванной поляризации |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US8138763B2 (ru) |
| EP (1) | EP2352041A4 (ru) |
| CN (1) | CN102265188B (ru) |
| MX (1) | MX2011002163A (ru) |
| RU (1) | RU2399931C2 (ru) |
| WO (1) | WO2010047613A1 (ru) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2491580C1 (ru) * | 2012-02-07 | 2013-08-27 | Общество с ограниченной ответственностью "Сибирская геофизическая научно-производственная компания" | Способ измерения геофизических характеристик с применением последующей инверсии геоэлектрических данных с дополнительным временным фильтром |
| RU2494419C1 (ru) * | 2012-03-11 | 2013-09-27 | Открытое акционерное общество Научно-производственное предприятие "Научно-исследовательский и проектно-конструкторский институт геофизических исследований геологоразведочных скважин" (ОАО НПП "ВНИИГИС") | Способ геоэлектроразведки |
| RU2595327C1 (ru) * | 2015-03-03 | 2016-08-27 | Общество с ограниченной ответственностью "ЛУКОЙЛ-Приморьенефтегаз" | Устройство сейсморазведки 2d или 3d, электроразведки и гис для повышения надежности картирования кровли соли и для прогноза нефтегазоносности подсолевых отложений в районах с развитой солянокупольной тектоникой |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120099396A1 (en) * | 2010-10-22 | 2012-04-26 | Chevron U.S.A. Inc. | System and method for characterization with non-unique solutions of anisotropic velocities |
| EA025583B1 (ru) * | 2012-09-05 | 2017-01-30 | Сергей Иванович Говорков | Способ геоэлектроразведки и устройство для его осуществления |
| BR112017023431B1 (pt) * | 2015-04-30 | 2021-12-28 | Schlumberger Technology B.V. | Método para estimar saturação de água em medições eletromagnéticas |
| CN105068130B (zh) * | 2015-07-22 | 2018-02-16 | 中石化石油工程地球物理有限公司江汉分公司 | 全频段电性参数的统一数学表征形式及参数的计算方法 |
| CN105549101B (zh) * | 2016-01-28 | 2018-01-05 | 中国矿业大学 | 一种瞬变电磁数据微分电导解释方法 |
| CN108241781B (zh) * | 2018-01-04 | 2021-07-30 | 国家电网公司 | 一种实验室标准电场发生器的参数扫描优化设计方法 |
| CN110058316B (zh) * | 2019-05-10 | 2020-09-11 | 成都理工大学 | 一种基于电阻率等值性原理的电磁测深约束反演方法 |
| CN110673218A (zh) * | 2019-09-05 | 2020-01-10 | 中国科学院地质与地球物理研究所 | 一种接地导线源瞬变电磁响应中ip信息的提取方法 |
| CN114791633B (zh) * | 2022-06-21 | 2022-09-30 | 西安石油大学 | 监测页岩气压裂的方法、系统及介质 |
| CN115829001B (zh) * | 2022-11-08 | 2023-06-20 | 中国科学院地质与地球物理研究所 | 一种瞬变电磁-激电场分离及多参数信息提取方法和系统 |
| CN116699724B (zh) * | 2023-08-01 | 2023-10-10 | 北京建工环境修复股份有限公司 | 一种时间域激发极化数据质量评价方法、体系及系统 |
| CN118068445B (zh) * | 2024-04-18 | 2024-07-09 | 山东大学 | 一种时间域激发极化法极化率参数标准化方法及系统 |
| CN118348598B (zh) * | 2024-06-18 | 2024-08-09 | 山东大学 | 一种隧道超前地质预报全时域激发极化探测系统及方法 |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU1122998A1 (ru) * | 1983-06-03 | 1984-11-07 | Всесоюзное Морское Научно-Производственное Объединение Инженерной Геологии "Союзморинжгеология" | Способ морской геоэлектроразведки |
| WO1999018454A1 (en) | 1997-10-08 | 1999-04-15 | Shell Internationale Research Maatschappij B.V. | Resistivity log correction method |
| DE10006083B4 (de) * | 2000-02-11 | 2004-01-22 | INSTITUT FüR MIKROTECHNIK MAINZ GMBH | Verfahren zur quantitativen und/oder qualitativen Bestimmung von Schichtdicken sowie ein Mikroreaktionsgefäß und eine Titerplatte |
| CA2299743A1 (en) * | 2000-02-22 | 2001-08-22 | Kenneth Duckworth | A method and apparatus for the detection of induced polarization and electrical resistivity in rocks and soils free of the effect of the electromagnetic coupling phenomenon |
| FR2815124A1 (fr) | 2000-09-30 | 2002-04-12 | Schlumberger Services Petrol | Procede de determination de la saturation en hydrocarbure d'une formation |
| RU2253137C2 (ru) | 2003-07-02 | 2005-05-27 | Федеральное государственное унитарное предприятие Забайкальский комплексный научно-исследовательский институт (ФГУП ЗабНИИ) | Способ анализа геоэлектроразведочного сигнала |
| RU2231089C1 (ru) | 2003-07-08 | 2004-06-20 | РЫХЛИНСКИЙ Николай Иванович | Способ геоэлектроразведки (варианты) |
| RU2235347C1 (ru) | 2003-08-11 | 2004-08-27 | РЫХЛИНСКИЙ Николай Иванович | Способ геоэлектроразведки (варианты) |
| RU2301431C2 (ru) * | 2005-03-24 | 2007-06-20 | Общество с ограниченной ответственностью "Сибирская геофизическая научно-производственная компания" | Способ электроразведки с использованием пространственного дифференцирования поля становления на нескольких разносах |
| RU2284555C1 (ru) | 2005-06-01 | 2006-09-27 | Николай Иванович РЫХЛИНСКИЙ | Способ морской геоэлектроразведки с фокусировкой электрического тока (варианты) |
| US8447521B2 (en) * | 2009-03-13 | 2013-05-21 | Schlumberger Technology Corporation | Subsurface electrical conductivity contrast reconstruction using narrow-band electromagnetic data and fullwave imaging |
-
2008
- 2008-10-23 RU RU2008141889/28A patent/RU2399931C2/ru active
- 2008-11-17 CN CN200880132434.XA patent/CN102265188B/zh not_active Expired - Fee Related
- 2008-11-17 WO PCT/RU2008/000711 patent/WO2010047613A1/ru active Application Filing
- 2008-11-17 MX MX2011002163A patent/MX2011002163A/es active IP Right Grant
- 2008-11-17 EP EP08877586.1A patent/EP2352041A4/en not_active Withdrawn
-
2009
- 2009-10-19 US US12/588,507 patent/US8138763B2/en not_active Expired - Fee Related
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2491580C1 (ru) * | 2012-02-07 | 2013-08-27 | Общество с ограниченной ответственностью "Сибирская геофизическая научно-производственная компания" | Способ измерения геофизических характеристик с применением последующей инверсии геоэлектрических данных с дополнительным временным фильтром |
| RU2494419C1 (ru) * | 2012-03-11 | 2013-09-27 | Открытое акционерное общество Научно-производственное предприятие "Научно-исследовательский и проектно-конструкторский институт геофизических исследований геологоразведочных скважин" (ОАО НПП "ВНИИГИС") | Способ геоэлектроразведки |
| RU2595327C1 (ru) * | 2015-03-03 | 2016-08-27 | Общество с ограниченной ответственностью "ЛУКОЙЛ-Приморьенефтегаз" | Устройство сейсморазведки 2d или 3d, электроразведки и гис для повышения надежности картирования кровли соли и для прогноза нефтегазоносности подсолевых отложений в районах с развитой солянокупольной тектоникой |
Also Published As
| Publication number | Publication date |
|---|---|
| EP2352041A1 (en) | 2011-08-03 |
| CN102265188A (zh) | 2011-11-30 |
| MX2011002163A (es) | 2011-10-06 |
| US20100164500A1 (en) | 2010-07-01 |
| RU2008141889A (ru) | 2010-04-27 |
| EP2352041A4 (en) | 2017-08-02 |
| US8138763B2 (en) | 2012-03-20 |
| WO2010047613A1 (ru) | 2010-04-29 |
| CN102265188B (zh) | 2014-07-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2399931C2 (ru) | Способ количественного разделения эффектов электромагнитной индукции и вызванной поляризации | |
| RU2301431C2 (ru) | Способ электроразведки с использованием пространственного дифференцирования поля становления на нескольких разносах | |
| Fiandaca et al. | Resolving spectral information from time domain induced polarization data through 2-D inversion | |
| Karaoulis et al. | Time-lapse three-dimensional inversion of complex conductivity data using an active time constrained (ATC) approach | |
| Strahser et al. | Dependence of seismoelectric amplitudes on water content | |
| Colombo et al. | Quantifying surface-to-reservoir electromagnetics for waterflood monitoring in a Saudi Arabian carbonate reservoir | |
| Zeng | Thickness imaging for high-resolution stratigraphic interpretation by linear combination and color blending of multiple-frequency panels | |
| US7737699B2 (en) | Method of marine electromagnetic survey using focusing electric current | |
| Qin et al. | Fast prediction method of Archie’s cementation exponent | |
| US20160084980A1 (en) | Method and apparatus for detecting and mapping subsurface anomalies | |
| Souza et al. | Estimation of reservoir fluid saturation from 4D seismic data: effects of noise on seismic amplitude and impedance attributes | |
| Newgord et al. | Improved interpretation of electrical resistivity measurements in mixed-wet rocks: An experimental core-scale application and model verification | |
| Kang et al. | The effect of power lines on time-domain airborne electromagnetic data | |
| Attwa et al. | Resistivity characterization of aquifer in coastal semiarid areas: an approach for hydrogeological evaluation | |
| Dalgaard et al. | Enhancing SNMR model resolution by selecting an optimum combination of pulse moments, stacking, and gating | |
| Römhild et al. | Imaging hydraulic conductivity in near-surface aquifers by complementing cross-borehole induced polarization with hydraulic experiments | |
| Strack et al. | CCUS plume monitoring: Verifying surface CSEM measurements to log scale | |
| Houzé et al. | Combining experimental and modelling approaches to monitor the transport of an artificial tracer through the hyporheic zone | |
| Comina et al. | Joint inversion of VES and surface wave data | |
| Newgord et al. | A New Workflow for Joint Interpretation of Electrical Resistivity and NMR Measurements to Simultaneously Estimate Wettability and Water Saturation | |
| Butler | Research Note: The mean sensitivity depth of the electrical resistivity method | |
| Bobachev et al. | Features of instability in solving the inverse problem of vertical electrical sounding for precision monitoring | |
| Ageenkov et al. | Manifestation of induction and induced polarization in the case of axial and symmetrical electrical arrays | |
| RU2235347C1 (ru) | Способ геоэлектроразведки (варианты) | |
| RU2491580C1 (ru) | Способ измерения геофизических характеристик с применением последующей инверсии геоэлектрических данных с дополнительным временным фильтром |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| QB4A | Licence on use of patent |
Free format text: LICENCE Effective date: 20130221 |