RU2595327C1 - Device for seismic survey 2d or 3d, electrical exploration and gis to increase reliability of mapping the roof of salt and for prediction of oil and gas content subsalt deposits in areas with developed salt-dome tectonics - Google Patents
Device for seismic survey 2d or 3d, electrical exploration and gis to increase reliability of mapping the roof of salt and for prediction of oil and gas content subsalt deposits in areas with developed salt-dome tectonics Download PDFInfo
- Publication number
- RU2595327C1 RU2595327C1 RU2015107361/28A RU2015107361A RU2595327C1 RU 2595327 C1 RU2595327 C1 RU 2595327C1 RU 2015107361/28 A RU2015107361/28 A RU 2015107361/28A RU 2015107361 A RU2015107361 A RU 2015107361A RU 2595327 C1 RU2595327 C1 RU 2595327C1
- Authority
- RU
- Russia
- Prior art keywords
- seismic
- exploration
- model
- logging
- data
- Prior art date
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
Изобретение относится к геологоразведке, в частности к техническим средствам проведения региональных и поисково-разведочных геологических работ при прямых поисках и разведке нефтегазовых месторождений.The invention relates to geological exploration, in particular to technical means of conducting regional and exploration geological work in direct searches and exploration of oil and gas fields.
Назначением изобретения является использование комплекса для проведения поисково-разведочных работ на нефть и газ в сложно построенных районах с развитой соляно-купольной тектоникой.The purpose of the invention is the use of a complex for prospecting and exploration for oil and gas in complex areas with developed salt-dome tectonics.
В частности, комплекс геофизической разведки обеспечивает проведение комплексных полевых работ по единой сети профилей для реализации комплексных измерений сейсмического и электромагнитного полей с использованием дифференциально-нормированного метода электроразведки (ДНМЭ), совместного анализа полученных данных в единой информационной среде. При этом разведка с применением данного Комплекса применяется с целью повышения надежности картирования кровли соляных структур и изучения их внутреннего строения, что существенно влияет на достоверность структурных построений подсолевых отложений и прогноз их нефтегазоносности.In particular, the complex of geophysical exploration provides comprehensive field work on a single network of profiles for the implementation of complex measurements of seismic and electromagnetic fields using the differential-normalized method of electrical exploration (DNME), a joint analysis of the data in a single information environment. At the same time, exploration using this Complex is used to increase the reliability of mapping the roof of salt structures and study their internal structure, which significantly affects the reliability of structural constructions of subsalt deposits and the forecast of their oil and gas potential.
Известно изобретение "Способ геоэлектроразведки", патент RU №2210092, опубл. 10.08.2003, МПК G01V 3/06, в котором возбуждение электрического поля осуществляют с помощью питающих электродов, регистрируют и измеряют на выходе коммутирующего устройства разности потенциалов между каждой парой соседних приемных гальваноемкостных электродов, применяется устройство, последовательно коммутирующее на две выходные клеммы соседних приемных гальваноемкостных электродов. Изобретение относится к электроразведке малых глубин и может быть использовано при изучении геоэлектрической неоднородности верхней части разреза при инженерно-геологических изысканиях в сложных условиях заземлении (мерзлый грунт, сухие пески, твердые искусственные покрытия). Однако в отличие от предложенного комплекса, электроды не имеют заземления. В данном способе определяют только кажущееся сопротивление, использующееся для проведения инженерно-геологических изысканий, в отличие от предложенного комплекса, в котором определяют не только кажущееся сопротивление, но и послойное распределение удельного электрического сопротивления и коэффициента поляризуемости после разделения индукционной и гальванической составляющей поля DU. Кроме того, возбуждают электрическое поле частотой от первых сотен Гц до первых кГц, в то время, как в предложенном комплексе диапазон частот больше - от 0 до первых кГц.The invention is known "Method of geoelectrical exploration", patent RU No. 22210092, publ. 08/10/2003, IPC
Известно изобретение "Способ определения свойств подземных формаций", Патент RU 2294547, опубл. 27.02.2007, МПК G01V 3/38, в котором оптимизируют положение пластов и оценивают физические/нефтефизические свойства формации на основе измеренных электрических свойств, способ позволяет объединить параметры неизвестных формаций действием под названием «группировка параметров» и позволяет проверить совместимость данных сопротивления с указанными допусками. Однако данным способом измеряют электрические свойства формации в буровой скважине и разрабатывают модель параметров формации, а положение распределенных пластов моделируют измерениями электрических свойств при помощи определенных положений пластов, принимая, что один или более из указанных электрических свойств равны на всем протяжении некоторых из указанных распределенных пластов. Т.е. способ не позволяет изучать латеральное распределение продольного сопротивления и параметров вызванной поляризации, по которым осуществляется оконтуривание залежи углеводородов на больших расстояниях от скважины.The invention is known "Method for determining the properties of underground formations", Patent RU 2294547, publ. 02/27/2007, IPC
Известно изобретение "Способ измерения геофизических характеристик с применением последующей инверсии геоэлектрических данных с дополнительным временным фильтром", патент RU 2491580, опубл. 27.08.2013, МПК G01V 3/08, в котором определяют послойное распределение удельного электрического сопротивления (проводимости) геологического разреза с выбором количества слоев и их толщин. Однако используют 1D-инверсию, которая не позволяет картировать сложные протяженные трехмерные объекты, характерные для районов со сложной соляно-купольной тектоникойThe invention is known "A method of measuring geophysical characteristics using the subsequent inversion of geoelectric data with an additional time filter", patent RU 2491580, publ. 08/27/2013, IPC
Известно изобретение "Способ электроразведки с использованием пространственного дифференцирования поля становления на нескольких разносах", патент RU 2301431, опубл. 20.06.2007, МПК G01V 3/06, в котором возбуждают электромагнитное поле в среде, посылая в нее низкочастотную периодическую последовательность знакопеременных прямоугольных импульсов тока при прохождении зондирующей установки вдоль профиля исследований, измеряют и регистрируют мгновенные значения первых и вторых разностей потенциалов переходных процессов электрического поля на нескольких приемных линиях, применяют зондирующую установку с гибкой геометрией, состоящую из двух заземленных питающих линий, расположенных по обе стороны от нескольких приемных линий. Однако для достижения единственности решения необходимо привлечение дополнительной информации - комбинирование зондирования становлением поля с дистанционными зондированиями, использование различных способов возбуждения и приема поля и разных форм возбуждающих токовых импульсов, что не позволяет за счет 1D-инверсии картировать сложные протяженные трехмерные объекты, характерные для районов со сложной соляно-купольной тектоникой.The invention is known "Method of electrical exploration using spatial differentiation of the formation field at several spacings", patent RU 2301431, publ. 06/20/2007, IPC
Наиболее близким техническим решением является изобретение "Способ количественного разделения эффектов электромагнитной индукции и вызванной поляризации" патент RU 2399931, опубл. 20.09.2010, МПК G01V 3/38, G01V 3/06, в котором функцию измеренных параметров формируют как комбинацию временных и пространственных производных поля становления. Способ позволяет осуществить количественное определение интенсивности полей электромагнитной индукции и вызванной поляризации постоянного тока с погрешностью не более 0,5%, обеспечить большой временной интервал разделения полей от 1 мс до нескольких секунд, в том числе и над разрезами с высокой проводимостью. Однако 1D-инверсия не позволяет картировать сложные протяженные трехмерные объекты, характерные для районов со сложной соляно-купольной тектоникой, что, в свою очередь, не позволяет достичь повышения надежности картирования кровли соли, структурных построений подсолевых горизонтов и прогноза их нефтегазоносности в районах с развитой соляно- купольной тектоникой.The closest technical solution is the invention "Method for the quantitative separation of the effects of electromagnetic induction and induced polarization" patent RU 2399931, publ. 09/20/2010, IPC G01V 3/38, G01V 3/06, in which the function of the measured parameters is formed as a combination of temporal and spatial derivatives of the formation field. The method allows for the quantitative determination of the intensity of the fields of electromagnetic induction and the induced polarization of the direct current with an error of not more than 0.5%, to provide a large time interval for the separation of the fields from 1 ms to several seconds, including over sections with high conductivity. However, the 1D inversion does not allow the mapping of complex extended three-dimensional objects characteristic of areas with complex salt-dome tectonics, which, in turn, does not allow to increase the reliability of the mapping of the roof of salt, structural constructions of subsalt horizons and the forecast of their oil and gas potential in areas with developed salt - dome tectonics.
При разработке предложенного комплекса требуется определить глубинное положение кровли соленосных отложений, выделенное по данным электро- и сейсморазведки, и преобразовать в масштаб временного сейсмического разреза на основе закона взаимосвязи между временами регистрации сейсмического и электромагнитного сигналов (параметр «а»), (см. Смилевец Н.П. «Комплексирование геофизических полей на основе их адекватного представления в едином координатном пространстве // Диссертация на соискание ученой степени доктора технических наук // г. Москва 1999 г. с. 255). Для этого необходимо параметр «а» рассчитать по данным сейсмического (СК) каротажа или вертикального сейсмического профилирования (ВСП) и точкам электромагнитных зондирований, полученным вблизи параметрических скважин, а затем по данным годографов для одних и тех же глубин с определенным шагом дискретизации снять значения и , по которым в дальнейшем построить графики зависимости tв(tэ) и рассчитать параметр «а» по формуле:When developing the proposed complex, it is required to determine the deep position of the roof of saline sediments, selected according to electrical and seismic data, and convert it to the scale of a temporary seismic section based on the law of the relationship between the times of registration of seismic and electromagnetic signals (parameter “a”), (see Smilevec N .P. “Integration of geophysical fields on the basis of their adequate representation in a single coordinate space // The dissertation for the degree of Doctor of Technical Sciences // Mr. Mo Squa 1999, p. 255). For this, parameter “a” must be calculated using seismic (SC) logging data or vertical seismic profiling (VSP) and electromagnetic sounding points obtained near parametric wells, and then using values of hodographs for the same depths with a certain sampling step to take values and by which to further plot the dependences of t in (t e ) and calculate the parameter "a" according to the formula:
где - время регистрации сейсмического сигнала, а - это время регистрации электромагнитного сигнала. Where - time of registration of the seismic signal, and - this is the time of registration of the electromagnetic signal.
Тогда для преобразования геоэлектрических характеристик в масштаб временного сейсмического разреза можно использовать формулу, устанавливающую взаимосвязь параметра "а" - параметра взаимосвязи между временами регистрации сейсмического и электромагнитного сигналов, параметра "V" - средней скорости продольной волны, и параметра "σ" - кажущуюся продольную проводимость (или сопротивление).Then, to convert the geoelectric characteristics to the time scale of the seismic section, you can use the formula that establishes the relationship between the parameter "a" - the relationship between the recording times of the seismic and electromagnetic signals, the parameter "V" - the average velocity of the longitudinal wave, and the parameter "σ" - the apparent longitudinal conductivity (or resistance).
При этих условиях необходимо создать комплекс, в котором возможно измерить параметр взаимосвязи между временами регистрации сейсмического и электромагнитного сигналов, который будет являться функцией средней скорости продольной волны и кажущейся продольной проводимости (или сопротивления) разреза. Возможность достижения такой взаимосвязи обусловлена тем, что при пересечении профилем геологических условий с резко меняющимися глубинами кровли соли в каждой зоне (мульда, купол, склоновая часть), изученной сейсмокаротажем (СК) или вертикальным сейсмическим профилированием (ВСП), рассчитывают каждый раз свой параметр связи «а», после чего осуществляют его интерполяцию и экстраполяцию по всему профилю (ссылку см. там же).Under these conditions, it is necessary to create a complex in which it is possible to measure the parameter of the relationship between the recording times of seismic and electromagnetic signals, which will be a function of the average longitudinal wave velocity and the apparent longitudinal conductivity (or resistance) of the section. The possibility of achieving this relationship is due to the fact that when the profile crosses geological conditions with sharply changing depths of the salt roof in each zone (trough, dome, slope) studied by seismic logging (SC) or vertical seismic profiling (VSP), each time they calculate their communication parameter “A”, after which it is interpolated and extrapolated over the entire profile (see link there).
Таким образом, на основе закона взаимосвязи сейсмических и электромагнитных времен, измеренных на базе комплекса параметров сейсмо- и электроразведки и ГИС, комплекс должен обеспечить расчет прогнозной модели интервальных скоростей. При этом необходимо получить уточненную конфигурацию соляных структур, их внутреннюю неоднородность, согласованную корреляцию подсолевых отложений и прогнозную скоростную модель, которые должны учитываться в дальнейшем при глубинных преобразованиях временных сейсмоэлектроразведочных разрезов. Прогнозная на основе Комплекса глубинно-скоростная модель может использоваться в качестве исходной в процедурах сейсмической инверсии, в расчетах скоростей суммирования путем спрямления годографа отраженной волны во временной и глубинной областях.Thus, based on the law of the relationship of seismic and electromagnetic times, measured on the basis of a set of parameters of seismic and electrical exploration and well logging, the complex should provide a calculation of the predicted model of interval speeds. In this case, it is necessary to obtain an updated configuration of salt structures, their internal heterogeneity, consistent correlation of subsalt deposits and a predictive velocity model, which should be taken into account in the future during deep transformations of temporary seismic-electrical exploration sections. The deep-speed model predicted on the basis of the Complex can be used as the initial model in seismic inversion procedures, in calculating the summation speeds by straightening the travel time curve of the reflected wave in the time and depth regions.
Посредством предложенного комплекса на временном сейсмическом разрезе требуется уточнять положение кровли соли по комплексным измеренным данным сейсмоэлектроразведки и изучить внутреннее строение соляных куполов с выделением в них рапоносных горизонтов, сложенных, преимущественно, терригенно-карбонатными отложениями, проявляющимися в электромагнитном поле как проводящие (низкоомные) горизонты, что создает благоприятные условия для их обнаружения в высокоомных соляных структурах.Using the proposed complex at a temporary seismic section, it is required to clarify the position of the salt roof using the complex measured seismic and electrical exploration data and study the internal structure of salt domes with the separation of rape horizons composed mainly of terrigenous-carbonate deposits, which manifest themselves in the electromagnetic field as conducting (low resistance) horizons, which creates favorable conditions for their detection in high-resistance salt structures.
Техническим результатом предложенного комплекса является расширение технических средств, в частности разработка схемы комплекса, за счет которого достигается повышение надежности картирования кровли соли и который позволяет осуществить структурные построения подсолевых горизонтов и прогноз их нефтегазоносности в районах с развитой соляно- купольной тектоникой.The technical result of the proposed complex is the expansion of technical means, in particular, the development of a complex scheme, due to which an increase in the reliability of mapping the roof of salt is achieved and which allows structural construction of subsalt horizons and a forecast of their oil and gas potential in areas with developed salt-dome tectonics.
Заявленный технический результат достигается за счет того, что Комплекс включаетThe claimed technical result is achieved due to the fact that the Complex includes
- систему геофизического исследования скважин (ГИС), скважины глубокого бурения (скважины глубокого бурения - параметрические скважины) для получения комплексной сейсмоэлектроразведочной информации в масштабе временного сейсмического разреза с помощью данных, полученных от сейсмокаротажа (СК) и электрокаротажа (ЭК), и единую сеть сейсморазведочных профилей для снятия показаний (измеренных данных) измерительными установками сейсмо- и электроразведки. Комплекс отличается тем, что он снабжен следующими измерительными установками с использованием систем наблюдений: установкой 2D или 3D сейсморазведки в широкоазимутальном исполнении с использованием блоковой, псевдодиагональной системы наблюдения, в которой все линии возбуждения (ЛВ) и линии приема (ЛП) расположены вкрест простирания соляных структур, образуя сеть сейсморазведочных профилей, и взаимно встречной соосной электроразведочной установкой AB-MN с использованием дифференциально-нормированного метода электроразведки (ДНМЭ), измерительные установки устанавливают вблизи параметрических скважин, изученных на основе данных ГИС с учетом априорной информации о развитии соляной тектоники на изучаемой территории, при этом в установке электроразведки измеренные данные рассчитывают вдоль сейсмических профилей 2D или вдоль электроразведочных профилей в рамках сети сейсморазведки 3D. Комплекс дополнительно снабжен телеметрической аппаратурой типа 428 XL и программно-технологическим комплексом (ПТК), связанным с сейсморазведочным (БСР) и электроразведочным (БЭР) блоками, которые связаны соответственно с сейсмокаротажем (СК) и электрпокаротажем (ЭК) ГИС, а посредством программно- технологического комплекса осуществляют комплексную интерпретацию полученных функций сейсмогеоэлектрических параметров в масштабе временного сейсмического разреза, и на основе интерпретации получают (строят) сейсмогеоэлектрическую модель, согласованную с данными ГИС. Например, в установке сейсморазведки используют 3D сейсмокаротаж (СК) или вертикальное сейсмическое профилирование (ВСП) или акустический каротаж (АК). Так в установке 3D сейсморазведки расстояние между линиями возбуждения и приема составляет 200 м; шаг наблюдений по профилям ЛП и ЛВ - 50 м, а номинальная кратность наблюдения составляет 324. В рамках сети 3D сейсморазведки выполняют многокомпонентную электроразведку с использованием модификации ДНМЭ, которая является чувствительной к наличию высокоомных объектов, а при (в рамках сети) 2D сейсморазведки, выполняют электроразведку с использованием модификации ДНМЭ вдоль сейсмических профилей 2D. В установке электроразведки ДНМЭ используют электрокаротаж: стандартный (КС) или боковой каротаж (БК) или боковое электрическое зондирование (БКЗ). Данные сейсмических и электрических профилей, измеренных установками 2D или 3D сейсморазведки и электроразведки с использованием модификации ДНМЭ, рассчитывают вблизи скважин, изученных ГИС посредством сейсмокаротажа (СК), акустического каротажа (АК) и электрокаротажа типа КС или БК или БКЗ. При этом система установок электроразведки с использованием модификации ДНМЭ может представлять собой систему ортогональных профилей, длина которых увеличена по сравнению с участком, где требуется определить глубины кровли соли, по меньшей мере, на 4 км в обе стороны. И, например, с помощью ПТК осуществляют расчет полученных от измерительной установки электроразведки, использующей модификацию ДНМЭ, измеренные величины (характеристики) и формирование набора параметров:- a system of geophysical exploration of wells (GIS), deep-hole wells (deep-hole wells — parametric wells) to obtain complex seismic and electrical exploration information on the scale of a temporary seismic section using data obtained from seismic logging (SC) and electric logging (EC), and a single network of seismic exploration profiles for taking readings (measured data) by measuring installations of seismic and electrical exploration. The complex is characterized in that it is equipped with the following measuring installations using observation systems: a 2D or 3D seismic survey installation in wide-azimuth design using a block, pseudo-diagonal observation system in which all excitation lines (LC) and reception lines (LP) are located across the strike of salt structures forming a network of seismic profiles, and a reciprocal coaxial AB-MN electrical prospecting installation using the differential-normalized method of electrical exploration (DNME), rev The test rigs are installed near parametric wells studied on the basis of well log data taking into account a priori information on the development of salt tectonics in the study area, while the measured data are calculated along the 2D seismic profiles or along the electric prospecting profiles within the 3D seismic network in the electrical exploration installation. The complex is additionally equipped with telemetry equipment of type 428 XL and a software and technology complex (PTK) associated with seismic (BSR) and electrical exploration (BER) blocks, which are associated with seismic logging (SC) and electric logging (EC) GIS, and through software and technological of the complex carry out a comprehensive interpretation of the obtained functions of seismic-geoelectric parameters on the scale of a temporary seismic section, and on the basis of the interpretation, a seismic-geoelectric model is obtained (built), voiced with GIS data. For example, 3D seismic logging (SC) or vertical seismic profiling (VSP) or acoustic logging (AK) is used in a seismic survey installation. So in a 3D seismic survey installation, the distance between the excitation and reception lines is 200 m; the observation step along the profiles of PL and LV is 50 m, and the nominal magnification of observation is 324. Within the 3D seismic network, multicomponent electrical exploration is performed using the DNME modification, which is sensitive to the presence of high-impedance objects, and when (within the network) 2D seismic exploration is performed electrical exploration using a modification of DNME along 2D seismic profiles. In the DNME electrical exploration installation, electric logging is used: standard (KS) or lateral logging (BK) or lateral electric sounding (BKZ). The data of seismic and electrical profiles measured by 2D or 3D seismic and electrical exploration installations using the DNME modification are calculated near wells studied by well logging using seismic logging (SC), acoustic logging (AK) and electric logging of the KS or BK or BKZ type. In this case, the system of electrical exploration installations using the DNME modification can be a system of orthogonal profiles, the length of which is increased compared to the site where it is necessary to determine the depth of the roof of the salt by at least 4 km in both directions. And, for example, using the PTC, they calculate the electrical exploration obtained from the measuring installation using the DNME modification, the measured values (characteristics) and the formation of a set of parameters:
- DU(t)=ΔU(t)/ΔU0, где ΔU0 - разность потенциалов ΔU, измеренная во время пропускания тока;- DU (t) = ΔU (t) / ΔU 0 , where ΔU 0 is the potential difference ΔU, measured during the passage of current;
- Ps(t)=Δ2U(t)/ΔU(t), который рассчитан в режиме включения тока;- Ps (t) = Δ 2 U (t) / ΔU (t), which is calculated in the current on mode;
- P1(t)=Δ2U(t)/ΔU(t), который рассчитан в режиме выключения тока;- P 1 (t) = Δ 2 U (t) / ΔU (t), which is calculated in the current shutdown mode;
и Dφ(t)=It(t)-P1(t), где It(t)=Δt(Δ2U)/Δt(ΔU), поступающих на вход программного блока инверсии 3D для разделения исходного сигнала ДНМЭ на электродинамическую составляющую, используемую для расчета сопротивления среды, и поляризационную составляющую, используемую для расчета параметров вызванной поляризации, и обе составляющие, в которых учтено искажающее влияние на сигнал ДНМЭ элементов соляной тектоники, используют для интерпретации в комплексе с данными сейсморазведки и ГИС в едином информационном пространстве. Таким образом, ПТК обеспечивает проведение обработки и комплексной интерпретации разнородной геолого-геофизической интерпретации.and Dφ (t) = I t (t) -P 1 (t), where I t (t) = Δt (Δ 2 U) / Δt (ΔU), which are input to the 3D inversion software block to divide the initial DNME signal into the electrodynamic component used to calculate the resistance of the medium, and the polarization component used to calculate the parameters of the induced polarization, and both components that take into account the distorting effect of salt tectonics on the DNME signal, are used for interpretation in conjunction with seismic data and well logging in a single information space . Thus, the PTK provides processing and a comprehensive interpretation of a heterogeneous geological and geophysical interpretation.
Например, с помощью ПТК осуществляют расчет полученных от измерительной установки электроразведки, использующей модификацию ДНМЭ, измеренные величины (характеристики) и формирование набора параметров:For example, using the PTC, they calculate the electrical exploration obtained from the measuring installation using the DNME modification, the measured values (characteristics) and the formation of a set of parameters:
- DU(t)=ΔU(t)/ΔU0, где ΔU0 - разность потенциалов ΔU, измеренная во время пропускания тока;- DU (t) = ΔU (t) / ΔU 0 , where ΔU 0 is the potential difference ΔU, measured during the passage of current;
- Ps(t)=Δ2U(t)/ΔU(t), который рассчитан в режиме включения тока;- Ps (t) = Δ 2 U (t) / ΔU (t), which is calculated in the current on mode;
- P1(t)=Δ2U(t)/ΔU(t), который рассчитан в режиме выключения тока;- P 1 (t) = Δ 2 U (t) / ΔU (t), which is calculated in the current shutdown mode;
- и DΔ(t)=It(t)-P1(t), где It(t)=Δt(Δ2U)/Δt(ΔU),- and DΔ (t) = I t (t) -P 1 (t), where I t (t) = Δt (Δ 2 U) / Δt (ΔU),
поступающих на вход программного блока инверсии 3D для разделения исходного сигнала ДНМЭ на электродинамическую составляющую, используемую для расчета сопротивления среды, и поляризационную составляющую, используемую для расчета параметров вызванной поляризации, и обе составляющие, в которых учтено искажающее влияние на сигнал ДНМЭ элементов соляной тектоники, используют для интерпретации в комплексе с данными сейсморазведки и ГИС в едином информационном пространстве. Так Программный технологический комплекс (ПТК) для реализации комплексной интерпретации данных сейсмоэлектроразведки и ГИС может состоять из: Блока сейсмической инверсии (БСИ), на вход которого подают совокупность характеристик нескольких глубинно-скоростных разрезов, получая на выходе набор глубинно - скоростных моделей, Преобразователя (конвертора) геоэлектрических параметров и данных ГИС в масштаб временного сейсмического разреза, осуществляющего также прогнозирование скоростных параметров среды по комплексу данных сейсмоэлектроразведки и ГИС, Блока формирования согласованной сейсмогеоэлектрической модели по кровле соляных структур, их внутреннего строения и подсолевых горизонтов с использованием данных ГИС и Блока построения уточненных рельефов кровли соли и подсолевых горизонтов и прогноза нефтегазоносности подсолевых горизонтов (Блок построения и прогноза нефтегазоносности) с использованием прогнозных скоростных моделей, полученных по комплексу данных сейсмоэлектроразведки и ГИС. При этом Блок сейсмической инверсии может включать набор модулей, например 3 модуля: первый модуль расчета сейсмической модели на основании глубинно-скоростного разреза; второй модуль сравнения построенной сейсмической модели с наблюденным волновым полем; третий модуль множителя вариаций глубинно-скоростных разрезов в пределах заданного диапазона параметров - глубина, скорость, плотность, на выходе его получают набор глубинно-скоростных моделей в заданном количестве, рассчитанных вероятностным способом в пределах заданного диапазона. Все три модуля связаны между собой и с БСР прямой и обратной связью, и каждую глубинно-скоростную модель подают на вход модуля расчета сейсмической модели, а полученную сейсмическую модель, рассчитанную по каждой глубинно-скоростной модели, передают на модуль сравнения построенной сейсмической модели БСИ вместе с наблюденным волновым полем по сейсмическим параметрам, что обеспечивает расчет оптимальной сейсмической модели с наименьшей погрешностью расхождения с параметрами наблюденного волнового поля, полученной на выходе. При этом последовательное взаимодействие блоков и модулей, например. Обеспечивает после осуществления требуемых взаимодействий блоков попадание оптимальной сейсмической модели на вход блока формирования согласованной сейсмогеоэлектрической временной модели, и полученную в каждой точке измерения модификацией ДНМЭ на выходе этого блока сейсмоэлектроразведочную временную модель сопоставляют в БЭР с наблюденным электромагнитным сигналом на основе 3D инверсии до получения максимального совпадения с наблюденным электромагнитным полем на выходе БЭР, и которая поступает на вход Преобразователя, на выходе которого получают конечный рассчитанный параметр «а», который поступает на вход Преобразователя, на выходе которого получают конечный рассчитанный параметр «а», который в дальнейшем используют при корректировке прогнозной скоростной модели, а скорректированная прогнозная скоростная модель поступает на вход Блока построения и прогноза нефтегазоносности. Причем Блок сейсмической инверсии, Преобразователь, Блок формирования согласованной сейсмогеоэлектрической модели и Блок построения, и прогноза нефтегазоносности связаны прямой и обратной связью между собой. На вход Преобразователя поступают данные (сигналы) от БСР, БЭР и данные ГИС, а на каждый вход Блок формирования согласованной сейсмогеоэлектрической модели и на вход Блока построения и прогноза нефтегазоносности поступают данные (сигналы) от БСР, БЭР и Преобразователя. Блок построения и прогноза нефтегазоносности имеет выход на графопостроитель или интерфейс.received at the input of the 3D inversion program block to separate the initial DNME signal into the electrodynamic component used to calculate the medium resistance, and the polarization component used to calculate the induced polarization parameters, and both components that take into account the distorting effect of salt tectonics on the DNM signal, for interpretation in conjunction with seismic and GIS data in a single information space. So, the Program Technological Complex (PTC) for implementing a comprehensive interpretation of seismic data and well logging can consist of: A seismic inversion unit (BSI), the input of which serves as a set of characteristics of several deep-speed sections, receiving at the output a set of deep-speed models, a Converter (converter ) of geoelectric parameters and GIS data to the scale of a temporary seismic section, which also predicts the velocity parameters of the medium using the seismoelectrics data set exploration and well logging, the formation of a consistent seismic geoelectric model for the roof of salt structures, their internal structure and subsalt horizons using data from a well survey and the construction block of refined reliefs of the roof of salt and subsalt horizons and the prediction of oil and gas potential of subsalt horizons (block for constructing and predicting oil and gas content) using predicted velocity models obtained from a set of seismic electrical exploration and well logging data. In this case, the Seismic inversion unit may include a set of modules, for example, 3 modules: the first module for calculating the seismic model based on the depth-velocity section; the second module comparing the constructed seismic model with the observed wave field; the third module of the multiplier of variations of deep-speed sections within a given range of parameters - depth, speed, density, at the output it receives a set of deep-speed models in a given quantity, calculated by a probabilistic method within a given range. All three modules are connected with each other and with direct and feedback BSS, and each deep-speed model is fed to the input of the seismic model calculation module, and the obtained seismic model, calculated for each deep-speed model, is transmitted to the comparing module of the constructed BSM seismic model together with the observed wave field by seismic parameters, which ensures the calculation of the optimal seismic model with the smallest error of discrepancy with the parameters of the observed wave field obtained at the output. In this case, the sequential interaction of blocks and modules, for example. After performing the required interactions of the blocks, it ensures that the optimal seismic model hits the input of the matched seismic-geoelectric time model, and the seismic-electrical exploration time model obtained at each measurement point by the DNME output at this block is compared with the observed electromagnetic signal based on 3D inversion to obtain the maximum match the observed electromagnetic field at the output of the BER, and which is fed to the input of the Converter, to the output which give the final calculated parameter "a" that is input to the converter, at whose output is obtained a final calculated parameter 'a', which is then used for adjusting the target velocity model, and the corrected predictive velocity model is input to Block construction and oil and gas prediction. Moreover, the block of seismic inversion, the converter, the block for the formation of a consistent seismic-geoelectric model and the block for constructing and forecasting oil and gas content are connected by direct and feedback. Data (signals) from the BSR, BER and GIS data are received at the input of the Converter, and data (signals) from the BSR, BER and the Converter are received at each input of the Coordinated seismic-geoelectric model formation unit and at the input of the Oil and gas potential construction and forecast block. The block for constructing and forecasting oil and gas content has access to a plotter or interface.
Предложенное техническое решение иллюстрируется чертежами, которые не отражают всех возможных вариантов реализации, но демонстрируют работу комплекса и результаты этой работы.The proposed technical solution is illustrated by drawings, which do not reflect all possible implementation options, but demonstrate the operation of the complex and the results of this work.
На Фиг. 1 показана схема расположения блоков комплекса.In FIG. 1 shows the arrangement of blocks of the complex.
На Фиг. 2 показана функциональная схема Комплекса.In FIG. 2 shows a functional diagram of the Complex.
На Фиг. 3 показана двухсторонняя дипольная криволинейная установка с использованием модификации ДНМЭ с двумя питающими линиями (линиями возбуждения) AB и приемной линией MN с точкой заземления О.In FIG. Figure 3 shows a two-sided dipole curvilinear installation using a DNME modification with two supply lines (excitation lines) AB and a receiving line MN with ground point O.
На Фиг. 4 показано внутреннее строение соляных структур, построенное на базе сейсмоэлектроразведочного временного разреза, полученного в пределах АГКМ в результате использования Комплекса. а) показано выделение кровли соли на схеме интервального сопротивления Rн; б) показан расчет коэффициента «а» - график и фрагмент профиля с пикетами ДНМЭ, участвующими в расчете коэффициента «а»; в) показана схема изменения интервального сопротивления в масштабе (x, t0) без фильтрации и г) показана схема изменения интервального сопротивления в масштабе (x, t0) с фильтрацией.In FIG. Figure 4 shows the internal structure of salt structures built on the basis of a seismic-electrical exploration time section obtained within the AGKM as a result of using the Complex. a) shows the allocation of the roof of salt on the scheme of interval resistance Rн; b) the calculation of the coefficient “a” is shown — a graph and a fragment of the profile with pickets of DNME participating in the calculation of the coefficient “a”; c) a diagram of the variation of the interval resistance in the scale (x, t0) without filtering is shown; d) a diagram of the variation of the interval resistance in the scale (x, t0) with filtration is shown.
На Фиг. 5 показана схема площадного распределения аномалий вызванной поляризации ВП (ДНМЭ), полученная на основе 3D электромагнитной инверсии в пределах АГКМ.In FIG. Figure 5 shows a diagram of the areal distribution of anomalies of induced polarization of the VP (DNME), obtained on the basis of 3D electromagnetic inversion within the AGCM.
Комплекс устроен следующим образом.The complex is arranged as follows.
Имеются система геофизических исследований скважин (ГИС) (1) единая сеть сейсморазведочных профилей (2). В сеть сейсморазведочных профилей входят скважины глубокого бурения (параметрические скважины) (3) для сейсмокаротажа (СК) (4) и для электрокаротажа (КС) (5). Скважины (3) могут быть как несколько, так и одна скважина, а также для каждого каротажа своя. Скважины глубокого бурения (параметрические скважины) (3) устанавливают (используют имеющиеся) и используют для получения комплексной сейсмоэлектроразведочной информации в масштабе временного сейсмического разреза с помощью данных, полученных от сейсмокаротажа (СК) (4) и электрокаротажа (ЭК) (5), которые хранятся в ГИС (1). Это могут быть как априорные данные, снятые заранее, так и данные, полученные в процессе описываемой геофизической разведки.There is a system of geophysical research of wells (GIS) (1) a unified network of seismic profiles (2). The network of seismic profiles includes deep-hole wells (parametric wells) (3) for seismic logging (SK) (4) and for electric logging (KS) (5). Wells (3) can be several, or one well, and also for each logging own. Deep-hole boreholes (parametric wells) (3) are installed (use available) and used to obtain complex seismic-electrical exploration information in the scale of a temporary seismic section using data obtained from seismic logging (SC) (4) and electric logging (EC) (5), which stored in GIS (1). This can be a priori data taken in advance, or data obtained during the described geophysical exploration.
В единой сейсморазведочной сети профилей (2) осуществляют отработку сейсмоэлектроразведки по единой сети профилей измерительными установками (6, 7): сейсмо- (4) и электро- (5) каротажей (разведки).In a single seismic exploration network of profiles (2), seismic exploration is carried out according to a single network of profiles by measuring devices (6, 7): seismic (4) and electric (5) logs (exploration).
При этом используют скважины глубокого бурения (3), в которых изучают характеристики указанными выше видами каротажа: сейсмокаротажем (СК) (8), или вертикальным сейсмическим профилированием (ВСП), или акустическим каротажем (АК); и электрокаротажем (ЭК) (5) всех видов: стандартным (КС), боковым каротажем (БК), боковым электрическим зондированием (БКЗ). В установках 2D или 3D сейсморазведки (6) в широкоазимутальном исполнении используют блоковую, псевдодиагональную систему наблюдения, в которой все линии возбуждения (ЛВ) (8) и линии приема (ЛП) (9) расположены в крест простирания соляных структур, образуя сеть профилей (2) сейсморазведки. Так, в измерительных установках 2D или 3D сейсморазведки (6) в широкоазимутальном исполнении используют блоковую, псевдодиагональную систему наблюдения, в которую входят приемные линии (ЛП) (9). Измерительные установки (6) в сейсморазведке 2D или 3D используют следующие параметры: расстояниями между линиями возбуждения (8) и приема (9) - 200 м; шаг наблюдений по профилям ЛП (9) и ЛВ (8) - 50 м, номинальная кратность наблюдения 324.In this case, deep-hole drilling (3) is used, in which the characteristics of the above types of logging are studied: seismic logging (SK) (8), or vertical seismic profiling (VSP), or acoustic logging (AK); and electric logging (EC) (5) of all kinds: standard (CS), side logging (BK), side electric sounding (BKZ). In 2D or 3D seismic exploration installations (6) in wide-azimuth design, a block, pseudo-diagonal observation system is used, in which all excitation lines (LC) (8) and reception lines (LP) (9) are located in the strike direction of salt structures, forming a network of profiles ( 2) seismic exploration. So, in 2D or 3D seismic surveying installations (6) in wide-azimuth design, a block, pseudo-diagonal observation system is used, which includes receiving lines (PL) (9). Measuring installations (6) in 2D or 3D seismic exploit the following parameters: the distance between the excitation lines (8) and reception (9) is 200 m; the observation step along the profiles of PL (9) and drug (8) is 50 m, the nominal magnification of observation is 324.
В качестве измерительной установки электроразведки (7) используют взаимно встречную соосную установку электроразведки AB-MN с использованием дифференциально-нормированного метода электроразведки (модификация ДНМЭ), например двухстороннюю дипольную криволинейную установку (см. Фиг. 2) с двумя питающими линиями (линиями возбуждения) AB (10) и приемной линией MN (11) с точкой заземления О. Предложенный Комплекс использует расчет для установки электроразведки (7) на основе параметров (измеренных данных), полученных вдоль сейсмических профилей 2D.As a measurement installation of electrical exploration (7), a mutually counter-coaxial coaxial electrical exploration installation AB-MN is used using the differential-normalized electrical exploration method (modification of DNME), for example, a two-sided dipole curvilinear installation (see Fig. 2) with two supply lines (excitation lines) AB (10) and the receiving line MN (11) with a grounding point O. The proposed complex uses calculation to install electrical exploration (7) based on parameters (measured data) obtained along 2D seismic profiles.
Все изученные характеристики скважин в процессе работы Комплекса подают в базу данных ГИС (1). Измерительные установки (6, 7) на основе априорных данных, позволяющих определить территорию с развитой соляной тектоникой, установлены вблизи параметрических скважин (3), изученных с помощью ГИС (1).All studied well characteristics during the operation of the Complex are submitted to the GIS database (1). Measurement facilities (6, 7) based on a priori data that allow to determine the territory with developed salt tectonics, were installed near parametric wells (3), studied using well logging data (1).
Комплекс дополнительно снабжен телеметрической аппаратурой (12) типа 428 XL и программным технологическим комплексом (ПТК) (13), включающим в себя следующие структурные блоки (См. Фиг. 1).The complex is additionally equipped with telemetry equipment (12) of type 428 XL and a software technological complex (PTC) (13), which includes the following structural units (See. Fig. 1).
Программно-технологический комплекс (ПТК) (13) связан с сейсморазведочным блоком (БСР) (14) и электроразведочным блоком (БЭР) (15), которые в свою очередь получают данные из ГИС (1) и связаны соответственно с сейсмокаротажем (СК) (4) и электрокаротажем (ЭК) (5), с которых передают данные в ГИС (1). В ПТК (13) наряду с необходимыми электронными блоками для функционирования ПТК включены: Преобразователь (16) (конвертор) в масштаб временного сейсмического разреза геоэлектрических параметров среды, использующий данные из БСР (17) и БЭР (18), и ГИС (1); Блок формирования согласованной сейсмогеоэлектрической модели по кровле соляных структур, их внутреннего строения и подсолевых горизонтов (19) с использованием данных ГИС (1), Блок сейсмической инверсии (БСИ) (20) и Блок построения уточненных рельефов кровли соли и подсолевых горизонтов и прогноза нефтегазоносности подсолевых горизонтов (Блок построения и прогноза) (21) Блок сейсмической инверсии (20), Преобразователь (16), Блок формирования согласованной сейсмогеоэлектрической модели (19) и Блок построения уточненных рельефов кровли соли и подсолевых горизонтов и прогнозирования нефтегазоносности подсолевых горизонтов (21) (показаны условно).The software and technological complex (PTK) (13) is connected with the seismic block (BSS) (14) and the electric prospecting block (BER) (15), which in turn receive data from the GIS (1) and are associated with seismic logging (SC) ( 4) and electric logging (EC) (5), from which data are transmitted to the GIS (1). In addition to the necessary electronic components for the functioning of the PTC, the PTC (13) includes: A transducer (16) (converter) to the time scale of a seismic section of the geoelectric parameters of the medium, using data from the BSS (17) and the BER (18), and the GIS (1); A unit for the formation of a consistent seismic geoelectric model for the roof of salt structures, their internal structure and subsalt horizons (19) using GIS data (1), a block of seismic inversion (BSI) (20) and a unit for constructing refined reliefs of the roof of salt and subsalt horizons and a prediction of oil and gas potential of subsalt horizons (Building and forecasting block) (21) Seismic inversion block (20), Transducer (16), Coordinated seismic-geoelectric model formation block (19) and Refined salt roofing and under construction relief block OLEV horizon and the prediction horizon petrogas-salt (21) (shown in phantom).
В БСР 2D или 3D (14) осуществляют стандартную обработку данных, специализированную обработку, направленную на выделение кровли соли, и на выходе получают суммарный разрез общей глубинной точки (ОГТ) ГИС (1) получает данные с СК, ВСП, АК сейсмокаротажей и с КС, БК, БКЗ электрокаротажей. В БЭР с использованием модификации ДНМЭ (15) осуществляют стандартную обработку данных, 3D инверсию и на выходе получают разрез интегральных сопротивлений с выделением высокоомных соляных структур и параметры вызванной поляризации (ВП).In BDS 2D or 3D (14), standard data processing is carried out, specialized processing aimed at extracting the salt roof, and at the output, a total section of the common deep point (CDP) of the GIS is obtained (1) receives data from the SK, VSP, AK seismic logs and from the CS , BK, BKZ electric logging. In the BER, using the DNME modification (15), standard data processing, 3D inversion are performed and, at the output, a section of integrated resistances is obtained with the allocation of high-resistance salt structures and parameters of induced polarization (VP).
Преобразователь (16) или конвектор осуществляет расчет параметра взаимосвязи «а» и прогноза скоростной модели на основе параметра «а».The converter (16) or convector calculates the relationship parameter “a” and the forecast of the velocity model based on the parameter “a”.
Блок формирования согласованной сейсмогеоэлектрической модели (19) формирует модель по кровле соляных структур.The coordinated seismic-geoelectric model formation unit (19) forms the model along the roof of salt structures.
Блок сейсмической инверсии (20) включает 3 модуля (см. Фиг. 2). Первый модуль расчета сейсмической модели (22) на основании глубинно-скоростного разреза; второй модуль сравнения построенной сейсмической модели с наблюденным волновым полем (23); третий модуль множителя вариаций глубинно-скоростных разрезов (24) в пределах заданного диапазона параметров - глубина, скорость, плотность, на выходе его получают набор глубинно-скоростных моделей в заданном количестве, рассчитанных вероятностным способом в пределах заданного диапазона. Все три модуля связаны между собой и с БСР прямой и обратной связью.The seismic inversion unit (20) includes 3 modules (see Fig. 2). The first module for calculating the seismic model (22) based on the depth-velocity section; the second module comparing the constructed seismic model with the observed wave field (23); the third module of the multiplier of variations of deep-speed sections (24) within a given range of parameters - depth, speed, density, at the output it receives a set of deep-speed models in a given number, calculated by a probabilistic method within a given range. All three modules are connected with each other and with BSS direct and feedback.
Блок построения и прогноза (21) осуществляет построение уточненных рельефов кровли соли и подсолевых горизонтов и прогнозирует нефтегазоносность подсолевых горизонтов.The construction and prediction unit (21) constructs refined reliefs of the salt roof and subsalt horizons and predicts the oil and gas potential of subsalt horizons.
Все вышеперечисленные в ПТК (13) операционные блоки связаны между собой посредством прямой и обратной связи, что позволяет проводить комплексный анализ сейсмогеоэлектрических характеристик в интерактивном режиме с целью уточнения, согласования положения кровли соли и подсолевых отложений на временном сейсмическом разрезе. Такой комплексный подход в значительной степени повышает достоверность структурных построений кровли соли и подсолевых отложений и прогноз их нефтегазоносности в районах с развитой соляной тектоникой.All of the operational blocks listed in the PTC (13) are interconnected via direct and feedback, which allows a comprehensive analysis of seismic-geoelectric characteristics in an interactive mode to clarify and coordinate the position of the salt roof and subsalt deposits in a temporary seismic section. Such an integrated approach significantly increases the reliability of the structural structures of the salt roof and subsalt deposits and the forecast of their oil and gas potential in areas with developed salt tectonics.
С помощью Комплекса выполняют следующие операции, реализующие комплексную интерпретацию данных сейсмоэлектроразведки и ГИС:With the help of the Complex, the following operations are performed that implement a comprehensive interpretation of seismic data and well logging data:
- корреляцию основных сейсмических границ с учетом данных ГИС и построение глубинных и временных сейсмических моделей;- correlation of the main seismic boundaries taking into account GIS data and the construction of deep and temporary seismic models;
- сейсмическую инверсию для расчета сейсмической модели на основании глубинно-скоростного разреза и сравнения построенной сейсмической модели с наблюденным волновым полем, на выходе которого получается набор глубинно-скоростных моделей в заданном количестве, рассчитанных вероятностным способом в пределах заданного диапазона Полученная оптимальная сейсмическая модель подается на вход блока построения согласованной сейсмогеоэлектрической временной модели. Полученный сейсмоэлектроразведочный временной разрез в каждой точке электроразведки сопоставляется с наблюденным электромагнитным сигналом и в результате ряда итераций с целью максимального совпадения с наблюденным электромагнитным полем осуществляется конечная настройка параметра «а».- seismic inversion for calculating a seismic model based on a depth-velocity section and comparing the constructed seismic model with the observed wave field, the output of which is a set of depth-velocity models in a given number, calculated by a probabilistic method within a given range The resulting optimal seismic model is input block building a consistent seismic geoelectric time model. The obtained seismic-electrical exploration time section at each point of electrical exploration is compared with the observed electromagnetic signal and, as a result of a series of iterations, in order to maximize coincidence with the observed electromagnetic field, the final setting of parameter “a” is carried out.
- В интерактивном режиме - специализированную обработку электродинамической составляющей сигнала, полученного с помощью модификации ДНМЭ на основе согласования с данными электрокаротажа и с учетом геометрии глубинной сейсмической модели. Данная обработка направлена на выделение высокоомных соляных структур в интегральных характеристиках электромагнитного поля, с формированием глубинной геоэлектрической модели с рельефом кровли соленосных отложений;- In interactive mode - specialized processing of the electrodynamic component of the signal obtained by modifying the DNME based on coordination with the electric logging data and taking into account the geometry of the deep seismic model. This treatment is aimed at highlighting high-resistance salt structures in the integrated characteristics of the electromagnetic field, with the formation of a deep geoelectric model with a relief of the roof of saline deposits;
- прогноз скоростной модели по комплексу данных сейсмоэлектроразведки и ГИС (СК или ВСП, АК) с использованием закона взаимосвязи между временами регистрации сейсмического или электромагнитного сигналов (Преобразователь);- forecasting a speed model based on a set of data from seismic electrical exploration and well logging (SC or VSP, AK) using the law of the relationship between the times of registration of seismic or electromagnetic signals (Converter);
- формирование согласованной сейсмогеоэлектрической модели на временном сейсмическом разрезе с уточнением положения кровли соляных структур и подсолевых отложений.- formation of a coordinated seismic-geoelectric model on a temporary seismic section with clarification of the roof position of salt structures and subsalt deposits.
Таким образом, предложенный Комплекс обеспечивает комплексную интерпретацию полученных параметров (характеристик) геологической среды в плоскости (в характеристиках) временного сейсмического разреза, а также построение сейсмогеоэлектрических моделей, согласованных с данными ГИС, и прогноз нефтегазоносности подсолевых отложений.Thus, the proposed complex provides a comprehensive interpretation of the obtained parameters (characteristics) of the geological environment in the plane (in the characteristics) of the temporary seismic section, as well as the construction of seismic-geoelectric models consistent with GIS data and the prediction of oil and gas content of subsalt deposits.
Комплекс работает следующим образом.The complex works as follows.
Каждая из набора глубинно-скоростных моделей подается на вход в модуья расчета сейсмической модели (22), а полученная сейсмическая модель по каждой из реализаций глубинно-скоростной модели передается далее в модуль сравнения построенной сейсмической модели с наблюденным волновым полем (23) и таким образом рассчитывается сейсмическая модель с наименьшей погрешностью расхождения с параметрами наблюденного волнового поля. Полученная оптимальная сейсмическая модель подается на вход блока формирования согласованной сейсмогеоэлектрической временной модели (19). Полученный сейсмоэлектроразведочный временной разрез в каждой точке электроразведки сопоставляется с наблюденным электромагнитным сигналом и в результате ряда итераций с целью максимального совпадения с наблюденным электромагнитным полем осуществляется конечная настройка параметра «а» в Преобразователе (20), который в дальнейшем используется при корректировке прогнозной скоростной модели, а скорректированная прогнозная скоростная модель поступает на вход Блока построения и прогноза нефтегазоносности.Each of the set of deep-speed models is fed to the input to the calculation modules of the seismic model (22), and the obtained seismic model for each of the implementations of the deep-speed model is transferred further to the module for comparing the constructed seismic model with the observed wave field (23) and thus is calculated seismic model with the smallest error of discrepancy with the parameters of the observed wave field. The obtained optimal seismic model is fed to the input of the unit for the formation of a consistent seismic-geoelectric time model (19). The obtained seismic-electrical exploration time section at each point of the electrical exploration is compared with the observed electromagnetic signal and, as a result of a series of iterations, in order to maximize coincidence with the observed electromagnetic field, the final setting of parameter “a” in the Converter (20) is carried out, which is further used to adjust the predicted velocity model, and The adjusted forecast velocity model is fed to the input of the Block for constructing and forecasting oil and gas potential.
Для решения задач выделения кровли соли и подсолевых горизонтов в комплексе с сейсморазведкой используют многокомпонентную установку электроразведки (7) с применением модификации ДНМЭ, результаты которой отрабатывают вдоль линий сейсмопрофилей 2D или в рамках съемки сейсморазведки 3D в крест простирания соляных структур. Сейсмоэлектроразведочные профили (2) отрабатывают вблизи скважин (3), изученных ГИС (1), например, сейсмокаротажем (6) - СК, акустическим каротажем - АК и электрокаротажем (7), например, КС, БК, БКЗ.To solve the problems of separating the roof of salt and subsalt horizons in combination with seismic exploration, a multicomponent electrical exploration installation is used (7) using the DNME modification, the results of which are worked out along the lines of 2D seismic profiles or as part of 3D seismic surveying in the cross of the strike of salt structures. Seismic-electrical exploration profiles (2) are practiced near wells (3) studied by well logging (1), for example, seismic logging (6) - SK, acoustic logging - AK and electric logging (7), for example, KS, BK, BKZ.
В качестве измерительной установки электроразведки (7) берут двухстороннюю дипольную криволинейную установку с использованием модификации ДНМЭ (см. Фиг. 2). Двусторонняя дипольная взаимовстречная установка типа A1B1MONA2B2 с произвольной геометрией относительно оси (источника) ЛВ (10) со следующими параметрами:As a measuring installation of electrical exploration (7), a two-sided dipole curvilinear installation is used using the DNME modification (see Fig. 2). Two-way dipole reciprocal installation of type A 1 B 1 MONA 2 B 2 with arbitrary geometry relative to the axis (source) of the drug (10) with the following parameters:
Первая и вторая разности потенциалов поля становления и первичного поля регистрируются одновременно разными каналами, с использованием приемной аппаратуры со следующими характеристиками:The first and second potential differences of the formation field and the primary field are recorded simultaneously by different channels, using receiving equipment with the following characteristics:
Формирование второй разности потенциалов осуществляется аналоговой схемой перед поступлением на вход ПТК (13), что позволяет подавить синфазные помехи в диапазоне до 160 дБ. При этом используется телеметрическая аппаратура (12) и типы измерителей разности потенциалов геофизических ИРПГ, зарегистрированных в Федеральном агентстве по техническому регулированию и метрологии под номером 42970-09 RU.C.34.002.A или аналогичных. Аппаратура должна быть сертифицирована и своевременно пройти все поверки, указанные в паспорте на сертификацию.The second potential difference is formed by an analog circuit before entering the PTC input (13), which makes it possible to suppress common-mode interference in the range up to 160 dB. At the same time, telemetry equipment (12) and types of meters of potential difference of geophysical IRPGs registered with the Federal Agency for Technical Regulation and Metrology under the number 42970-09 RU.C.34.002.A or similar are used. The equipment must be certified and timely pass all the verification indicated in the certification passport.
Система наблюдений должна представлять собой систему ортогональных профилей (2), при этом длина профилей должна быть увеличена по сравнению с участком, где требуется определение глубины кровли соли, как минимум, на 4 км в обе стороны.The observation system should be a system of orthogonal profiles (2), while the length of the profiles should be increased in comparison with the site where it is necessary to determine the depth of the roof of the salt by at least 4 km in both directions.
В результате обработки сигналов ДНМЭ с помощью Комплекса формируется следующий набор параметров, поступающий на вход ПТК (13) для инверсии 3D:As a result of processing the DNME signals with the help of the Complex, the following set of parameters is generated, which is input to the PTC (13) for 3D inversion:
- DU(t)=ΔU(t)/ΔU0, где ΔU0 - разность потенциалов ΔU, измеренная во время пропускания тока;- DU (t) = ΔU (t) / ΔU 0 , where ΔU 0 is the potential difference ΔU, measured during the passage of current;
- Ps(t)=Δ2U(t)/ΔU(t), который рассчитан в режиме включения тока;- Ps (t) = Δ 2 U (t) / ΔU (t), which is calculated in the current on mode;
- P1(t)=Δ2U(t)/ΔU(t), который рассчитан в режиме выключения тока;- P 1 (t) = Δ 2 U (t) / ΔU (t), which is calculated in the current shutdown mode;
- Dϕ(t)=It(t)-P1(t), где It(t)=Δt(Δ2U)/Δt(ΔU).- Dϕ (t) = I t (t) -P 1 (t), where I t (t) = Δt (Δ 2 U) / Δt (ΔU).
Погрешность расчета данных параметров (как в режиме пропускания тока, так и при регистрации переходных процессов) не должна превышать следующих величин: 0.5% для DU(t), 0.008 - для Ps(t), P1(t) и Dϕ(t).The error in calculating these parameters (both in current transmission mode and during transient recording) should not exceed the following values: 0.5% for DU (t), 0.008 for Ps (t), P 1 (t) and Dϕ (t) .
Затем на основе использования данного набора параметров проводится количественное разделение полей электромагнитной индукции (EM) и вызванной поляризации (IP) во всем временном диапазоне [3], 3D-моделирование и нелинейная 3D-инверсия электромагнитных зондирований [4] по технологии ДНМЭ.Then, based on the use of this set of parameters, a quantitative separation of the fields of electromagnetic induction (EM) and induced polarization (IP) in the entire time range [3], 3D modeling and nonlinear 3D inversion of electromagnetic sounding [4] using DNME technology is carried out.
В Программно-технологическом комплексе (ПТК) (13) для реализации комплексной интерпретации данных сейсмоэлектроразведки (6, 7) и ГИС (1) посредством Блока сейсмической инверсии (2), имеющего в качестве исходного материала глубинно-скоростной разрез и включающего в себя, например, 3 модуля: модуль расчета сейсмической модели на основании глубинно-скоростного разреза (22); модуль сравнения построенной сейсмической модели с наблюденным волновым полем (23) и модуль множителя вариаций глубинно-скоростных разрезов (24) в пределах заданного диапазона параметров (глубина, скорость, плотность), получают на выходе модуля (24) набор глубинно-скоростных моделей в заданном количестве, рассчитанных вероятностным способом в пределах заданного диапазона. Все три модуля (22, 23, 24) связаны между собой и с БСР (6) прямой и обратной связью.In the Program and Technology Complex (PTC) (13) for the implementation of a comprehensive interpretation of seismic exploration data (6, 7) and GIS (1) by means of the Seismic Inversion Unit (2), which has a depth-velocity section as source material and includes, for example , 3 modules: a module for calculating a seismic model based on a depth-velocity section (22); a module for comparing the constructed seismic model with the observed wave field (23) and a module for the multiplier of variations of depth-velocity sections (24) within a given range of parameters (depth, speed, density), receive at the output of module (24) a set of depth-velocity models in a given the number calculated in a probabilistic manner within a given range. All three modules (22, 23, 24) are connected with each other and with BSS (6) direct and feedback.
Каждая из полученного набора глубинно-скоростных моделей подается на вход в модуль расчета сейсмической модели (22), а полученная сейсмическая модель по каждой из реализаций глубинно-скоростной модели передается далее в модуль сравнения построенной сейсмической модели (23) с наблюденным волновым полем и таким образом рассчитывается сейсмическая модель с наименьшей погрешностью расхождения с параметрами наблюденного волнового поля.Each of the obtained set of deep-speed models is fed to the input to the module for calculating the seismic model (22), and the obtained seismic model for each of the implementations of the deep-speed model is transmitted further to the module for comparing the constructed seismic model (23) with the observed wave field and thus the seismic model is calculated with the smallest error of discrepancy with the parameters of the observed wave field.
Полученная оптимальная сейсмическая модель подается на вход блока формирования (построения) согласованной сейсмогеоэлектрической временной модели (19). Полученный сейсмоэлектроразведочный временной разрез в каждой точке ДНМЭ сопоставляется с наблюденным электромагнитным сигналом и в результате ряда итераций с целью максимального совпадения с наблюденным электромагнитным полем осуществляется конечная настройка параметра «а» (20) Преобразователя (конвертора) (20) геоэлектрических параметров и данных ГИС в масштаб временного сейсмического разреза, осуществляющего также прогнозирование скоростных параметров среды по комплексу данных сейсмоэлектроразведки и ГИС (1), блока формирования согласованной сейсмогеоэлектрической модели (19) по кровле соляных структур, их внутреннего строения и подсолевых горизонтов с использованием данных ГИС (1), и блока построения уточненных рельефов кровли соли и подсолевых горизонтов и прогноза нефтегазоносности (21) подсолевых горизонтов с использованием прогнозных скоростных моделей, полученных по комплексу данных сейсмоэлектроразведки и ГИС (1).The resulting optimal seismic model is fed to the input of the formation (construction) block of the coordinated seismic-geoelectric time model (19). The obtained seismic-electrical exploration time section at each point of the DNME is compared with the observed electromagnetic signal and, as a result of a series of iterations, in order to maximize coincidence with the observed electromagnetic field, the parameter “a” (20) of the Converter (converter) (20) of the geoelectric parameters and GIS data is scaled of a temporary seismic section, which also predicts the velocity parameters of the medium using a set of seismic electrical exploration and well logging data (1), a formation unit the agreed seismic geoelectric model (19) for the roof of salt structures, their internal structure and subsalt horizons using GIS data (1), and the block for constructing refined reliefs of the roof of salt and subsalt horizons and the forecast of oil and gas potential (21) of subsalt horizons using predictive velocity models obtained according to a set of seismic electrical exploration and well logging data (1).
Пример работы Комплекса для представления результатов работы с применением измерительной установки электроразведки с использованием модификации ДНМЭ на Пойменном лицензионном участке в 2013 г. приведен ниже.An example of the Complex’s work for presenting the results of work with the use of a measuring installation of electrical exploration using a modification of DNME at the Poymenny license area in 2013 is given below.
Исследования проводились в Астраханской области на Пойменном лицензионном участке. Работы Комплекса осуществлялись с использованием ДНМЭ в условиях суши. Использовалась двусторонняя дипольная криволинейная установка с произвольной геометрией относительно оси источника. Использовалось гальваническое заземление.Research was conducted in the Astrakhan region on the Floodplain licensed site. The work of the Complex was carried out using DNME on land. A two-sided dipole curvilinear setup with arbitrary geometry relative to the axis of the source was used. Used galvanic grounding.
Разности потенциалов первого и второго порядка - AU и A2U замерялись непосредственно в ходе работ.The potential differences of the first and second order - AU and A 2 U were measured directly during the work.
С помощью ПТК (13) формируются следующие дифференциально-нормированные параметры (ДНП): DU=ΔU/ΔUпр,, D2U=Δ2U/ΔUпр, P1=Δ2U/ΔU, PS=(Δ2U-Δ2Uпр)/(ΔU-ΔUпр), P1пр=Δ2Uпр/ΔUпр. Посредством численного дифференцирования рассчитывают временные производные d(ΔU)/dt и d(Δ2U)/dt, из их отношения получают параметр Im=, который необходим для формирования параметра: Dφ=Im-P1. Проводят обработку данных по спадам, и после устранения тренда выполняют робастную 2D фильтрацию. Сглаженные кривые используют для расчета всего набора кривых ДНП и их доверительных интервалов: P1, P1пр, PS, Dφ или DφS. Кривые ДНП заносят в базу данных, где вычисляют осредненные параметры для взаимовстречных установок - блок (16). В конечном итоге этот набор параметров используют в качестве входных данных при дальнейшем решении задачи инверсии и построении геоэлектрического разреза в блоках (19, 20).With the PTC (13) are formed following differential normalized parameters (DNP): DU = ΔU / ΔU pr ,, D2U = Δ 2 U / ΔU etc., P1 = Δ 2 U / ΔU, PS = (Δ 2 2 U-Δ U pr ) / (ΔU-ΔU pr ), P1 pr = Δ 2 U pr / ΔU pr Using numerical differentiation, the time derivatives d (ΔU) / dt and d (Δ 2 U) / dt are calculated, from their ratio the parameter Im = is obtained, which is necessary for the formation of the parameter: Dφ = Im-P1. Processing data on recessions is carried out, and after eliminating the trend, robust 2D filtering is performed. Smooth curves are used to calculate the entire set of DNP curves and their confidence intervals: P1, P1 pr , PS, Dφ or DφS. The DNP curves are entered into the database, where the averaged parameters for mutual settings are calculated - block (16). Ultimately, this set of parameters is used as input for further solving the inversion problem and constructing a geoelectric section in blocks (19, 20).
Поскольку электромагнитное поле переходного процесса характеризуется пространственно-временной структурой и данный набор параметров является эффективным инструментом для ее изучения, который представлен в виде предложенного Комплекса, использование ДНП позволяет уверенно разделить - начиная с некоторого времени - влияние вызванной поляризации и электродинамических процессов. В результате из суммарного поля DU получаем поля EM и IPg [3]. Обязательной процедурой при проведении инверсии данных является проверка воспроизводимости полученных результатов.Since the electromagnetic field of the transition process is characterized by a spatio-temporal structure and this set of parameters is an effective tool for studying it, which is presented in the form of the proposed Complex, the use of DNP allows us to reliably separate, starting from some time, the effects of polarization and electrodynamic processes. As a result, from the total field DU we obtain the fields EM and IPg [3]. A mandatory procedure for conducting data inversion is to verify the reproducibility of the results.
На территории Пойменного участка для учета искажающего влияния соляно-купольной тектоники проведена 3D-инверсия, для осуществления которой использовалась программа GEMIR v0.0 (Новосибирский государственный технический университет), реализующая конечно-элементное 3D-моделирование и нелинейную 3D-инверсию электромагнитных зондирований [4] по технологии ДНМЭ. 3D-инверсия данных ДНМЭ выполнялась на основе сигналов DU и D2U.On the territory of the Floodplain site, 3D inversion was carried out to take into account the distorting effect of salt-dome tectonics, which was implemented using the GEMIR v0.0 program (Novosibirsk State Technical University), which implements finite element 3D modeling and non-linear 3D inversion of electromagnetic soundings [4] according to DNME technology. 3D inversion of the DNME data was performed based on the DU and D2U signals.
На карте (см. Фиг. 4) показаны этапы получения окончательной схемы изменения интервального сопротивления в масштабе (x, t0) с фильтрацией и без фильтрации и выделение кровли соли на схеме интервального сопротивления (Rн). При этом показан график расчета коэффициента «а».The map (see Fig. 4) shows the steps for obtaining the final scheme for changing the interval resistance in the scale (x, t0) with and without filtration and the allocation of the salt roof on the interval resistance scheme (Rн). In this case, a graph for calculating the coefficient "a" is shown.
По результатам 3D-инверсии данных ДНМЭ получена карта распределения аномалий ВП, согласующаяся с картой распределения аномалий, выделенных на основе инверсии 1-D, и карта положения кровли высокоомного (соленосного) горизонта (см. Фиг. 5).Based on the 3D inversion of the DNME data, a map of the distribution of VP anomalies is obtained, which is consistent with the map of the distribution of anomalies identified on the basis of 1-D inversion, and a map of the position of the roof of the high-resistance (saline) horizon (see Fig. 5).
Таким образом, за счет использования представленного комплекса сейсморазведки 2D или 3D, электроразведки и ГИС, направленного на повышение надежности картирования кровли соли и прогноза нефтегазоносности подсолевых отложений в районах с развитой соляно-купольной тектоникой, достигается заявленный технический результат.Thus, through the use of the presented complex of 2D or 3D seismic exploration, electrical exploration and well logging, aimed at improving the reliability of the mapping of the salt roof and prediction of oil and gas content of subsalt deposits in areas with developed salt-dome tectonics, the claimed technical result is achieved.
1. Смилевец Н.П. «Комплексирование геофизических полей на основе их адекватного представления в едином координатном пространстве // Диссертация на соискание ученой степени доктора технических наук // г. Москва 1999 г., с. 255.1. Smilevets N.P. “Integration of geophysical fields on the basis of their adequate representation in a single coordinate space // The dissertation for the degree of Doctor of Technical Sciences // Moscow 1999, p. 255.
2. Патент №2155977 (РФ) Способ геофизической разведки / Смилевец Н.П., Соколова И.П. // приор. От 19.05.1999 г., зарегистрирован в Государственном реестре изобретений Российской Федерации г. Москва 10 сентября 2000 г.2. Patent No. 2155977 (RF) Method of geophysical exploration / Smilevets NP, Sokolova IP // prior. From 05/19/1999, registered in the State Register of Inventions of the Russian Federation, Moscow, September 10, 2000
3. П.Ю. Легейдо, С.А. Иванов, Е.В. Агеенков и др. Способ количественного разделения эффектов электромагнитной индукции и вызванной поляризации. Патент RU 2399931 - Бюл. №12. - 2010.3. P.Yu. Legeido, S.A. Ivanov, E.V. Ageenkov et al. Method for the quantitative separation of the effects of electromagnetic induction and induced polarization. Patent RU 2399931 - Bull. No. 12. - 2010.
4. М.Г. Персова, Ю.Г. Соловейчик, Г.М. Тригубович, М.Г. Токарева. Методы и алгоритмы восстановления трехмерной структуры проводимости и поляризуемости среды по данным электромагнитных зондирований на основе конечно-элементного 3D-моделирования // Физика Земли, - 2013. - №3, - С. 30-45.4. M.G. Persova, Yu.G. Soloveichik, G.M. Trigubovich, M.G. Tokarev. Methods and algorithms for reconstructing the three-dimensional structure of the conductivity and polarizability of the medium according to electromagnetic sounding data based on finite element 3D modeling // Physics of the Earth, - 2013. - No. 3, - P. 30-45.
Claims (9)
- DU(t)=ΔU(t)/ΔU0, где ΔU0 - разность потенциалов ΔU, измеренная во время пропускания тока;
- Ps(t)=Δ2U(t)/ΔU(t), который рассчитан в режиме включения тока;
- P1(t)=Δ2U(t)/ΔU(t), который рассчитан в режиме выключения тока;
и Dφ(P(t)=It(t)-P1(t), где It(t)=Δt(Δ2U)/Δt(ΔU),
и которые поступают на вход программного блока инверсии 3D для разделения исходного сигнала ДНМЭ на электродинамическую составляющую, используемую для расчета сопротивления среды, и поляризационную составляющую, используемую для расчета параметров вызванной поляризации, и обе составляющие, в которых учтено искажающее влияние на сигнал ДНМЭ элементов соляной тектоники, используют для интерпретации в комплексе с данными сейсморазведки и ГИС в едином информационном пространстве.8. The device according to p. 1, characterized in that using the block of the software and technology complex (PTC), they calculate the received data from the measuring installation of electrical exploration using the DNME modification, which are measured and formed with the following set of parameters:
- DU (t) = ΔU (t) / ΔU 0 , where ΔU 0 is the potential difference ΔU, measured during the passage of current;
- Ps (t) = Δ 2 U (t) / ΔU (t), which is calculated in the current on mode;
- P 1 (t) = Δ 2 U (t) / ΔU (t), which is calculated in the current shutdown mode;
and Dφ ( P (t) = I t (t) -P 1 (t), where I t (t) = Δt (Δ 2 U) / Δt (ΔU),
and which are fed to the input of the 3D inversion program block for dividing the initial DNME signal into the electrodynamic component used to calculate the medium resistance, and the polarization component used to calculate the induced polarization parameters, and both components that take into account the distorting effect of salt tectonics on the DNME signal , used for interpretation in combination with seismic data and GIS in a single information space.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2015107361/28A RU2595327C1 (en) | 2015-03-03 | 2015-03-03 | Device for seismic survey 2d or 3d, electrical exploration and gis to increase reliability of mapping the roof of salt and for prediction of oil and gas content subsalt deposits in areas with developed salt-dome tectonics |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2015107361/28A RU2595327C1 (en) | 2015-03-03 | 2015-03-03 | Device for seismic survey 2d or 3d, electrical exploration and gis to increase reliability of mapping the roof of salt and for prediction of oil and gas content subsalt deposits in areas with developed salt-dome tectonics |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2014136106 Substitution | 2014-09-04 |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2595327C1 true RU2595327C1 (en) | 2016-08-27 |
Family
ID=56891916
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2015107361/28A RU2595327C1 (en) | 2015-03-03 | 2015-03-03 | Device for seismic survey 2d or 3d, electrical exploration and gis to increase reliability of mapping the roof of salt and for prediction of oil and gas content subsalt deposits in areas with developed salt-dome tectonics |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2595327C1 (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU172989U1 (en) * | 2016-12-28 | 2017-08-03 | Общество с ограниченной ответственностью "Геологический центр СПбГУ" (ООО "Геологический центр СПбГУ") | DEVICE FOR THREE-DIMENSIONAL SEISMIC AND ELECTROMAGNETIC PHYSICAL SIMULATION |
CN107422386A (en) * | 2017-04-20 | 2017-12-01 | 上海艾都能源科技有限公司 | A kind of intelligent electrical method physical prospecting instrument |
CN112630842A (en) * | 2020-10-14 | 2021-04-09 | 中国地质调查局南京地质调查中心(华东地质科技创新中心) | Dipole difference device and method applied to direct-current resistivity method exploration |
CN118549995A (en) * | 2024-07-29 | 2024-08-27 | 中国地质科学院地球物理地球化学勘查研究所 | Gravity data-based geological body boundary vertical information identification method and system |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5841280A (en) * | 1997-06-24 | 1998-11-24 | Western Atlas International, Inc. | Apparatus and method for combined acoustic and seismoelectric logging measurements |
RU2154847C1 (en) * | 1999-02-02 | 2000-08-20 | Нижне-Волжский научно-исследовательский институт геологии и геофизики | Process of geophysical prospecting while searching for oil and gas fields |
RU2155977C1 (en) * | 1999-05-19 | 2000-09-10 | Нижне-Волжский научно-исследовательский институт геологии и геофизики | Technique for geophysical prospecting |
RU37418U1 (en) * | 2004-01-28 | 2004-04-20 | Сочельников Виктор Васильевич | SYSTEM FOR SEISMOELECTRIC LOGGING |
RU50008U1 (en) * | 2005-05-13 | 2005-12-10 | Коноплев Михаил Юрьевич | SYSTEM FOR EXPLORATION OF OIL AND GAS DEPOSITS |
RU2301431C2 (en) * | 2005-03-24 | 2007-06-20 | Общество с ограниченной ответственностью "Сибирская геофизическая научно-производственная компания" | Mode of electrical exploration with using of spacing differentiation of the field of formation on several spreads |
RU2399931C2 (en) * | 2008-10-23 | 2010-09-20 | Общество с ограниченной ответственностью "Сибирская геофизическая научно-производственная компания" | Method for quantitative dividing effects of electromagnetic induction and induced polarisation |
-
2015
- 2015-03-03 RU RU2015107361/28A patent/RU2595327C1/en active IP Right Revival
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5841280A (en) * | 1997-06-24 | 1998-11-24 | Western Atlas International, Inc. | Apparatus and method for combined acoustic and seismoelectric logging measurements |
RU2154847C1 (en) * | 1999-02-02 | 2000-08-20 | Нижне-Волжский научно-исследовательский институт геологии и геофизики | Process of geophysical prospecting while searching for oil and gas fields |
RU2155977C1 (en) * | 1999-05-19 | 2000-09-10 | Нижне-Волжский научно-исследовательский институт геологии и геофизики | Technique for geophysical prospecting |
RU37418U1 (en) * | 2004-01-28 | 2004-04-20 | Сочельников Виктор Васильевич | SYSTEM FOR SEISMOELECTRIC LOGGING |
RU2301431C2 (en) * | 2005-03-24 | 2007-06-20 | Общество с ограниченной ответственностью "Сибирская геофизическая научно-производственная компания" | Mode of electrical exploration with using of spacing differentiation of the field of formation on several spreads |
RU50008U1 (en) * | 2005-05-13 | 2005-12-10 | Коноплев Михаил Юрьевич | SYSTEM FOR EXPLORATION OF OIL AND GAS DEPOSITS |
RU2399931C2 (en) * | 2008-10-23 | 2010-09-20 | Общество с ограниченной ответственностью "Сибирская геофизическая научно-производственная компания" | Method for quantitative dividing effects of electromagnetic induction and induced polarisation |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU172989U1 (en) * | 2016-12-28 | 2017-08-03 | Общество с ограниченной ответственностью "Геологический центр СПбГУ" (ООО "Геологический центр СПбГУ") | DEVICE FOR THREE-DIMENSIONAL SEISMIC AND ELECTROMAGNETIC PHYSICAL SIMULATION |
CN107422386A (en) * | 2017-04-20 | 2017-12-01 | 上海艾都能源科技有限公司 | A kind of intelligent electrical method physical prospecting instrument |
CN112630842A (en) * | 2020-10-14 | 2021-04-09 | 中国地质调查局南京地质调查中心(华东地质科技创新中心) | Dipole difference device and method applied to direct-current resistivity method exploration |
CN112630842B (en) * | 2020-10-14 | 2024-04-12 | 中国地质调查局南京地质调查中心(华东地质科技创新中心) | Dipole differential device and method applied to direct-current resistivity method exploration |
CN118549995A (en) * | 2024-07-29 | 2024-08-27 | 中国地质科学院地球物理地球化学勘查研究所 | Gravity data-based geological body boundary vertical information identification method and system |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101438176B (en) | Time lapse analysis with electromagnetic data | |
US7062072B2 (en) | Methods of producing images of underground formations surrounding a borehole | |
US8812237B2 (en) | Deep-reading electromagnetic data acquisition method | |
EP3039460B1 (en) | Mapping resistivity distribution within the earth | |
Sharlov et al. | Transient electromagnetic surveys for highresolution near-surface exploration: basics and case studies | |
EA033836B1 (en) | Systems for locating and imaging proppant in induced fractures | |
RU2595327C1 (en) | Device for seismic survey 2d or 3d, electrical exploration and gis to increase reliability of mapping the roof of salt and for prediction of oil and gas content subsalt deposits in areas with developed salt-dome tectonics | |
RU2722861C1 (en) | Static corrections calculation method | |
Riddle et al. | Subsurface tunnel detection using electrical resistivity tomography and seismic refraction tomography: A case study | |
CN104834006A (en) | Processing method and apparatus of structural cause deposition parameter diagram | |
CN106199708B (en) | A kind of city Geophysical Method for Prediction containing structure of coal and device | |
Hu et al. | Quantitative characterization of deep fault-karst carbonate reservoirs: A case study of the Yuejin block in the Tahe oilfield | |
RU2594112C2 (en) | System for exploration of oil and gas in complex-structure areas with developed salt-dome tectonics with mapping of roof of salt and subsalt deposits and computer-process system therefor | |
Cho | Stochastic discrete fracture network modeling in shale reservoirs via integration of seismic attributes and petrophysical data | |
RU2700836C1 (en) | Method of reservoirs saturation prediction based on complex analysis of data cpp, stm, gis | |
Commer et al. | Three-dimensional fracture continuum characterization aided by surface time-domain electromagnetics and hydrogeophysical joint inversion—proof-of-concept | |
Aminzadeh et al. | Fundamentals of Petroleum Geophysics | |
RU2540216C1 (en) | Method of forecast of capacitive parameters and type of fluid saturation of reservoirs | |
CN108072901B (en) | Method and system for obtaining accurate hydrostatic pressure and overlying formation pressure | |
Jahani et al. | Limits of 3D detectability and resolution of LWD deep-sensing borehole electromagnetic measurements acquired in the Norwegian Continental Shelf | |
Aktürk et al. | Integration of electrical resistivity imaging (ERI) and ground-penetrating radar (GPR) methods to identify soil profile around Necatibey Subway Station, Ankara, Turkey | |
CN109507741A (en) | Dual induction log fracture porosity analysis method and system | |
Ezersky et al. | Integrated study of the sinkhole development site on the Western shores of the Dead Sea using geophysical methods | |
Mooney et al. | Investigation of geoelectric-while-tunneling methods through numerical modeling | |
Xiao et al. | Characterizing fractured-vuggy carbonate reservoirs quantitatively based on the integration of well and seismic |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20190304 |
|
PD4A | Correction of name of patent owner | ||
NF4A | Reinstatement of patent |
Effective date: 20210119 |
|
PC43 | Official registration of the transfer of the exclusive right without contract for inventions |
Effective date: 20210125 |