RU2491580C1 - Method of measuring geophysical characteristics using successive inversion of geoelectric data with additional time filter - Google Patents
Method of measuring geophysical characteristics using successive inversion of geoelectric data with additional time filter Download PDFInfo
- Publication number
- RU2491580C1 RU2491580C1 RU2012104785/28A RU2012104785A RU2491580C1 RU 2491580 C1 RU2491580 C1 RU 2491580C1 RU 2012104785/28 A RU2012104785/28 A RU 2012104785/28A RU 2012104785 A RU2012104785 A RU 2012104785A RU 2491580 C1 RU2491580 C1 RU 2491580C1
- Authority
- RU
- Russia
- Prior art keywords
- geoelectric
- relaxation time
- filter
- time
- section
- Prior art date
Links
Images
Abstract
Description
Изобретение применяется в геологоразведке для поиска месторождений нефти и газа путем выделения аномальных зон вызванной поляризации. В частности, используют в дифференциально-нормированном методе электроразведки (ДНМЭ) с разделением полей ЕМ (электромагнитной индукции) и ВП (вызванной поляризации).The invention is used in geological exploration for the search for oil and gas deposits by identifying anomalous zones caused by polarization. In particular, they are used in the differential-normalized method of electrical exploration (DNME) with separation of the fields EM (electromagnetic induction) and VP (induced polarization).
Назначением предложенного способа является осуществление измерительных действий в процессе поиска углеводородов путем выделения аномальных зон вызванной поляризации (ВП) над залежами углеводородов (УВ).The purpose of the proposed method is the implementation of measuring actions in the process of searching for hydrocarbons by identifying anomalous zones of induced polarization (VP) over hydrocarbon deposits (HC).
Известно изобретение «Способ коррекции сейсмического разреза», заявка RU 2006139541, опубл. 20,05.2008, МПК G01V 1/28, в котором вводят статические поправки в сейсмические трассы, на сейсмических пикетах определяют временные сдвиги и осуществляют разделение гравитационного поля и аномалий потенциального поля методом частотной фильтрации, в результате чего получают временные аномалии разных порядков, которые преобразуют их во временные сдвиги. Однако данный метод не использует метод ДНМЭ с разделением полей ЕМ и ВП. Следовательно, не обеспечивает моделирование электромагнитных измерений в рамках горизонтально-слоистой поляризующейся среды, а также не учитывает появление аномалий над залежами пирита в зонах геохимических барьеров.The invention is known "Method of correction of seismic section", application RU 2006139541, publ. 20.05.2008, IPC
Известно изобретение «Способ контроля процесса гидроразрыва пласта залежи углеводородов», патент RU 2319177, опубл. 19.06.2006, МПК G01V 1/00, включающий обработку сейсмических сигналов в реальном времени в каждый дискретный момент времени и вычисляют координаты источников сейсмических колебаний для всех пар точек наблюдения. Способ позволяет контролировать процесс гидроразрыва пласта залежи углеводородов путем регистрации сейсмических сигналов. Однако определяет местоположение источника информационного каждого сигнала, а не в целом зоны микросейсмической активности, что обусловлено спектральным представлением сигналов. Способ не позволяет проводить фильтрацию исходного поля ВП, выделяя из него составляющую, связанную с узким диапазоном малых значений τ. Это, в свою очередь, не обеспечивает применения настройки временного фильтра на эталонный объект и не позволяет выделить отклик ВП, связанный с залежами УВ над пиритом.The invention is known "Method for controlling the process of hydraulic fracturing of hydrocarbon deposits", patent RU 2319177, publ. 06/19/2006, IPC
Известно изобретение «Сбор и фильтрация сейсмических данных (способ определения цифрового фильтра для сейсмических сигналов)», патент GB 0400409.9, патент RU 2364894 опубл. 20.08.2009, МПК G01V 1/36, в котором поэтапно определяют ограничения, представляющие фильтр для сохранения сигналов, используют итерационный процесс, при этом каждая итерация дополнительно содержит этапы, и налагают в указанной выборочной области вторые ограничения для получения итерированного фильтра. Способ позволяет при возможности разместить и изолировать волны-помехи в области преобразования, их удалить из преобразованных данных путем фильтрации или мьютинга, то есть путем задания значений в области, соответствующей волнам-помехам. Однако фильтры используют для ослабления шума, но коррелированный шум обычно не может быть исключен простым суммированием, а пространственно-временные способы работают по временному ряду зарегистрированных данных и не обеспечивают фильтрацию исходного поля ВП, выделяя из него составляющую, связанную с узким диапазоном малых значений τ. Иными словами, не используют возможность за счет применения настройки временного фильтра на эталонный объект, что позволяет выделить отклик ВП, связанный с залежами УВ над залежами пирита в зонах геохимических барьеров.The invention is known "Collection and filtering of seismic data (a method for determining a digital filter for seismic signals)", patent GB 0400409.9, patent RU 2364894 publ. 08/20/2009, IPC
Известно изобретение «Способ геофизической разведки и устройство для регистрации параметров естественного импульсного электромагнитного поля земли», заявка RU 2009100934, опубл. 20.07.2011, МПК G01V 3/08, в котором синхронные измерения интенсивности т в диапазоне частот используют для построения графиков пространственных изменений интенсивности полей вдоль профиля работ и затем настраивают чувствительность измерительных каналов в графики зависимостей, сглаживают изменения скользящим окном выбранной длительности, определяют вариации измеренных параметров вдоль профиля и дают геологическую интерпретацию полученных результатов при оценке структурных и литологических неоднородностей по изменению интенсивности сигнала. Однако при данном способе уточнять полученные данные можно только последующим бурением выделенных аномальных зон, что очень дорого. Кроме того, способ не обеспечивает настройки временного фильтра на эталонный объект, что позволяет выделить отклик ВП, связанный с залежами УВ над залежами пирита в зонах геохимических барьеров. Следовательно, не позволяет применить настройки временного фильтра на эталонный объект.The invention is known "Method of geophysical exploration and a device for recording parameters of the natural pulsed electromagnetic field of the earth", application RU 2009100934, publ. 07/20/2011, IPC
Известно изобретение «Способ устранения искажающего влияния верхней части разреза в сейсморазведке», патент RU 2381529, опубл. 10.02.2010, МПК G01V 1/36, G01V 1/28, включающий использование эталонных величин амплитуд прямой волны в разных породах относительно глубины залегания кровли коренных пород, определяют коэффициент затухания упругих колебаний, рассчитывают запаздывания времен, формируют импульсную характеристику фильтра, описанного последовательностью взятых с учетом знаков единичных импульсов с амплитудами, распределенных во времени. Способ позволяет: более полно устранить искажающее влияние верхней части разреза на данные сейсморазведки, и повысить точности измерения и интерпретации. Однако в зарегистрированные сейсмозаписи вводят равные статическим поправкам временные сдвиги, что на малых расстояниях затрудняет прослеживание прямой волны с повышенной частотой и неустойчивой полярностью. Не позволяет уменьшить влияние эквивалентных зависимостей между расчетными поляризационными параметрами, и, следовательно, повысить надежность прогноза. Не позволяет выделить аномалии ВП, связанные с невскрытыми бурением залежами УВ, если таковые имеются на изучаемой площади, а, следовательно, не позволяет за счет применения настройки временного фильтра на эталонный объект выделить отклик ВП, связанный с залежами УВ над залежами пирита.The invention is known "Method of eliminating the distorting effect of the upper part of the section in seismic exploration", patent RU 2381529, publ. 02/10/2010, IPC
Наиболее близким к предложенному способу является изобретение «Сбор и фильтрация сейсмических данных», патент RU 2364894, опубл. 20.08.2009, GB 0400409.9, МПК G01V 1/36, включающий осуществление возбуждения и регистрацию в среде неустановившегося электромагнитного поля, измерение минимального времени регистрации на стадиях, соответствующих однократным возмущениям ЭМ поля, выявление интервала времен и обработку, которую осуществляют путем выделения разностных сигналов в выявленном интервале времен, и последующее построение пространственных годографов. Однако, при данном способе имеет место слабая дифференциация разреза на отдельные слои, то есть невозможно расчленить весь разрез на отдельные слои при зондировании мощных геологических толщ. Данный способ обеспечивает невысокую точность результатов и разрешающую способность, которая увеличивается за счет осреднения геоэлектрических параметров в пределах выделяемых толщин из-за больших, без интегрирования, сигналов в пределах расстановок приемников поля. А, следовательно, не обеспечивает за счет применения настройки временного фильтра на эталонный объект выделения отклика ВП, связанного с залежами УВ над залежами пирита. Не обеспечивает устойчивости получаемых результатов при изменении входных данных, и не уменьшает влияние эквивалентных зависимостей между расчетными поляризационными параметрами, т.е. не обеспечивает повышение надежности прогноза.Closest to the proposed method is the invention "the Collection and filtering of seismic data", patent RU 2364894, publ. 08/20/2009, GB 0400409.9, IPC
Одной из наиболее существенных причин появления аномалий вызванной поляризации (ВП) над залежами углеводородов (УВ) является наличие эпигенетического пирита в верхней части разреза в зонах геохимических барьеров, обусловленное наличием восстановительной обстановки над залежами УВ, что является известным фактом (Моисеев B.C. Метод вызванной поляризации при поиске нефтеперспективных площадей. Новосибирск: Наука, 2000. - 136 с.; Моисеев B.C., Тараторкин Б.Ф., Шлепнев В.Б. Результативность прогноза залежей углеводородов методом вызванной поляризации в Западной Сибири. Тез. Докл. - Международная геофизическая конференция, С.-Петербург, 1995; Комаров В.А. Электроразведка методом вызванной поляризации. 2-е изд., перераб. и доп. «Недра», Ленинград, 1980 - 391 с.).One of the most significant reasons for the appearance of induced polarization (VP) anomalies over hydrocarbon (HC) deposits is the presence of epigenetic pyrite in the upper part of the section in zones of geochemical barriers, due to the presence of a reducing environment over hydrocarbon deposits, which is a known fact (Moiseev BC The method of induced polarization at Search for Oil-Prospective Areas, Novosibirsk: Nauka, 2000. - 136 pp .; Moiseev BC, Taratorkin BF, Shlepnev VB Efficiency of forecasting hydrocarbon deposits by the method of induced polarization and in Western Siberia, Abstracts: International Geophysical Conference, St. Petersburg, 1995; Komarov, V.A., Electrical Exploration by Evoked Polarization Method, 2nd ed., revised and supplemented "Nedra", Leningrad, 1980 - 391 p.).
Эпигенетический пирит характеризуется электронной проводимостью и преимущественно является тонкозернистым (менее 0.1 мм), т.е. присутствуют эпигенетические тонкозернистые вкрапления пирита в зонах геохимических барьеров над залежами углеводородов.Epigenetic pyrite is characterized by electronic conductivity and is predominantly fine-grained (less than 0.1 mm), i.e. epigenetic fine-grained pyrite inclusions are present in zones of geochemical barriers over hydrocarbon deposits.
Процесс вызванной поляризации электронных проводников связан с различными типами перенапряжения - перенапряжением разряда, перенапряжением адсорбции и перенапряжением диффузии. Опытным путем для образцов пирита из полиметаллической руды было установлено время релаксации τ (синоним - постоянная времени) для данных типов перенапряжения: для перенапряжения разряда τ составляет 0.002-0.01 с, для перенапряжения адсорбции - 0.1-0.2 с, для перенапряжения диффузии - от единиц секунд до нескольких минут (Комаров, 1980). Поскольку в ДНМЭ время спада, в течение которого производятся замеры DU, составляет 2-8 с, значимыми остаются два первых типа перенапряжения.The process of induced polarization of electronic conductors is associated with various types of overvoltage — discharge overvoltage, adsorption overvoltage, and diffusion overvoltage. Experimentally, for pyrite samples from polymetallic ore, the relaxation time τ (synonym - time constant) was established for these types of overvoltage: for discharge overvoltage, τ is 0.002-0.01 s, for adsorption overvoltage - 0.1-0.2 s, for diffusion overvoltage - from units of seconds up to several minutes (Komarov, 1980). Since in DNME the decay time during which DU measurements are made is 2–8 s, the first two types of overvoltage remain significant.
Кроме того, по данным В.А. Комарова, при постоянном токе время достижения максимального значения вызванной поляризации для больших сфер асимптотически приближается к постоянной времени τ и не зависит от радиуса, в то время как для маленьких электронопроводящих сфер оно пропорционально радиусу сферы и будет меньше постоянной времени τ. Эпигенетический пирит является тонкозернистым (менее 0.1 мм). Следовательно, время релаксации, связанное с переходными процессами в перекрывающей залежь УВ толще, включающей эпигенетический пирит, будет не выше, чем определенное опытным путем для образцов пирита (т.е. десятые доли секунды и меньше). Таким образом, диапазон выбора значений времени релаксации значительно сужается и находится в области малых значений. Проводить фильтрацию исходного поля ВП, выделяя из него составляющую, связанную с узким диапазоном малых значений τ, необходимо с опорой на эталонные объекты - продуктивные и непродуктивные скважины.In addition, according to V.A. Komarova, at constant current the time to reach the maximum value of the induced polarization for large spheres asymptotically approaches the time constant τ and does not depend on the radius, while for small electronically conducting spheres it is proportional to the radius of the sphere and will be less than the time constant τ. Epigenetic pyrite is fine-grained (less than 0.1 mm). Consequently, the relaxation time associated with the transient processes in the thicker overlapping hydrocarbon deposit, including epigenetic pyrite, will not be higher than determined experimentally for pyrite samples (i.e., tenths of a second or less). Thus, the range of choice of relaxation time values is significantly narrowed and lies in the region of small values. It is necessary to carry out the filtration of the initial field of the IP, extracting from it a component associated with a narrow range of small values of τ, relying on reference objects - productive and unproductive wells.
Пределы изменения времени релаксации τ будут зависеть от количества пирита и состава пород, содержащих пирит, следовательно, для разных площадей они могут отличаться.The limits of variation of the relaxation time τ will depend on the amount of pyrite and the composition of the rocks containing pyrite; therefore, for different areas, they may differ.
Поле становления, замеряемое в ходе полевых работ, включает поля разной природы: поле вызванной поляризации (ВП), обусловленное гальваническими процессами (поле IPg), и поле электромагнитной индукции, вызванное вихревыми токами (поле ЕМ), при этом поля ЕМ и ВП не являются аддитивными. Изучение поля ЕМ позволяет определять послойное распределение удельного электрического сопротивления (проводимости), а, следовательно, определять литологические характеристики разреза. При этом в методе ДНМЭ поиски месторождений нефти и газа осуществляют путем выделения аномальных зон вызванной поляризации, приуроченных к верхней части разреза и обусловленных наличием восстановительной обстановки над скоплениями углеводородов. Аномальное повышение поля ВП над залежами УВ связано в значительной мере с наличием эпигенетического тонкозернистого вкрапленного пирита в зонах геохимических барьеров (в районе расположения покрышки над залежью УВ и верхнего регионального водоупора). Это не позволяет без использования дополнительных временных фильтров наиболее полно выделить аномальные зоны, указывающие на залежи УВ над этими зонами, а, следовательно, возникают «слепые» зоны, где наличие УВ не определено, что существенно снижает эффективность геологической разведки.The field of formation, measured during field work, includes fields of a different nature: the field of induced polarization (VP) due to galvanic processes (field IPg), and the field of electromagnetic induction caused by eddy currents (field EM), while the fields EM and VP are not additive. Studying the EM field allows one to determine the layered distribution of electrical resistivity (conductivity), and, therefore, determine the lithological characteristics of the section. Moreover, in the DNME method, the search for oil and gas deposits is carried out by identifying anomalous zones of induced polarization, confined to the upper part of the section and due to the presence of a reducing environment over hydrocarbon accumulations. An abnormal increase in the VP field over hydrocarbon deposits is associated largely with the presence of epigenetic fine-grained disseminated pyrite in the zones of geochemical barriers (in the region where the tire is located above the hydrocarbon pool and the upper regional confinement). This does not allow without the use of additional time filters to most fully identify the anomalous zones indicating hydrocarbon deposits above these zones, and, therefore, “blind” zones arise where the presence of hydrocarbons is not determined, which significantly reduces the effectiveness of geological exploration.
В условиях терригенного разреза с хорошей проводимостью (сотни См и выше), величина индуктивной составляющей уменьшается медленно в течение спада ВП, что затрудняет количественное разделение индукционного поля ЕМ и поля вызванной поляризации гальванического происхождения ВП. В общепринятой методике инверсии данных с использованием формулы Cole-Cole и в методике ДНМЭ с разделением полей IPg и ЕМ, поляризационные параметры рассчитываются в широких пределах. В условиях хорошо проводящего разреза, исходя из изменения соотношения полей ЕМ и ВП в течение времени спада, необходимо использовать дополнительные временные фильтры, которые бы позволили выделить полезный сигнал ВП, позволяющий устанавливать наличие углеводородов.Under the conditions of a terrigenous section with good conductivity (hundreds of S and higher), the magnitude of the inductive component decreases slowly during the decay of the VP, which makes it difficult to quantify the separation of the EM induction field and the field caused by polarization of the galvanic origin of the VP. In the generally accepted method of data inversion using the Cole-Cole formula and in the DNME method with separation of the IPg and EM fields, polarization parameters are calculated over a wide range. In conditions of a well-conducting section, based on the change in the ratio of the EM and VP fields during the recession time, it is necessary to use additional time filters that would allow us to select a useful VP signal, which allows us to determine the presence of hydrocarbons.
Сущностью предложенного способа является то, что осуществление измерения процесса становления над поляризующейся средой и дальнейшее моделирование электромагнитных измерений в рамках горизонтально-слоистой поляризующейся среды осуществляют с использованием временных фильтров. Настройки фильтров по времени релаксации τ проводятся на точках наблюдения ДНМЭ в районе эталонных геофизических объектов. Такая настройка по τ позволяет выделить аномалии ВП, связанные с невскрытыми бурением залежами УВ, если таковые имеются на изучаемой площади.The essence of the proposed method is that the measurement of the formation process over a polarizing medium and further modeling of electromagnetic measurements within a horizontally layered polarized medium is carried out using temporary filters. The filter settings for the relaxation time τ are made at the observation points of the DNME in the region of the reference geophysical objects. Such a tuning with respect to τ makes it possible to identify airspace anomalies associated with undiscovered drilling of hydrocarbon deposits, if any, on the study area.
Технический результат, достигаемый в результате предложенного способа состоит в расширении технических средств измерения за счет применения настройки временного фильтра на эталонный объект, что позволяет выделить отклик ВП, связанный с залежами УВ, обеспечить устойчивость получаемых результатов при изменении входных данных, уменьшить влияние эквивалентных зависимостей между расчетными поляризационными параметрами, и, следовательно, повысить надежность прогноза.The technical result achieved as a result of the proposed method consists in expanding the technical means of measurement by applying the time filter setting to the reference object, which allows one to distinguish the response of the EP associated with hydrocarbon deposits, to ensure the stability of the results obtained when the input data changes, to reduce the effect of equivalent dependencies between the calculated polarization parameters, and, therefore, increase the reliability of the forecast.
Данный технический результат достигается за счет того, что в способе измерения геофизических характеристик, применяя последующую инверсию геоэлектрических данных с дополнительным временным фильтром, используют определение послойного распределения удельного электрического сопротивления (проводимости) геологического разреза, после чего строят геоэлектрическую характеристику разреза с использованием временных фильтров, получая геоэлектрическую модель разреза. Предложенный способ отличается тем, что при построении геоэлектрической характеристики разреза вводят дополнительные временные фильтры для времени релаксации τ, предварительно осуществляют выбор количества слоев и пределов изменения их толщин в выбранной геоэлектрической модели на основе имеющейся априорной информации о глубинах расположения стратиграфических горизонтов, полученных по результатам глубокого бурения, и о глубинах расположения отражающих горизонтов, полученных на основе данных сейсморазведки, затем осуществляют определение пределов изменения удельного электрического сопротивления в каждом измеряемом слое геоэлектрической модели на основе имеющихся данных каротажа, послойно выбирают пределы времени релаксации τ для осуществления инверсии таким образом, чтобы они позволяли получить аномальный отклик ВП в районе выбранного эталонного объекта. При этом по меньшей мере в одном геоэлектрическом слое выбирают диапазон с низкими значениями τ (менее 0.1 с), после чего настраивают измерительную систему на диапазон времени релаксации τ, равный эталонному интервалу времени релаксации τ таким образом, чтобы минимальное и максимальное значения τ в выбранном диапазоне отличались менее чем в 10 раз; осуществляют инверсию геоэлектрических данных путем подбора измеренных с помощью ДНМЭ кривых дифференциально-нормированных параметров, характеризующих поляризационные процессы в геологической среде по следующим характеристикам: коэффициент поляризуемости η, время релаксации τ, ширина релаксационного спектра с и осуществляют проверку временного фильтра в районе остальных эталонных геофизических объектов. При этом если в районе эталонного объекта имеется продуктивная скважина, то при использовании данного фильтра должен быть получен аномальный отклик ВП, если в районе эталонного объекта имеется непродуктивная скважина, должен быть получен только фоновый отклик ВП. Если с помощью дополнительного временного фильтра выявлены все аномалии ВП над известными залежами на исследуемой площади, и дополнительный временной фильтр не создает ложных аномалий там, где достоверно известно, что залежей нет, делают вывод о работоспособности выбранного временного фильтра. Если выбранный фильтр не позволяет получить данные, соответствующие априорной информации, выбранные пределы постоянной времени τ корректируют и проводят повторную проверку. После проверки работоспособности фильтра на эталонных объектах осуществляют моделирование геологического разреза методом ДНМЭ по всей исследуемой площади. В частном случае для выделения полезного сигнала ВП по времени релаксации τ используют пределы времени релаксации τ с применением двух видов дополнительных фильтров τ:This technical result is achieved due to the fact that in the method of measuring geophysical characteristics, using the subsequent inversion of geoelectric data with an additional time filter, the determination of the layer-by-layer distribution of the electrical resistivity (conductivity) of the geological section is used, after which the geoelectric characteristic of the section is constructed using temporary filters, obtaining geoelectric sectional model. The proposed method is characterized in that when constructing the geoelectric characteristics of the section, additional time filters are introduced for the relaxation time τ, and the number of layers and the limits of changes in their thickness are preliminarily selected in the selected geoelectric model based on a priori information on the depths of the stratigraphic horizons obtained from deep drilling , and the depths of the reflecting horizons obtained on the basis of seismic data, then Determining the limits of changes in the electrical resistivity in each measured layer of the geoelectrical model based on the available logging data, the limits of the relaxation time τ for the inversion are selected in layers so that they allow an abnormal response of the VP in the region of the selected reference object. In this case, a range with low values of τ (less than 0.1 s) is selected in at least one geoelectric layer, after which the measuring system is tuned to the range of relaxation time τ equal to the reference interval of relaxation time τ so that the minimum and maximum values of τ in the selected range differed less than 10 times; the geoelectric data are inverted by selecting the differential-normalized parameters curves measured using DNME that characterize the polarization processes in the geological environment according to the following characteristics: polarizability coefficient η, relaxation time τ, width of the relaxation spectrum s and check the time filter in the region of other reference geophysical objects. At the same time, if there is a productive well in the area of the reference object, then using this filter an abnormal response of the EW should be obtained, if in the area of the reference object there is an unproductive well, only the background response of the EAP should be obtained. If with the help of an additional time filter all the anomalies of the airspace over known deposits in the studied area are detected, and the additional time filter does not create false anomalies where it is reliably known that there are no deposits, they conclude that the selected temporary filter is working. If the selected filter does not allow obtaining data corresponding to a priori information, the selected limits of the time constant τ are corrected and a repeated check is carried out. After checking the filter's operability at the reference objects, the geological section is modeled by the DNME method over the entire studied area. In a particular case, to extract a useful VP signal from the relaxation time τ, the relaxation time limits τ are used using two types of additional filters τ:
- - - с относительно высокими значениями τ=0,1-5 сек- - - with relatively high values of τ = 0.1-5 sec
- - - с относительно низкими значениями τ=0,01-0,1 сек.- - - with relatively low values of τ = 0.01-0.1 sec.
Чертежи демонстрируют и поясняют предложенный способ, однако не охватывают всех частных случаев предложенного способа.The drawings demonstrate and explain the proposed method, but do not cover all the particular cases of the proposed method.
На Фиг.1 - показано распределение эпигенетической пиритизации над залежью Северо-Гуляевского месторождения.Figure 1 - shows the distribution of epigenetic pyritization over the reservoir of the North-Gulyaevskoye field.
На Фиг.2 - показаны теоретические расчеты концентрации пирита на уровне регионального водоупора при наличии залежи УВ и без нее (Северо-Гуляевское месторождение).Figure 2 - shows theoretical calculations of the concentration of pyrite at the level of regional confinement in the presence of hydrocarbon deposits and without it (Severo-Gulyaevskoye field).
На Фиг.3 - показана временная зависимость ВП образца пирита из полиметаллической руды.Figure 3 - shows the time dependence of the VP of a sample of pyrite from polymetallic ore.
На Фиг.4 - показано сопоставление распределения поля поляризуемости, полученного при использовании различных временных фильтров (район Причерноморско-Северо-Кавказской нефтегазоносной провинции (НГП))Figure 4 shows a comparison of the distribution of the field of polarizability obtained using various time filters (the region of the Black Sea-North-Caucasian oil and gas province (NGP))
На Фиг.5 - показано сопоставление распределения отклика вызванной поляризации (охарактеризован комплексным поляризационным параметром N) с данными сейсморазведки и глубокого бурения, демонстрирующее достаточно хорошее совпадение планового положения аномалий ВП с контурами нефтяных залежей (район Волго-Уральской НГП).Figure 5 shows a comparison of the response distribution of the induced polarization (characterized by the complex polarization parameter N) with seismic and deep drilling data, which demonstrates a fairly good agreement between the planned position of the airspace anomalies and the contours of the oil deposits (the Volga-Ural oil and gas field).
На Фиг.6 - показано сопоставление модельных кривых, полученных для разных временных фильтров, демонстрирующее зависимость полученного отклика ВП от выбранных фильтров постоянной времени τ (район Волго-Уральской НГП).Figure 6 shows a comparison of the model curves obtained for different time filters, demonstrating the dependence of the obtained VP response on the selected filters with a time constant τ (region of the Volga-Ural oil and gas condensate field).
На Фиг.7 - показан отклик ВП в районе структурного поднятия, полученный по данным ДНМЭ на близко расположенных профилях без использования временного фильтра.Figure 7 - shows the response of the VP in the region of the structural uplift, obtained according to DNME on closely spaced profiles without the use of a time filter.
На Фиг.8 - показан отклик ВП в районе того же структурного поднятия, что и на Фиг.7, полученный по данным ДНМЭ на близко расположенных профилях с использованием временного фильтра.On Fig - shows the response of the VP in the region of the same structural uplift as in Fig.7, obtained according to DNME on closely spaced profiles using a time filter.
Предложенный способ реализовывают следующим образом.The proposed method is implemented as follows.
Для осуществления измерений требуется учитывать данные о наличии известных месторождений УВ. При этом, применяя метод измерения ДНМЭ, следует учитывать, что:To perform measurements, it is necessary to take into account data on the presence of known hydrocarbon deposits. In this case, using the method of measuring DNME, it should be borne in mind that:
τ - время релаксации - которая есть суть относительная скорость протекания электрохимических реакций. Поскольку широкий диапазон используемых в ходе инверсии значений времени релаксации (τ=0,01-5 с) - не обеспечивает устойчивость получаемых результатов (определено опытным путем), то надо его сузить.τ - relaxation time - which is the essence of the relative rate of electrochemical reactions. Since the wide range of relaxation times used during the inversion (τ = 0.01-5 s) does not ensure the stability of the results obtained (determined experimentally), it is necessary to narrow it down.
Из общего диапазона можно выбрать два диапазона - диапазон относительно высоких значений (τ=0,1-5 с) и диапазон относительно низких значений (τ=0,01-0,1 с).Two ranges can be selected from the general range — the range of relatively high values (τ = 0.1–5 s) and the range of relatively low values (τ = 0.01–0.1 s).
Первый диапазон - диапазон относительно высоких значений времени релаксации (τ=0,1-5 сек) - не позволяет выделить все известные нефтяные месторождения на заданной изучаемой площади.The first range - the range of relatively high values of the relaxation time (τ = 0.1-5 sec) - does not allow you to select all known oil fields in a given study area.
Второй диапазон - диапазон относительно низких значений времени релаксации (τ=0,01-0,1 сек.) - позволяет получить аномальный отклик ВП над районами известных, т.е. эталонных залежей УВ, и выделить новые перспективные участки (получение аномального отклика над известными месторождениями является показателем работоспособности данного фильтра и позволяет говорить о надежности прогноза новых перспективных участков).The second range — the range of relatively low values of the relaxation time (τ = 0.01–0.1 sec.) –– allows one to obtain an abnormal response of VP over known regions, i.e. reference hydrocarbon deposits, and identify new promising areas (obtaining an abnormal response over known deposits is an indicator of the performance of this filter and allows us to talk about the reliability of the forecast of new promising areas).
Следовательно, правильно выбранный для проведения инверсии диапазон τ позволяет по всей площади моделируемого геологического участка получить аномальный отклик ВП над залежами углеводородов, связанный, главным образом, с присутствием тонкозернистого пирита (доли мм) в зонах геохимических барьеров. Для решения поставленной задачи необходимо выбрать узкий диапазон τ, изменения τ в котором не превышают одного порядка, т.е. минимальное значение τ в выбранном диапазоне отличается от максимального значения τ менее, чем в 10 раз.Therefore, the correctly selected range of inversions τ allows one to obtain an abnormal response of VP over hydrocarbon deposits over the entire area of the simulated geological section, mainly associated with the presence of fine-grained pyrite (fraction of mm) in the zones of geochemical barriers. To solve this problem, it is necessary to choose a narrow range of τ, the changes of τ in which do not exceed one order, i.e. the minimum value of τ in the selected range differs from the maximum value of τ by less than 10 times.
Последовательность измерений при этом следующая:The measurement sequence is as follows:
- сначала выбираются пределы изменения мощности слоев и их сопротивления на основе данных каротажа и сейсморазведки- first, the limits of the thickness variation of the layers and their resistance are selected based on the logging and seismic data
- потом послойно выбирают узкий интервал времени релаксации τ- then a narrow relaxation time interval τ is selected in layers
- затем осуществляют проверку измерений в районе эталона и в его окрестностях- then carry out verification of measurements in the area of the standard and in its vicinity
- если результат не может быть геологически обоснован, тогда выбирают другие узкие интервалы времени релаксации и повторяют проверку- if the result cannot be geologically substantiated, then select other narrow relaxation time intervals and repeat the test
- если результат геологически обоснован и соответствует априорной информации, тогда осуществляют моделирование по всей площади- if the result is geologically justified and corresponds to a priori information, then modeling is performed over the entire area
Аномальный отклик ВП над залежами УВ в значительной мере обусловлен наличием тонкозернистого пирита, характеризующегося электронной проводимостью. Связь аномального отклика ВП с наличием эпигенетического пирита над залежами УВ подтверждается опытным путем. Для примера на Фиг.1 и 2 показано распределение эпигенетической пиритизации над залежью Северо-Гуляевского месторождения и показаны теоретические расчеты концентрации пирита на уровне регионального водоупора при наличии залежи УВ и без нее. На Фиг.3 (взят из монографии В.А. Комарова, 1980) показаны временные зависимости отклика ВП образца пирита из полиметаллической руды. Из данного графика видно, что при постоянном токе время достижения максимального значения вызванной поляризации Tmax для маленьких электронопроводящих сфер (С<<1) пропорционально радиусу сферы, в то время как для больших сфер (С>>1) асимптотически приближается к постоянной времени τ (обозначена как Т0) и не зависит от радиуса. Процесс вызванной поляризации электронных проводников связан:The abnormal response of VP over hydrocarbon deposits is largely due to the presence of fine-grained pyrite, characterized by electronic conductivity. The relationship of the abnormal response of VP with the presence of epigenetic pyrite over hydrocarbon deposits is confirmed experimentally. For example, Figs. 1 and 2 show the distribution of epigenetic pyritization over the reservoir of the Severo-Gulyaevskoye deposit and show theoretical calculations of the concentration of pyrite at the level of regional confinement in the presence of hydrocarbon deposits and without it. Figure 3 (taken from a monograph by V.A. Komarov, 1980) shows the time dependence of the response of the VP of a pyrite sample from polymetallic ore. It can be seen from this graph that, at constant current, the time to reach the maximum value of the induced polarization T max for small electronically conducting spheres (C << 1) is proportional to the radius of the sphere, while for large spheres (C >> 1) it asymptotically approaches the time constant τ (designated as T 0 ) and is independent of radius. The process of induced polarization of electronic conductors is related:
с перенапряжением разряда (время релаксации 0.002-0.01 с);with discharge overvoltage (relaxation time 0.002-0.01 s);
с перенапряжением адсорбции (время релаксации 0.1-0.2 с);with adsorption overvoltage (relaxation time 0.1-0.2 s);
с перенапряжением диффузии (время релаксации от единиц секунд до нескольких минут).with diffusion overvoltage (relaxation time from units of seconds to several minutes).
Хорошо известен тот факт, что скорость спада ВП для маленьких электронопроводящих сфер пропорциональна радиусу сферы (Комаров В.А. Электроразведка методом вызванной поляризации. 2-е изд., перераб. и доп. Л., Недра, 1980, 391 с.).It is well known that the VP decay rate for small electronically conducting spheres is proportional to the radius of the sphere (Komarov V.A. Electrical exploration by the method of induced polarization. 2nd ed., Revised and additional L., Nedra, 1980, 391 pp.).
Для выбора необходимых пределов времени релаксации τ, позволяющих установить аномальный отклик над скоплениями УВ в условиях хорошо проводящего разреза, для одного из районов Причерноморско-Северо-Кавказской НГП, где имеется ряд известных нефтяных месторождений, были проведены расчеты с использованием трех различных временных фильтров: с широким диапазоном изменения значений времени релаксации τ (0.01-5 с), с относительно высокими значениями τ (0.1-5 с) и с относительно низкими значениями τ (0.01-0.1 с) (Фиг.4). На рисунке 4 представлены временные разрезы поля поляризуемости IPg по профилю ДНМЭ, пересекающему два месторождения нефти, которые приурочены к структурным поднятиям и подтверждены результатами бурения. По геологическим данным предполагается, что зона сброса не является перспективной на поиски УВ. Установлено, что использование широкого диапазона значений τ (Фиг.4а) не обеспечивает устойчивости получаемых результатов, и при изменении входных данных в выбранных пределах может получиться прямо противоположный результат. Применение относительно высоких значений времени релаксации не позволяет выделить все известные на данной площади нефтяные месторождения (Фиг.4б). Оптимальным для решения поставленной задачи является использование малых значений τ. При этом, опираясь на данные о местоположении известных месторождений и с учетом необходимости устойчивости получаемого результата, методом подбора диапазон значений τ был максимально сужен (Фиг.4в). Выбранная модель позволила получить аномальный отклик ВП над всеми известными в районе залежами нефти, а также выделить новые участки, перспективные на поиски УВ.To select the necessary limits of the relaxation time τ, which allow us to establish an anomalous response over hydrocarbon accumulations in a well-conducting section, for one of the regions of the Black Sea-North-Caucasian oil and gas field, where there are a number of known oil fields, calculations were performed using three different time filters: a wide range of changes in the values of the relaxation time τ (0.01-5 s), with relatively high values of τ (0.1-5 s) and with relatively low values of τ (0.01-0.1 s) (Figure 4). Figure 4 shows temporary sections of the polarizability field IPg along the DNME profile intersecting two oil fields, which are confined to structural uplifts and confirmed by drilling results. According to geological data, it is assumed that the discharge zone is not promising for the search for hydrocarbons. It was found that the use of a wide range of values of τ (Fig. 4a) does not ensure the stability of the results obtained, and when the input data changes within the selected limits, the opposite result can be obtained. The use of relatively high values of the relaxation time does not allow to isolate all known in this area of oil fields (Fig.4b). The best way to solve this problem is to use small values of τ. Moreover, relying on data on the location of known deposits and taking into account the need for stability of the result, the range of τ values was narrowed as much as possible using the selection method (Fig. 4c). The selected model made it possible to obtain an abnormal response of airspace over all known oil deposits in the region, as well as to identify new areas promising for the search for hydrocarbons.
Метод использования временных фильтров был апробирован на одной из площадей в пределах Волго-Уральской НГП, в районе известного месторождения. Здесь также был выбран узкий диапазон значений времени релаксации, причем разный в разных геоэлектрических слоях. Кроме того, были сужены пределы изменения значений ширины релаксационного спектра, что позволило дополнительно уменьшить влияние эквивалентных зависимостей между рассчитываемыми поляризационными параметрами. Контуры аномалий ВП по большей части совпали с контурами ВНК, выделенными по данным бурения и сейсморазведки. При этом полученные данные позволили уточнить прогноз для проектных скважин (Фиг.5). Применение значений τ более 0.1 с, как и на первой площади, приводит к получению иного результата (Фиг.6) и не позволяет выделить все известные залежи УВ на данной территории.The method of using temporary filters was tested on one of the areas within the Volga-Ural oil and gas field, in the region of a well-known field. A narrow range of relaxation times was also chosen here, and it is different in different geoelectric layers. In addition, the limits of the variation in the width of the relaxation spectrum were narrowed, which made it possible to further reduce the effect of equivalent dependences between the calculated polarization parameters. The contours of the airspace anomalies for the most part coincided with the contours of the KSS, identified according to drilling and seismic data. At the same time, the data obtained made it possible to clarify the forecast for design wells (Figure 5). The use of values of τ more than 0.1 s, as in the first area, leads to a different result (Fig.6) and does not allow you to select all the known hydrocarbon deposits in this territory.
На Фиг.7 показано распределение поляризационного отклика в районе структурного поднятия, полученное без использования временного фильтра, демонстрирующего наличие аномалии ВП только на небольшом участке на краю структуры (показан одни из участков, расположенных в Причерноморско-Северо-Кавказской НГП). На момент проведения работ ДНМЭ имелись данные только по сейсморазведочным работам 2D, не обеспечивающим надежные контуры структуры. Позже на данном участке были проведены сейсморазведочные работы 3D. Повторные работы ДНМЭ, уже с использованием временных фильтров, позволили получить аномальный отклик в районе данной структуры, подтвержденный данными бурения (Фиг.8).Figure 7 shows the distribution of the polarization response in the region of structural elevation, obtained without the use of a temporary filter, demonstrating the presence of VP anomalies only in a small area at the edge of the structure (one of the areas located in the Black Sea-North Caucasus OGP is shown). At the time of the DNME work, only 2D seismic data were available that did not provide reliable structure contours. Later, 3D seismic surveys were performed at this site. Repeated work of DNME, already using temporary filters, allowed to obtain an abnormal response in the region of this structure, confirmed by drilling data (Fig. 8).
Таким образом, использование временных фильтров, подобранных для конкретных геолого-геофизических условий, повышает надежность прогноза нефтегазоперспективных участков.Thus, the use of temporary filters selected for specific geological and geophysical conditions increases the reliability of the forecast of oil and gas prospective areas.
Выбор узких диапазонов времени релаксации осуществляется путем произвольного перебора разных диапазонов относительно низких значений τ отдельно для каждого слоя модели. Под произвольным перебором понимается следующее действие: берется произвольный узкий диапазон низких значений τ, производится моделирование, определение правильности выбора основано на сравнении полученных данных с априорной информацией. Если результат удовлетворяет, данный диапазон используется для дальнейших расчетов, если результат не удовлетворяет, берется другой диапазон низких значений τ, производятся повторные расчеты и сравнение с априорной информацией.The narrow ranges of relaxation time are selected by randomly sorting different ranges of relatively low values of τ separately for each layer of the model. By arbitrary enumeration, the following action is understood: an arbitrary narrow range of low values of τ is taken, modeling is carried out, determination of the correctness of the choice is based on a comparison of the data obtained with a priori information. If the result is satisfactory, this range is used for further calculations, if the result is not satisfactory, another range of low values of τ is taken, repeated calculations and comparison with a priori information are performed.
Правильность выбора значений τ определяется:The correctness of the choice of values of τ is determined by:
1) хорошей сходимостью полевых и модельных кривых, при этом погрешность сходимости не должна превышать допустимую погрешность измерения, которая составляет 5% (Инструкция по электроразведке: Наземная электроразведка, скважинная электроразведка, шахтно-рудничная электроразведка, аэроэлектроразведка, морская электроразведка / М-во геологии СССР. - Л.: Недра, 1984. - 352 с., страница 19, пункт 3.3.1.19);1) good convergence of field and model curves, while the convergence error should not exceed the permissible measurement error, which is 5% (Instructions for electrical exploration: electrical exploration, downhole electrical exploration, mine-electrical exploration, aerial exploration, marine electrical exploration / M-geology of the USSR . - L .: Nedra, 1984. - 352 p.,
2) получением аномального отклика ВП над достоверно известными залежами УВ (подтвержденными бурением).2) obtaining an abnormal response of airspace over reliably known hydrocarbon deposits (confirmed by drilling).
Далее осуществляется моделирование на всей площади, и, с учетом сопоставления с распределением удельного электрического сопротивления, полученного по результатам ДНМЭ, а также с имеющейся априорной информацией о геолого-геофизическом строении района и положении известных на основе сейсморазведки и данных бурения ловушек разного типа, делают вывод о наличии или отсутствии аномального отклика ВП, связанного с предполагаемыми залежами УВ.Next, modeling is carried out over the entire area, and, taking into account the comparison with the distribution of electrical resistivity obtained from the results of DNME, as well as with the available a priori information about the geological and geophysical structure of the region and the position of traps of various types known from seismic data and drilling data, we conclude the presence or absence of an abnormal response of the airspace associated with the alleged hydrocarbon deposits.
Claims (2)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2012104785/28A RU2491580C1 (en) | 2012-02-07 | 2012-02-07 | Method of measuring geophysical characteristics using successive inversion of geoelectric data with additional time filter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2012104785/28A RU2491580C1 (en) | 2012-02-07 | 2012-02-07 | Method of measuring geophysical characteristics using successive inversion of geoelectric data with additional time filter |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| RU2491580C1 true RU2491580C1 (en) | 2013-08-27 |
Family
ID=49163913
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2012104785/28A RU2491580C1 (en) | 2012-02-07 | 2012-02-07 | Method of measuring geophysical characteristics using successive inversion of geoelectric data with additional time filter |
Country Status (1)
| Country | Link |
|---|---|
| RU (1) | RU2491580C1 (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105116453A (en) * | 2015-08-14 | 2015-12-02 | 中国石油天然气股份有限公司 | Transient electromagnetic exploration method and apparatus for tundra natural gas hydrate |
| RU2630852C1 (en) * | 2016-07-15 | 2017-09-13 | Общество С Ограниченной Ответственностью "Сибирская Геофизическая Компания" | Method of forecast of effective capacity of collectors on basis of received polarization parameters and conductivity for selected type of environment |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2235347C1 (en) * | 2003-08-11 | 2004-08-27 | РЫХЛИНСКИЙ Николай Иванович | Method for geoelectrosurveying (variants) |
| US6900639B2 (en) * | 2001-08-07 | 2005-05-31 | Statoil Asa | Electromagnetic method and apparatus for determining the nature of subterranean reservoirs using refracted electromagnetic waves |
| 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 |
| US20100052688A1 (en) * | 2006-02-09 | 2010-03-04 | Electromagnetic Geoservices As | Electromagnetic surveying |
| RU2399931C2 (en) * | 2008-10-23 | 2010-09-20 | Общество с ограниченной ответственностью "Сибирская геофизическая научно-производственная компания" | Method for quantitative dividing effects of electromagnetic induction and induced polarisation |
-
2012
- 2012-02-07 RU RU2012104785/28A patent/RU2491580C1/en active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6900639B2 (en) * | 2001-08-07 | 2005-05-31 | Statoil Asa | Electromagnetic method and apparatus for determining the nature of subterranean reservoirs using refracted electromagnetic waves |
| RU2235347C1 (en) * | 2003-08-11 | 2004-08-27 | РЫХЛИНСКИЙ Николай Иванович | Method for geoelectrosurveying (variants) |
| 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 |
| US20100052688A1 (en) * | 2006-02-09 | 2010-03-04 | Electromagnetic Geoservices As | Electromagnetic surveying |
| RU2399931C2 (en) * | 2008-10-23 | 2010-09-20 | Общество с ограниченной ответственностью "Сибирская геофизическая научно-производственная компания" | Method for quantitative dividing effects of electromagnetic induction and induced polarisation |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105116453A (en) * | 2015-08-14 | 2015-12-02 | 中国石油天然气股份有限公司 | Transient electromagnetic exploration method and apparatus for tundra natural gas hydrate |
| RU2630852C1 (en) * | 2016-07-15 | 2017-09-13 | Общество С Ограниченной Ответственностью "Сибирская Геофизическая Компания" | Method of forecast of effective capacity of collectors on basis of received polarization parameters and conductivity for selected type of environment |
| WO2018013004A1 (en) * | 2016-07-15 | 2018-01-18 | Общество С Ограниченной Ответственностью "Сибирская Геофизическая Компания" | Method for predicting effective collector capacity |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chambers et al. | Testing the ability of surface arrays to monitor microseismic activity | |
| RU2301431C2 (en) | Mode of electrical exploration with using of spacing differentiation of the field of formation on several spreads | |
| Valoroso et al. | Active faults and induced seismicity in the Val d’Agri area (Southern Apennines, Italy) | |
| Baumann-Wilke et al. | P-wave traveltime tomography for a seismic characterization of black shales at shallow depth on Bornholm, Denmark | |
| Witten et al. | Microseismic image-domain velocity inversion: Marcellus Shale case study | |
| Eker et al. | Evaluation of site effect within the tectonic basin in the northern side of Ankara | |
| Dobróka et al. | Interval inversion of borehole data for petrophysical characterization of complex reservoirs | |
| Naseer et al. | Simulating the stratigraphy of meandering channels and point bars of Cretaceous system using spectral decomposition tool, Southwest Pakistan: Implications for petroleum exploration | |
| Mahgoub et al. | Seismic inversion as a predictive tool for porosity and facies delineation in Paleocene fluvial/lacustrine reservoirs, Melut Basin, Sudan | |
| RU2491580C1 (en) | Method of measuring geophysical characteristics using successive inversion of geoelectric data with additional time filter | |
| Li et al. | Characterization of interbedded thin beds using zero-crossing-time amplitude stratal slices | |
| Wang et al. | Seismic attributes for characterization and prediction of carbonate faulted karst reservoirs in the Tarim Basin, China | |
| Kwietniak et al. | Resolution enhancement with relative amplitude preservation for unconventional targets | |
| Maity et al. | Novel fracture zone identifier attribute using geophysical and well log data for unconventional reservoirs | |
| Chan | Subsurface geophysical characterization of the crystalline Canadian Shield in northeastern Alberta: implications for geothermal development | |
| Aminzadeh et al. | Fundamentals of Petroleum Geophysics | |
| Poletto et al. | Seismic interferometry experiment in a shallow cased borehole using a seismic vibrator source | |
| Saad | Sedimentary magnetic anomalies: Part 1—The validity of short-wavelength, low-amplitude SEDMAG anomalies | |
| Alshakhs | Shale play assessment of the Goldwyer formation in the Canning basin using property modelling | |
| Weinzierl et al. | Mass-balance threshold matching of geoelectric and seismic data: a case study from Ketzin | |
| Al-Khazraji | Depth prediction by using various velocity models of Khasib Reservoir in East Baghdad field, Iraq | |
| Paris et al. | Predicting reservoir quality in the Bakken Formation, North Dakota, using petrophysics and 3C seismic data | |
| Shreya et al. | Petrophysical reservoir characterization of Habiganj gas field, Surma Basin, Bangladesh | |
| Zeng | Thin-bed prediction by geomorphology-constrained waveform analysis | |
| Jiang et al. | Quantitative identification of pseudofaults underlying igneous rocks through wave-equation-based velocity modeling: A case study of the Shuntuoguole uplift, Tarim Basin, China |