RU2787967C1 - Способ измерения азимута горизонтальных компонент чувствительности скважинных сейсмоприемников - Google Patents

Способ измерения азимута горизонтальных компонент чувствительности скважинных сейсмоприемников Download PDF

Info

Publication number
RU2787967C1
RU2787967C1 RU2022105929A RU2022105929A RU2787967C1 RU 2787967 C1 RU2787967 C1 RU 2787967C1 RU 2022105929 A RU2022105929 A RU 2022105929A RU 2022105929 A RU2022105929 A RU 2022105929A RU 2787967 C1 RU2787967 C1 RU 2787967C1
Authority
RU
Russia
Prior art keywords
sensitivity
values
angle
azimuth
ksi
Prior art date
Application number
RU2022105929A
Other languages
English (en)
Inventor
Петр Михайлович Уткин
Андрей Юрьевич Кожевников
Original Assignee
Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" (Госкорпорация "Росатом")
Федеральное государственное унитарное предприятие "Всероссийский научно-исследовательский институт автоматики им. Н.Л. Духова" (ФГУП "ВНИИА")
Filing date
Publication date
Application filed by Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" (Госкорпорация "Росатом"), Федеральное государственное унитарное предприятие "Всероссийский научно-исследовательский институт автоматики им. Н.Л. Духова" (ФГУП "ВНИИА") filed Critical Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" (Госкорпорация "Росатом")
Application granted granted Critical
Publication of RU2787967C1 publication Critical patent/RU2787967C1/ru

Links

Abstract

Изобретение относится к измерительной технике и системам обработки информации и может быть использовано для измерения азимута горизонтальных компонент чувствительности скважинных сейсмоприемников. Способ измерения азимута горизонтальных компонент чувствительности скважинных сейсмоприемников заключается в том, что измерения проводятся методом сравнения их направления с измеренным направлением осей чувствительности контрольного сейсмоприемника, при этом отсутствует необходимость размещения дополнительного оборудования в скважине и определены диапазоны частот, каждый из которых используется при расчете азимута при обеспечении заданной точности и определяется как сумма значений угла между осями чувствительности проверяемого и контрольного сейсмоприемников, угла ориентации оси чувствительности контрольного сейсмоприемника относительно направления на магнитный север и магнитного склонения в месте установки сейсмоприемника. Техническим результатом является измерение азимута горизонтальных компонент чувствительности скважинных сейсмоприемников с приемлемой точностью без применения дополнительных технических средств, размещенных в скважине. 1 табл.

Description

Изобретение относится к измерительной технике и системам обработки информации и может быть использовано для измерения азимута горизонтальных компонент чувствительности скважинных сейсмоприемников.
Известны способы определения азимутов и зенитных углов скважин, основанные на эффекте гироскопа, т.е. способности инертной массы сохранять свою ориентацию при вращении вокруг своей оси. Исаченко В.X. Инклинометрия скважин. - М: Недра, 1987.
Недостатками этих способов являются:
- низкая виброустойчивость и, соответственно, большая погрешность измерений при внешних динамических воздействиях;
- необходимость включения такого устройства определения азимута в конструкцию скважинного сейсмоприемника, что приводит к повышению его стоимости и снижению надежности, при этом устройство используется однократно - при установке сейсмоприемника в скважину.
Известен способ определения азимута путем сравнения направления с направлением осей чувствительности контрольного сейсмоприемника при механическом повороте контрольного сейсмоприемника вокруг вертикальной оси на заданный угол с использованием поворотного стола. Патент РФ №2233459, МПК G01V 1/16, G01V 1/40, 27.04.2004.
Недостатками способа являются:
- необходимость использования дополнительных (кроме контрольного сейсмоприемника) приспособлений;
- высокая погрешность измерений, связанная с однократными измерениями азимута для установленного угла поворотного стола;
- отсутствие описания порядка расчета погрешности измерений и способов достижения нормируемых показателей точности при их наличии.
Известен способ определения азимута путем сравнения направления с направлением осей чувствительности контрольного сейсмоприемника при математическом повороте контрольного сейсмоприемника вокруг вертикальной оси на заданный угол. Guralp ЗТВ. Operator's Guide. Document No. MAN-BHO-0001, Issue G. Guralp Systems Ltd, Nov. 2019. Данное техническое решение принято в качестве прототипа.
Прототип обладает следующими недостатками:
- для расчетов используются данные регистрации во всем рабочем диапазоне частот, что приводит к увеличению погрешности однократных вычислений из-за некоррелированности сигналов скважинного и контрольного сейсмоприемников в области высоких частот;
- отсутствует описание порядка расчета азимута, погрешности измерений азимута и способов достижения нормируемых показателей точности при их наличии;
- расчет реализован в программном продукте, позволяющем работать только с одним уникальным типом исходных данных в соответствии со спецификацией разработчика программного продукта.
Техническим результатом является измерение азимута горизонтальных компонент чувствительности скважинных сейсмоприемников с приемлемой точностью без применения дополнительных технических средств, размещенных в скважине.
Технический результат достигается тем, что в скважине размещают проверяемый сейсмоприемник, на поверхности над ним размещают контрольный сейсмоприемник, измеряют угол ориентации оси чувствительности контрольного сейсмоприемника относительно направления на магнитный север; в течение определенного интервала времени регистрируют массивы наборов мгновенных значений откликов проверяемого и контрольного сейсмоприемников {ХСС}, {YCC}, {ХКС}, {УКС} соответственно, полученные с выходов каналов регистрации горизонтальных осей чувствительности с определенной частотой дискретизации; из рабочих диапазонов частот указанных сейсмоприемников выделяют диапазон частот от 0,1 Гц до 1,6 Гц, который разбивают на два набора октав (Δƒ1, Δƒ3, Δƒ5, Δƒ7) и (Δƒ2, Δƒ4, Δƒ6), при этом четные октавы смещены на ½ октавы относительно нечетных; в каждой октаве Δƒj, где j=1,…,7, из массивов {ХСС}, {ХКС}, {УКС} цифровой фильтрацией полосовым фильтром формируют наборы мгновенных значений
Figure 00000001
Figure 00000002
Figure 00000003
по которым создается ряд наборов значений {ХССi}r, {ХКСi}r, {YКСi}r, где r=1,…,Mj - номер выборки объемом Nj в полных выборках {ХСС}, {ХКС}, {УКС}, i=1,…,Nj - номер мгновенного значения в массивах {ХССi}r, {ХКСi}r, {YКСi}r соответственно; для каждого набора значений {ХССi}r, {ХКСi}r, {YКСi}r из набора углов {ϕk}, формируемого путем последовательного перебора значений угла в диапазоне от 0 до 360° с шагом Δϕ, определяется угол поворота оси координат ϕr, для которого в текущей октаве Δƒj будет максимальным коэффициент взаимной корреляции набора мгновенных значений {XSSi}j и набора {X'KSi}j, полученного по формуле
X'KSi=Xi⋅cosϕk+Yi⋅sinϕk,
где Xi, Yi - мгновенные значения из наборов {ХКСi}r, {YКСi}r;
для каждой октавы Δƒj определяют угол между осями чувствительности проверяемого и контрольного сейсмоприемников ϕj усреднением полученных значений углов; для полученного набора средних значений углов {ϕj} определяется набор средних значений {ϕ}L, L≤7, для которого среднее значение угла между осями чувствительности сейсмоприемников AzCCКC имеет минимальную абсолютную неопределенность, а само значение имеет минимальное смещение относительно средних значений угла {ϕj} в октавах Δƒ1-Δƒ5, определяют азимут AzN горизонтальных компонент чувствительности как сумму значений угла между осями чувствительности проверяемого и контрольного сейсмоприемников AzCCКC, угла ориентации оси чувствительности контрольного сейсмоприемника относительно направления на магнитный север AzКC и магнитного склонения М в месте установки сейсмоприемников по формуле
AzN=AzCCКC+AzКC+M.
При определении азимута оси X проверяемого сейсмоприемника выполняют следующие действия:
- в скважине размещают проверяемый сейсмоприемник;
- на поверхности над ним размещают контрольный сейсмоприемник;
- измеряют угол AzКC ориентации оси чувствительности контрольного сейсмоприемника относительно направления на магнитный север;
- регистрируют массивы наборов мгновенных значений откликов проверяемого и контрольного сейсмоприемников {ХСС}, {YCC}, {ХКС}, {УКС} соответственно, полученные с выходов каналов регистрации горизонтальных осей чувствительности за определенный интервал времени Т с определенной частотой дискретизации Fd;
- диапазон частот от 0,1 до 1,6 Гц разбивают на два набора октав (Δƒ1, Δƒ3, Δƒ5, Δƒ7) и (Δƒ2, Δƒ4, Δƒ6), при этом четные октавы смещены на ½ октавы относительно нечетных в соответствии с таблицей 1:
Figure 00000004
- на каждой октаве Δƒj, где j=1,…,7, из массивов {ХСС}, {ХКС}, {УКС} цифровой фильтрацией полосовым фильтром формируют наборы мгновенных значений
Figure 00000005
Figure 00000006
Figure 00000007
по которым создается ряд наборов значений {ХССi}r, {ХКСi}r, {YКСi}r, где r=1,…,Mj - номер выборки объемом Nj, в полных выборках {ХСС}, {ХКС}, {УКС}, i=1,…,Nj - номер мгновенного значения в массивах {ХССi}r, {ХКСi}r, {YКСi}r соответственно;
- для каждого набора значений {ХССi}r, {ХКСi}r, {YКСi}r из набора углов {ϕk}, формируемого путем последовательного перебора значений угла в диапазоне от 0 до 360° с шагом Δϕ=0,25°, определяют угол поворота оси координат ϕr для которого в текущей октаве Δƒj будет максимальным коэффициент взаимной корреляции набора мгновенных значений {XSSi}j и набора {X'KSi}j полученного по формуле
X'KSi = Xi ⋅ cosϕk + Yi ⋅ sinϕk,
где Xi, Yi - мгновенные значения из наборов {ХКСi}r, {YКСi}r;
- для каждой октавы Δƒj определяют угол между осями чувствительности проверяемого и контрольного сейсмоприемников ϕj усреднением полученных значений углов;
- для полученного набора средних значений {ϕj} определяют набор средних значений {ϕ}L, L≤7, для которого среднее значение угла между осями чувствительности сейсмоприемников AzCCКC имеет минимальную абсолютную неопределенность, а само значение имеет минимальное смещение относительно средних значений угла {ϕj} в октавах Δƒ1-Δƒ5;
- определяют азимут горизонтальных компонент чувствительности как сумму значений угла между осями чувствительности проверяемого и контрольного сейсмоприемников AzCCКC угла ориентации оси чувствительности контрольного сейсмоприемника относительно направления на магнитный север AzКC и магнитного склонения М в месте установки сейсмоприемников:
AzN = AzCCКC + AzКC + М.

Claims (5)

  1. Способ измерения азимута горизонтальных компонент чувствительности скважинных сейсмоприемников, заключающийся в том, что в скважине размещают проверяемый сейсмоприемник, на поверхности над ним размещают контрольный сейсмоприемник, измеряют угол ориентации оси чувствительности контрольного сейсмоприемника относительно направления на магнитный север; в течение определенного интервала времени регистрируют массивы наборов мгновенных значений откликов проверяемого и контрольного сейсмоприемников {ХСС}, {YCC}, {ХКС}, {YКС} соответственно, полученные с выходов каналов регистрации горизонтальных осей чувствительности с определенной частотой дискретизации; из рабочих диапазонов частот указанных сейсмоприемников выделяют диапазон частот от 0,1 Гц до 1,6 Гц, который разбивают на два набора октав (Δƒ1, Δƒ3, Δƒ5, Δƒ7) и (Δƒ2, Δƒ4, Δƒ6), при этом четные октавы смещены на ½ октавы относительно нечетных; в каждой октаве Δƒj, где j=1,…,7, из массивов {ХСС}, {ХКС}, {YКС} цифровой фильтрацией полосовым фильтром формируют наборы мгновенных значений
    Figure 00000008
    Figure 00000009
    Figure 00000010
    по которым создается ряд наборов значений {XCCi}r, {ХКСi}r, {YКСi}r, где r=1,…,Mj - номер выборки объемом Nj - в полных выборках {ХСС}, {ХКС}, {YКС}, i=1,…,Nj - номер мгновенного значения в массивах {XCCi}r, {ХКСi}r, {YКСi}r соответственно; для каждого набора значений {XCCi}r, {ХКСi}r, {YКСi}r из набора углов {ϕk}, формируемого путем последовательного перебора значений угла в диапазоне от 0 до 360° с шагом Δϕ, определяется угол поворота оси координат ϕr, для которого в текущей октаве Δƒj будет максимальным коэффициент взаимной корреляции набора мгновенных значений {XSSi}j и набора {X'KSi}j, полученного по формуле
  2. X'KSi=Xi ⋅ cosϕk + Yi ⋅ sinϕk,
  3. где Xi Yi - мгновенные значения из наборов {ХКСi}r, {YКСi}r;
  4. для каждой октавы Δƒj определяют угол между осями чувствительности проверяемого и контрольного сейсмоприемников ϕj усреднением полученных значений углов; для полученного набора средних значений углов {ϕj} определяется набор средних значений {ϕ}L, L≤7, для которого среднее значение угла между осями чувствительности сейсмоприемников AzCCКC имеет минимальную абсолютную неопределенность, а само значение имеет минимальное смещение относительно средних значений угла {ϕj} в октавах Δƒ1-Δƒ5, определяют азимут AzN горизонтальных компонент чувствительности как сумму значений угла между осями чувствительности проверяемого и контрольного сейсмоприемников AzCCКC, угла ориентации оси чувствительности контрольного сейсмоприемника относительно направления на магнитный север AzКC и магнитного склонения М в месте установки сейсмоприемников по формуле
  5. AzN = AzCCКC + AzКC + М.
RU2022105929A 2022-03-04 Способ измерения азимута горизонтальных компонент чувствительности скважинных сейсмоприемников RU2787967C1 (ru)

Publications (1)

Publication Number Publication Date
RU2787967C1 true RU2787967C1 (ru) 2023-01-13

Family

ID=

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1045187A1 (ru) * 1981-08-06 1983-09-30 Институт геологии и геофизики СО АН СССР Способ азимутальной ориентации сейсмоприемников многоприборного трехкомпонентного зонда в скважинах
SU1427313A1 (ru) * 1986-03-03 1988-09-30 Институт физики Земли им.О.Ю.Шмидта Способ определени ориентации сейсмометра
RU2022302C1 (ru) * 1990-03-27 1994-10-30 Вычислительный центр СО РАН Способ определения ориентации трехкомпонентного сейсмометра
RU2209449C1 (ru) * 2002-04-24 2003-07-27 Томский политехнический университет Способ определения ориентации сейсмоприемников и устройство для его осуществления
RU2233459C2 (ru) * 2002-11-01 2004-07-27 Федеральное государственное унитарное предприятие "Научно-исследовательский институт импульсной техники" Устройство для определения азимутальной ориентации сейсмометра
CN107607990A (zh) * 2017-08-07 2018-01-19 中国石油天然气集团公司 三分量检波器水平分量的方向检测方法和装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1045187A1 (ru) * 1981-08-06 1983-09-30 Институт геологии и геофизики СО АН СССР Способ азимутальной ориентации сейсмоприемников многоприборного трехкомпонентного зонда в скважинах
SU1427313A1 (ru) * 1986-03-03 1988-09-30 Институт физики Земли им.О.Ю.Шмидта Способ определени ориентации сейсмометра
RU2022302C1 (ru) * 1990-03-27 1994-10-30 Вычислительный центр СО РАН Способ определения ориентации трехкомпонентного сейсмометра
RU2209449C1 (ru) * 2002-04-24 2003-07-27 Томский политехнический университет Способ определения ориентации сейсмоприемников и устройство для его осуществления
RU2233459C2 (ru) * 2002-11-01 2004-07-27 Федеральное государственное унитарное предприятие "Научно-исследовательский институт импульсной техники" Устройство для определения азимутальной ориентации сейсмометра
CN107607990A (zh) * 2017-08-07 2018-01-19 中国石油天然气集团公司 三分量检波器水平分量的方向检测方法和装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Guralp ЗТВ. Operator's Guide. Document No. MAN-BHO-0001, Issue G. Guralp Systems Ltd, Nov. 2019, pp.42-49. *

Similar Documents

Publication Publication Date Title
RU2106000C1 (ru) Способ для оценки ламинарной структуры и других характеристик земли
Stachnik et al. Determination of New Zealand ocean bottom seismometer orientation via Rayleigh-wave polarization
Wang et al. Seismic sensor misorientation measurement using P‐wave particle motion: An application to the NECsaids array
RU2007137982A (ru) Калибровка по результатам измерений хх, yy и zz, полученным индукционным прибором
US11487038B2 (en) Operating method of a metal detector capable of measuring target depth
US7386942B2 (en) Method and apparatus for mapping the trajectory in the subsurface of a borehole
WO2014163536A1 (ru) Устройство для диагностики технического состояния металлических трубопроводов
Xu et al. Systemic comparison of seismometer horizontal orientations based on teleseismic earthquakes and ambient‐noise data
JPH0743260B2 (ja) 方位角設定機能を有する測量装置
RU2717566C1 (ru) Способ определения погрешностей инерциального блока чувствительных элементов на двухосном поворотном столе
RU2787967C1 (ru) Способ измерения азимута горизонтальных компонент чувствительности скважинных сейсмоприемников
RU2577806C1 (ru) Способ калибровки акселерометрического трехосевого инклинометра
US3149278A (en) Method of geophysical prospecting by measuring the earth's magnetic time transients simultaneously in two different directions
CN112649889A (zh) 一种六分量地震数据及绝对重力测量仪、测量方法
US4672752A (en) Method of determining the difference in borehole azimuth at successive points
US4614040A (en) Borehole survey system and method of determining the difference in borehole azimuth at successive points
RU2433262C1 (ru) Способ контроля азимутальной направленности скважины с использованием gps (варианты) и поверочная инклинометрическая установка для реализации способа контроля азимутальной направленности скважины с использованием gps
RU2357078C2 (ru) Способ определения траектории скважины
Park et al. Correct Off‐Site Determination of Seismic Sensor Orientation from Combined Analyses of Earthquake and Microseism Records
de Chambrier The microlog continuous dipmeter
Li et al. Methods for amplitude calibration and orientation discrepancy measurement: Comparing co‐located sensors of different types in the southern California seismic network
Tasič Interdependent quality control of collocated seismometer and accelerometer
GB2148012A (en) Induced magnetic field borehole surveying method and probe
Rentsch et al. Migration-based location of seismicity recorded with an array installed in the main hole of the San Andreas Fault Observatory at Depth (SAFOD)
RU2436132C1 (ru) Измерительный комплекс для проведения георазведки