RU2148843C1 - Способ измерения ядерно-магнитных свойств пород в скважине - Google Patents

Способ измерения ядерно-магнитных свойств пород в скважине Download PDF

Info

Publication number
RU2148843C1
RU2148843C1 RU99101414/28A RU99101414A RU2148843C1 RU 2148843 C1 RU2148843 C1 RU 2148843C1 RU 99101414/28 A RU99101414/28 A RU 99101414/28A RU 99101414 A RU99101414 A RU 99101414A RU 2148843 C1 RU2148843 C1 RU 2148843C1
Authority
RU
Russia
Prior art keywords
signal
nuclear
induction
well
rocks
Prior art date
Application number
RU99101414/28A
Other languages
English (en)
Inventor
А.И. Бураков
А.Б. Гутнер
А.А. Попов
Л.Г. Филиппычева
Original Assignee
ВНИИ разведочной геофизики-Рудгеофизика
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by ВНИИ разведочной геофизики-Рудгеофизика filed Critical ВНИИ разведочной геофизики-Рудгеофизика
Priority to RU99101414/28A priority Critical patent/RU2148843C1/ru
Application granted granted Critical
Publication of RU2148843C1 publication Critical patent/RU2148843C1/ru

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

Предлагается помехозащищенный способ комплексных измерений ядерно-магнитных и магнитных свойств пород в процессе каротажа. В каждой точке разреза последовательно возбуждают сигнал ядерной индукции от водородсодержащего вещества и измеряют его частоту. Затем производят поляризацию окружающих пород, подвергают принятый сигнал индукции от ядер водорода свободного флюида прямому преобразованию Фурье, измеряют амплитуду и ширину спектральной линии протонного резонанса, по которым судят об индексе свободного флюида и приведенном времени поперечной релаксации. После этого производят повторное возбуждение сигнала ядерной индукции от водородсодержащего вещества и измеряют его частоту, сравнивают частоты сигналов ядерной индукции и по соотношению частот при первом и третьем измерениях судят об индукции магнитного поля в скважине до и после поляризации, а по их разности судят об остаточной намагниченности пород, а также об истинном времени поперечной релаксации свободного флюида. Все указанные измерения производят в одной и той же точке пространства. В результате преобразования Фурье исходный сигнал переводится из временного в частотное представление, что позволяет применить эффективную цифровую фильтрацию сигнала и измерить амплитуду и ширину спектральной линии протонного резонанса. 1 з.п.ф-лы, 4 ил.

Description

Изобретение относится к скважинным методам изучения ядерно-магнитных и магнитных свойств горных пород.
Известен способ проверки каротажных сигналов, полученных с помощью ядерного магнитного резонанса (пат. США N 3484680, кл. G 01 N 27/00, приор. 25.05.67, публ. 16.12.69), согласно которому в скважинный прибор ядерного магнитного каротажа (ЯМК) помещается протонный образец, который служит для обнаружения ферромагнитных частиц в ближайшем слое породы. Если из окружающей породы не поступает сигнала ядерной индукции, но такой сигнал поступает от протонного образца, то в породе могут присутствовать ферромагнитные частицы, неоднородное магнитное поле которых способно подавить резонансный сигнал свободных флюидов породы. Недостатком способа являются невозможность количественной оценки магнитных свойств пород и одновременного изучения ядерно-магнитных и магнитных свойств, а также отсутствие возможности уменьшения или исключения влияния ферромагнитных включений породы на сигналы ядерной индукции свободных флюидов.
Известен также способ каротажа буровых скважин с использованием эффекта ядерного-магнитного резонанса (пат. США N 3775671, кл. G 01 N 27/78, приор. 3.07.72, публ. 27.11.73), в соответствии с которым для снижения влияния магнитных минералов, находящихся в горных породах, на сигналы ядерной индукции свободных флюидов предлагается перед поляризацией пропускать через катушку зонда ток обратного направления в течение определенного времени. Недостатком способа является невозможность одновременного изучения ядерно-магнитных и магнитных свойств горных пород.
Наиболее близким к предлагаемому является способ, положенный в основу работы устройства для ядерного магнитного каротажа (а.с. СССР N 441541, кл. G 01 V 3/14, приор. 15.12.72, публ. 30.08.74) и заключающийся в использовании однокатушечного зонда для определения ядерно-магнитных свойств пород и тороидального датчика ядерной индукции, заполненного водородсодержащей жидкостью и используемого для эталонировки аппаратуры. Способ включает поляризацию пород постоянным магнитным полем, выключение поля и измерение амплитуды сигнала ядерной индукции ядер водорода в свободном флюиде породы, возбуждение и измерение аналогичным образом ядерной индукции от водородсодержащего образца с тороидальной обмоткой, индуктивно не связанной с окружающими скважину породами. Недостатком способа является невозможность непрерывного последовательного проведения комплексного измерения ядерно-магнитных и магнитных свойств пород.
Кроме того, для всех перечисленных способов измерения ядерно-магнитных свойств пород характерна низкая помехозащищенность в условиях интенсивных помех и шумов, что приводит к понижению точности измерений.
В известных способах измерения параметров сигнала ядерной индукции свободного флюида горных пород используется процедура непосредственного измерения параметров сигнала во временной области. Суть этой процедуры сводится к следующему. После выключения тока поляризации через рамку зонда ядерного магнитного каротажа в последней возникает ЭДС, наведенная свободно прецессирующими магнитными моментами протонов, имеющихся в поровой жидкости пород, пересеченных скважиной. В идеальном случае наведенный в рамке сигнал ядерной индукции имеет вид (Аксельрод С.М., Неретин В.Д. Ядерный магнитный резонанс в нефтегазовой геологии и геофизике - М.: Недра, 1990, с. 15):
Figure 00000002

где Uo - начальная амплитуда сигнала,
T*2 - приведенное время поперечной релаксации,
ω = γB - круговая частота магнитного резонанса,
γ = гиромагнитное отношение протона,
B - индукция магнитного поля в скважине в объеме действия рамки зонда ЯМК.
Начальная амплитуда сигнала Uo пропорциональна числу протонов, создающих сигнал, т.е. пропорциональна индексу свободного флюида (ИСФ) или эффективной пористости пород. Как видно из выражения (1), U = Uo при t = 0, т.е. непосредственно измеряя начальную амплитуду, можно было бы получить информацию об ИСФ. Но в реальном устройстве для ЯМК сразу же после подачи команды на выключение тока поляризации в рамке ЯМК возникают переходные процессы выключения тока, там же могут появиться помехи, связанные с переключением рамки с процесса поляризации на прием слабых сигналов вида (1). Как правило, рамка для зонда ЯМК в режиме приемника указанных сигналов является контуром, настраиваемым в резонанс с частотой ω. В силу довольно высокой добротности контура происходит затягивание процесса установления колебаний в нем. Перечисленные выше явления приводят к тому, что форма полезного сигнала резко отличается от идеального случая (1), что не позволяет судить о величине Uo непосредственным ее измерением при t ≅ 0.
Из временных диаграмм основных процессов, протекающих при проведении ЯМК (фиг. 1), видно, что в момент времени tn начинается отключение тока поляризации в рамке ЯМК (фиг. 1.I), а через интервал времени tm, называемый мертвым временем, катушка работает в режиме приема сигнала свободной прецесси, ЭДС которого представлена на фиг. 1. II. Затягивание сигнала в резонансных цепях приемника приводит к тому, что реальная форма сигнала имеет вид, представленный на фиг. 1. III. Реально для известных зондов ЯМК величина tm составляет от 15 до 30 мс. При временах приведенной поперечной релаксации нефтей порядка 50-80 мс наблюдаемый сигнал может потерять до 50% своей первоначальной величины. ЭДС сигнала свободной прецессии в лучшем случае составляет единицы микровольт, а тепловые шумы рамки и усилителя - около 0,5 мкВ. Кроме того, на входе присутствуют помехи с весьма широкой полосой частот, дополнительно ухудшающие условия приема малых сигналов и измерение их параметров. Поэтому для улучшения помехозащищенности приемного тракта используется узкополосная фильтрация сигнала путем настройки в резонанс рамки ЯМК и элементов, определяющих полосу пропускания приемника. Узкополосная фильтрация сигнала в известных устройствах его обработки производится высокодобротными (Q ≅ 100-180) перестраиваемыми аналоговыми фильтрами, искажающими форму сигнала. Причем требуется перестройка этих фильтров в зависимости от индукции магнитного поля Земли и в скважинном приборе (входной контур), и на поверхности (блок обработки). Искажения формы сигнала тем больше, чем уже полоса пропускания приемного тракта. Таким образом, выделение слабых сигналов с помощью узкополосной фильтрации ведет к снижению точности определения параметров полезного сигнала, особенно в условиях сильных помех и шумов.
Для борьбы с сетевыми помехами и увеличения отношения сигнал/шум используется интегрирование отрезков огибающей сигнала (1) в определенные моменты времени t1 и t2 (фиг. 1. IV). При этом время интегрирования tи выбирается кратным периоду сетевой помехи.
Таким образом, в результате приема, усиления, детектирования и интегрирования определенных участков огибающей сигнала регистрируются средние значения сигнала U1 и U2, приписываемые соответственно средним временам их существования t1 и t2. Для определения искомого значения Uo и времени релаксации Т*2 достаточно решить систему уравнений
Figure 00000003

откуда
Figure 00000004

Решение исходной задачи в виде выражений (3) требует высокой точности определения амплитуд U1 и U2 в заданные интервалы времени t1 и t2.
Как правило, t1 ≥ (60 - 70) мс и t2 ≥ (90 - 100) мс. Для быстро затухающих сигналов (T*2 ≥ 60 мс) регистрируемые в эти времена значения U1 и U2 составляют соответственно 37% и 22% от первоначальной искомой амплитуды Uo, и в условиях интенсивных помех и шумов при отношении сигнал/шум менее 3 погрешность определения параметров по выражению (3) становится более 20%, а результаты наблюдений носят качественный характер. Известный способ борьбы с этим недостатком заключается в попытке проведения измерений сигнала в более ранние времена наблюдения t1 и t2, но они недоступны из-за конечных (и не малых!) значений tм и реакции приемного контура рамки ЯМК.
Помимо рассмотренных технических ограничений возможности определения параметров сигнала по стандартной процедуре существует известное физическое ограничение метода наблюдения сигнала свободной прецессии, связанное специфической зависимостью параметра T*2 от неоднородности индукции ΔB магнитного поля в объеме действия датчика ЯМК (там же, с.9).
Известно, что
Figure 00000005

где T2 - истинное время поперечной релаксации, не зависящее от магнитного поля.
Условие (4) приводит к тому, что в случае высокой неоднородности ΔB поля в объеме датчика время T*2 становится весьма малым, а сигнал затухает так быстро, что становится ненаблюдаемой величиной. Быстрое затухание сигнала происходит потому, что в разных частях исследуемого датчиком пространства прецессия протонов происходит с различной частотой и фазой, что приводит к разрушению ансамбля ядерных спинов, суммарный магнитный момент которых очень быстро стремится к нулю.
В то же время величина T2 определяется характером спин-спинового взаимодействия, а ее значение тесно связано с размером пор, содержащих резонирующее вещество. Поэтому корректное измерение значений T*2 и ΔB позволило бы с помощью (4) судить об истинном времени релаксации T2, свободном от влияния неоднородности поля.
В реальных условиях скважины неоднородность поля может быть связана с техногенным нарушением магнитных минералов (появлением нестабильной компоненты намагниченности стенок скважины), приобретенным в процессе бурения. В процессе бурения скважины наблюдается нарушение магнитной сплошности среды горной породы и баро-термическое воздействие бурового инструмента на магнитные зерна ферромагнитных минералов в стенках скважины. Вследствие этого повышается неоднородность магнитного поля пород, слагающих стенки скважины, что создает существенную помеху для применения метода ЯМК: снижается точность определения ИСФ. Определение вышеуказанной характеристики происходит на фоне магнитовозмущенного состояния среды, что физически некорректно.
Из изложенного видно, что недостатком известного технического решения является невозможность обеспечения требуемой точности комплексных измерений ядерно-магнитных и магнитных свойств флюидов горных пород в условиях повышенного уровня электромагнитных помех и шумов, а также в условиях техногенных нарушений магнитной сплошности горной породы в процессе бурения скважины.
Ставится задача создания помехозащищенного способа комплексного измерения ядерно-магнитных и магнитных свойств пород.
Задача решается за счет того, что в способе измерения ядерно-магнитных свойств пород в скважинах, включающем поляризацию постоянным магнитным полем окружающих скважину пород, последующее выключение этого поля и измерение амплитуды сигнала индукции ядер водорода в свободном флюиде породы, возбуждение и измерение аналогичным образом сигнала ядерной индукции от водородсодержащего вещества, в котором размещена тороидальная обмотка, индуктивно не связанная с окружающими скважину породами, в каждой точке разреза последовательно возбуждают сигнал ядерной индукции от водородсодержащего вещества и измеряют его частоту, затем производят поляризацию окружающих пород, подвергают принятый сигнал индукции от ядер водорода свободного флюида прямому преобразованию Фурье, измеряют амплитуду и ширину спектральной линии протонного резонанса, по которым судят об индексе свободного флюида и приведенном времени поперечной релаксации, после чего производят повторное возбуждение сигнала ядерной индукции от водородсодержащего вещества и измеряют его частоту, сравнивают частоты сигналов ядерной индукции и по соотношению частот при первом и третьем измерениях судят об индукции магнитного поля в скважине до и после поляризации, а по их разности судят об остаточной намагниченности пород, а также об истинном времени поперечной релаксации свободного флюида, причем все указанные измерения производят в одной и той же точке пространства.
Кроме того, дополнительное повышение точности может быть достигнуто за счет того, что перед поляризацией постоянным магнитным полем выполняют обработку пород стенок скважины переменным магнитным полем с убывающей амплитудой.
В результате преобразования Фурье исходный сигнал переводится из временного в частотное представление, что позволяет применить эффективную цифровую фильтрацию сигнала и измерить амплитуду и ширину спектральной линии протонного резонанса, по которым без обратного преобразования во временное представление судят об ИСФ и времени поперечной релаксации, причем для корректного измерения последней используют данные о разности частот при первом и третьем измерениях сигнала.
Магнитная обработка пород стенок скважины переменным магнитным полем убывающей амплитуды, проводимая перед поляризацией постоянным магнитным полем, позволяет провести размагничивание нестабильной компоненты намагниченности стенок скважины, что приводит к уменьшению неоднородности поля в ближайшей к датчику зоне. Уменьшение неоднородности поля увеличивает T*2, снимает дополнительное уширение спектральной линии протонного резонанса и увеличивает отношение сигнал/шум, что приводит к увеличению точности определения ядерно-магнитных и магнитных свойств пород в скважинах.
В результате измерений в одной и той же точке разреза скважины определяются: индукция магнитного поля до, во время и после измерения ядерно-магнитных свойств окружающей среды, степень неоднородности поля, ИСФ и время поперечной релаксации, а также начальная и остаточная намагниченность пород стенок скважины.
Использование при измерениях водородсодержащего образца в сочетании с тороидальным датчиком, ранее использовавшихся для измерения амплитуды сигналов ядерной индукции, позволило создать малогабаритный помехоустойчивый датчик поля, измеряющий модуль вектора магнитной индукции по частоте ядерного магнитного резонанса. Размещение двух таких датчиков выше и ниже датчика ЯМК позволяет в предложенном способе оценивать, с одной стороны, неоднородность поля в объеме рамки датчика ЯМК, а с другой, - индукцию поля до и после измерения ядерно-магнитных свойств среды, что, в свою очередь, позволяет оценить начальную и наведенную намагниченность стенок скважины.
В предлагаемом решении способ искусственного обратного намагничивания пород для снижения влияния магнитных минералов на сигнал ядерной индукции трансформировался в предложение о более эффективной обработке пород стенок скважины переменным магнитным полем с убывающей амплитудой, что повышает точность измерения не только ядерно-магнитных, но и магнитных свойств пород.
Таким образом, использование ряда известных технических решений и принципиально отличного от известных способа обработки сигнала ЯМК позволило осуществить помехозащищенный способ комплексных измерений ядерно-магнитных и магнитных свойств горных пород в скважинах.
На фиг. 1 приведены временные диаграммы процессов, протекающих при общепринятой процедуре проведения ЯМК: I - изменение тока поляризации; II - ЭДС сигнала свободной прецессии; III - форма усиленного сигнала свободной прецессии; IV - огибающая продетектированного сигнала.
На фиг. 2 приведена типичная спектральная линия протонного резонанса.
На фиг. 3 приведен пример временного (I) и частотного (II) представлений сигнала ЯМК с детализацией (III) участка спектра. Временное представление изображено тремя графиками: 1 - исходный сигнал ЯМК; 2 - обратное косинус-преобразование сигнала 1; 3 - обратное полное преобразование сигнала 1 в полосе свертки 10 Гц.
На фиг. 4 приведена структурная схема устройства для измерения ядерно-магнитных и магнитных свойств пород стенок скважины.
Способ реализуется при каротаже на подъеме скважинного прибора. Для этого необходимо произвести следующие операции.
1. В данной точке измерения возбуждают и регистрируют сигнал ядерной индукции, вызванный прецессией протонов жидкости в первом тороидальном датчике. По частоте прецессии определяется индукция магнитного поля в скважине.
2. В ту же точку по стволу скважины перемещают рамку датчика ЯМК, расположенную в скважинном приборе ниже первого тороидального датчика.
3. Производят поляризацию постоянным магнитным полем окружающей скважину среды и возбуждают сигналы ядерной индукции свободного флюида. Принятый сигнал подвергают прямому преобразованию Фурье и по его спектру определяют начальную амплитуду сигнала и приведенное время поперечной релаксации.
4. В ту же точку по стволу скважины перемещают второй тороидальный датчик, расположенный в скважинном приборе ниже рамки датчика ЯМК.
5. Возбуждают и измеряют сигнал ядерной индукции протонов жидкости во втором тороидальном датчике, что позволяет определить индукцию поля в скважине в той же точке пространства, где проводились предыдущие измерения.
6. По результатам измерений по пп. 1,3,5 определяют ИСФ (по начальной амплитуде сигнала), время поперечной релаксации (по времени приведенной поперечной релаксации и измеренной неоднородности магнитного поля с использованием формулы (4)), а также начальную и наведенную остаточную намагниченность пород стенок скважины.
Для того, чтобы измерения по пп. 1,3,5 производились в одной точке, необходимо согласовать расстояния между датчиками, время измерения и скорость каротажа.
Эффективность предлагаемого способа заключается в том, что он построен на использовании конкретного вида Фурье-спектра сигнала свободной прецессии. Нетрудно показать, что для сигнала свободной прецессии, представленного функцией U(t) в соответствии с выражением (1), комплексный Фурье-спектр G(ω) имеет вид
Figure 00000006

где
Figure 00000007

Figure 00000008

Figure 00000009

а энергетический спектр имеет вид
Figure 00000010

График спектра
Figure 00000011
представлен на фиг. 2. Кривая магнитного резонанса имеет колоколообразную форму, вершина которой
Figure 00000012
имеет координату
Figure 00000013
а ширина линии Δf по уровню 0,707 от максимума равна
Figure 00000014
Таким образом, непосредственно по спектру сигнала
Figure 00000015
можно определить U0 и T*2:
Figure 00000016

Магнитный резонанс обладает достаточно узкой спектральной линией с шириной
Figure 00000017
которая при 40 мс ≤ T2 ≤ 160 мс лежит в пределах 2 Гц ≤ Δf ≤ 8 Гц. При средней частоте сигнала fo≅ 2,3 кГц и T * 2 ≅ 100 мс добротность спектральной линии протонного резонанса составляет около 700. Последнее обстоятельство является определяющим для эффективности предлагаемого способа измерения ядерно-магнитных свойств.
Действительно, резонансная частота свободной прецесии флюида ωo= 2πfo достаточно точно определена по данным первого измерения сигнала тороидального датчика, форма сигнала также apriori известна (выражение (6)). Достаточно выделить из всего спектра смеси сигнала с помехами и шумами очень узкую полосу около fo и определить параметры Uo и T*2 по выражениям (7).
Для этого необходимо изменить общепринятую процедуру измерения сигнала путем введения специализированного аналого-цифрового преобразователя (АЦП) и использования прямого преобразования Фурье. Усиленный сигнал свободной прецессии длительностью Т поступает на вход быстродействующего программно-управляемого АЦП, имеющего время дискретизации Δt. На выходе АЦП получим цифровую последовательность Uk(tk), где
Figure 00000018
N - число записанных точек сигнала, k = 0, 1, 2 ... N-1. С помощью процедуры быстрого преобразования Фурье (БПФ) получаем дискретный спектр исходного сигнала, причем шаг по частоте равен 1/Т. Например, при длительности наблюдения сигнала свободной прецессии Т = 500 мс шаг в частотном представлении спектра составляет 2 Гц. Процедуры вычисления и регистрации спектра сигнала свободной прецессии существенно упрощаются, т. к. требуется вычисление и запись не полного спектра принятого аппаратурой сообщения, а лишь узкой полосы около частоты fo ожидаемого протонного резонанса. Вычисленная гребенка амплитуд спектральных составляющих в окрестности резонанса используется для определения Uo и Т*2 по формуле (7). Для более точного определения ширины спектральной линии, особенно для Т*2 > 100 мс, можно воспользоваться интерполяционными формулами и сгустить сеть спектральных составляющих.
Известный характер спектра (6) с квадратичным спадом начального участка резонансной кривой позволяет получить наиболее высокую точность при квадратичной интерполяции
Figure 00000019
Так, для выше рассмотренного примера с шагом БПФ, равным 2 Гц, возможно сгущение сети спектральных составляющих через 0,25 Гц, что позволяет вычислить ширину полосы с точностью 0,005 Гц. Нами было установлено, что с требуемой для эксперимента точностью в 10% достаточно применить БПФ с N=8192, чтобы вычислять параметры Uo и T*2 в реальном масштабе времени (время вычисления tв ≤ 0,8 с) и в условиях сильных помех.
Осуществляемая таким образом узкополосная цифровая фильтрация существенно отличается от аналоговой, используемой в стандартной процедуре известного способа: в нашем случае достаточно прямого перехода из временного представления сигнала в частотное для определения его искомых параметров, что не меняет характера принятого сообщения. При аналоговой фильтрации анализируемый после фильтра сигнал принципиально отличается от исходного из-за искажения формы.
Имеется существенное отличие также в виде регистрируемого первичного материала. В стандартном ЯМК производится регистрация нескольких (например, двух) проинтегрированных участков огибающей сигнала свободной прецессии, по которым в дальнейшем вычисляют Uo и T*2. В предлагаемом "частотном" способе помимо прямых определений Uo и T*2, выполняемых в реальном масштабе времени, имеется возможность записи небольшого интервала спектра, включающего полезный сигнал с помехами и шумами, по которому при анализе данных каротажа можно восстановить временное представление сигнала в любой точке записанного спектра. Такая возможность представляется исключительно полезной в случае, когда в непосредственной близости по частоте от полезного сигнала имеется помеха, например гармоника питающей сети. В этом случае с помощью процедуры обратного БПФ можно провести свертку записанного сообщения на заданной частоте и по характеру временного представления полученного сигнала судить о его происхождении.
Проиллюстрируем вышеизложенное данными фиг. 3. Сигнал ЯМК был смоделирован в виде суммы полезного сигнала и случайной помехи по формуле
Figure 00000020

где Uo = 400 мВ, Т*2 = 80 мс, fo = 2200 Гц, Us = 1600 мВ,
rad - функция генератора случайных чисел от 0 до 1.
График S(t) представлен на фиг. 3.I.1. Выбранные значения параметров сигнала соответствуют примерно ИСФ ≅ 20% при верхнем пределе отношения сигнал/шум, оцениваемом выражением
С/Ш=0,5e-t/100, (9)
где t - текущее время, мс.
Из (9) видно, что отношение сигнал/шум не превышает 0,5, а при "стандартном" для ЯМК времени наблюдения tн = 100 мс С/Ш составляет всего около 0,18.
Несмотря на столь жесткие условия по помехо-сигнальной обстановке, был рассчитан спектр сигнала (8), представленный на фиг. 3.II с шагом по частоте 2 Гц в диапазоне fo ± 64 Гц, а на фиг. 3.III приведен участок спектра полезного сигнала fo ± 8 Гц, детализированный с шагом 0,25 Гц. По данным фиг. 3. III были рассчитаны параметры исходного сигнала: Uop = 440 мВ и Т2p = 76,9 мс. Погрешность определения Uo равна 10%, а Т*2 - 3,9%. Для справки: при стандартной процедуре измерения параметров сигнала корректными считаются измерения при отношении С/Ш > 3. На фиг. 3.I.2 и 3.III.3 представлены соответственно обратные косинус-преобразование и полное преобразование сигнала при полосе свертки 10 Гц, т.е. результаты цифровой фильтрации исходного сигнала в узкой полосе частот.
Сравнение временного (фиг. 3.I.1) и частотного (фиг. 3.II) представлений одного и того же сигнала ЯМК показывает, что в частотном представлении отношение сигнал/шум в несколько раз выше, чем во временном, что позволяет существенно повысить точность определения ИСФ и Т*2 в условиях слабых сигналов и сильных помех.
При этом повышенная точность измерения Т*2 и прецизионное измерение ΔB с помощью двух тороидальных датчиков открывает возможность измерения Т1 - истинного времени поперечной релаксации, не зависящего от магнитного поля (см. выражение (4)).
Помимо очевидной практической эффективности предложенного частотного способа измерения ядерно-магнитных и магнитных свойств пород в скважинах существует и дополнительное техническое преимущество использования этого способа. При частотном представлении сигнала узкополосная фильтрация осуществляется цифровым способом, а при фильтрации во временном представлении - аналоговыми устройствами высокой добротности. В последнем случае на точность определения параметра сильно влияет нестабильность элементов частото-задающих цепей, что является дополнительным источником погрешности метода, который устраняется при использовании предлагаемого способа.
Предлагаемый способ измерения реализован в экспериментальном образце комплексной аппаратуры ядерно-магнитного каротажа и каротажа магнитного поля (ЯМК-КМП) для нефтяных скважин, разработанном в ВИРГ-Рудгеофизика в 1998 г. Структурная схема аппаратуры ЯМК-КМП изображена на фиг. 4. В состав аппаратуры входят скважинный прибор 1, блок измерения 2 и генератор импульсов тока поляризации 3. Скважинный прибор содержит датчики каротажа магнитного поля 4 и 5, датчик ЯМК 6, схему измерения, управления и настройки 7. Скважинный прибор 1 связан с наземной аппаратурой каротажным кабелем 8. В состав измерительного блока 2 входят процессорный модуль 9, плата АЦП 10 и ПЭВМ 11. Датчики 4 и 5 заполнены керосином.
Аппаратура работает следующим образом. Импульсы тока поляризации от генератора 3 по кабелю 8 поступают в скважинный прибор 1, где последовательно во времени подаются на датчики 4, 5 и 6. Возникающий сигнал ядерной индукции с каждого из датчиков поступает на схему 7, где происходит его прием и усиление. Сигнал в аналоговом виде по кабелю 8 поступает в процессорный модуль 9, где происходит его первичная обработка. Так, сигналы датчиков 4 и 5 используются для измерения частоты резонанса f0, по которой определяют модуль вектора геомагнитной индукции, соответственно B1 и В2. Эти же данные о частоте используются для настройки входного контура рамки ЯМК 6 и для выбора частотного диапазона f0±F при БПФ, как, например, для спектра сигнала на фиг. 3. II, где F = 64 Гц и f0 = 2200 Гц. Сигнал датчика ЯМК после его предварительной обработки в схеме 7 и в модуле 9 поступает на АЦП 10, откуда в цифровом виде следует в ПЭВМ 11. В ПЭВМ 11 рассчитывают по участку спектра полезного сигнала его параметры U0 и T*2. По разности показаний датчика поля определяют ΔB и с учетом зависимости (4) вычисляют значение T2.
В результате каротажа на экране монитора ПЭВМ 11 можно одновременно наблюдать в реальном масштабе времени каротажные диаграммы (графики изменения параметров в зависимости от глубины точки измерения) значений U0, T*2, B1 и В2.

Claims (2)

1. Способ измерения ядерно-магнитных свойств пород в скважинах, включающий поляризацию постоянным магнитным полем окружающих скважину пород, последующее выключение этого поля и измерение амплитуды сигнала индукции ядер водорода в свободном флюиде породы, возбуждение и измерение аналогичным образом сигнала ядерной индукции от водородсодержащего вещества, в котором размещена тороидальная обмотка, индуктивно не связанная с окружающими скважину породами, отличающийся тем, что в каждой точке разреза последовательно возбуждают сигнал ядерной индукции от водородсодержащего вещества и измеряют его частоту, затем производят поляризацию окружающих пород, подвергают принятый сигнал индукции от ядер водорода свободного флюида прямому преобразованию Фурье, измеряют амплитуду и ширину спектральной линии протонного резонанса, по которым судят об индексе свободного флюида и приведенном времени поперечной релаксации, после чего производят повторное возбуждение сигнала ядерной индукции от водородсодержащего вещества и измеряют его частоту, сравнивают частоты сигналов ядерной индукции и по соотношению частот при первом и третьем измерениях судят об индукции магнитного поля в скважине до и после поляризации, а по их разности судят об остаточной намагниченности пород, а также об истинном времени поперечной релаксации свободного флюида, причем все указанные измерения производят в одной и той же точке пространства.
2. Способ по п.1, отличающийся тем, что перед поляризацией постоянным магнитным полем выполняют обработку пород стенок скважины переменным магнитным полем с убывающей амплитудой.
RU99101414/28A 1999-01-25 1999-01-25 Способ измерения ядерно-магнитных свойств пород в скважине RU2148843C1 (ru)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU99101414/28A RU2148843C1 (ru) 1999-01-25 1999-01-25 Способ измерения ядерно-магнитных свойств пород в скважине

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU99101414/28A RU2148843C1 (ru) 1999-01-25 1999-01-25 Способ измерения ядерно-магнитных свойств пород в скважине

Publications (1)

Publication Number Publication Date
RU2148843C1 true RU2148843C1 (ru) 2000-05-10

Family

ID=20215056

Family Applications (1)

Application Number Title Priority Date Filing Date
RU99101414/28A RU2148843C1 (ru) 1999-01-25 1999-01-25 Способ измерения ядерно-магнитных свойств пород в скважине

Country Status (1)

Country Link
RU (1) RU2148843C1 (ru)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100373172C (zh) * 2004-10-21 2008-03-05 北京大学 用岩石核磁共振弛豫信号测量地层岩石物性的设备及测量方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100373172C (zh) * 2004-10-21 2008-03-05 北京大学 用岩石核磁共振弛豫信号测量地层岩石物性的设备及测量方法

Similar Documents

Publication Publication Date Title
US5212447A (en) Apparatus and technique for nmr diffusion measurement
US6570382B1 (en) Nuclear magnetic resonance method and logging apparatus
US6242912B1 (en) System and method for lithology-independent gas detection using multifrequency gradient NMR logging
US6069477A (en) Method for improving the accuracy of NMR relaxation distribution analysis with two echo trains
CA2246298C (en) Multiple frequency method for nuclear magnetic resonance longitudinal relaxation measurement and pulsing sequence for power use optimization
EP1397708B1 (en) System and method for determining oil, water and gas saturations for low-field gradient nmr logging tools
US5796252A (en) Nuclear magnetic resonance borehole logging apparatus and method for ascertaining a volume of hydrocarbons independent of a diffusion coefficient
CA2117291C (en) Nuclear magnetic resonance detection of geologic structures
US7622919B2 (en) Nuclear magnetic resonance measurement techniques in non-uniform fields
US6331775B1 (en) Gas zone evaluation by combining dual wait time NMR data with density data
US6023163A (en) Well logging method and apparatus for determining gas and diffusion coefficient using NMR
US20030132749A1 (en) Detecting tool motion effects on nuclear magnetic resonance measurements
Prammer et al. Lithology-lndependent Gas Detection by Gradient-NMR Logging
US6297632B1 (en) Detecting tool motion effects on spin echoes obtained with nuclear magnetic resonance measurements
US20040008027A1 (en) Method and apparatus for detecting diffusion sensitive phases with estimation of residual error in NMR logs
US20040041562A1 (en) Method for magnetic resonance fluid characterization
GB2407167A (en) Determining Properties of flowing fluids
US20030210043A1 (en) Processing NMR data in the presence of ringing
EP0852737B1 (en) Lithology-independent gradient nmr gas detection
US10890685B2 (en) Apparatus and methods for determining properties of hydrogen-containing samples using nuclear magnetic resonance
US6097184A (en) Nuclear magnetic resonance well logging to determine gas-filled porosity and oil-filled porosity of earth formations without a constant static magnetic field gradient
CA2436148C (en) Methods for detecting hydrocarbons from nmr data
RU2148843C1 (ru) Способ измерения ядерно-магнитных свойств пород в скважине
US6850060B2 (en) Method and apparatus for rapid characterization of diffusion
US3238446A (en) Nuclear magnetic resonance technique for measurement of mixtures

Legal Events

Date Code Title Description
MM4A The patent is invalid due to non-payment of fees

Effective date: 20060126