При работе со скважинными данными из las-файлов обычно единственными доступными данными о глубине является измеренная глубина (MD), которая учитывает отметку буровой платформы (поворотного стола).

При интерпретации скважин, когда анализируется только одна скважина, этим фактом можно в определенной степени пренебречь. Однако, когда анализ основан на нескольких скважинах одновременно, необходимо получить истинную вертикальную глубину под водой (TVDSS) каждой скважины, чтобы выполнить их корреляцию.

Чтобы сделать это, в MD нам придется вычесть из каждого значения шага глубины значение поворотного стола.

Это пример данных, над которым будут работать:

Сначала мы загружаем данные методом numpy loadtxt, так как первые две строки txt файла не имеют тип данных float, нам придется их пропустить с помощью параметра skiprows, а также нам нужно будет указать, что разделителями строк являются запятые с параметром-разделителем.

import numpy as np

data = np.loadtxt("example_data.txt", delimiter=",", skiprows=2)
print(data)

Теперь мы выбираем только данные MD, которые соответствуют первому значению в каждом из сгенерированных массивов.

measured_depth_data = list()
for item in data:
 measured_depth_data.append(item[0])

print(measured_depth_data)

Теперь, чтобы вычислить TVDSS, мы присваиваем значение поворотного стола переменной, а затем создаем функцию для вычитания из каждого значения глубины MD значения поворотного стола.

rotary_table_vale = 30
def correct_step_value(depth_step):
  return rotary_table_vale - depth_step
true_vertical_depth_sub_sea_data = [correct_step_value(i) for i in measured_depth_data]

print(true_vertical_depth_sub_sea_data)

И вот оно! Обычная интерпретация в программном обеспечении, которое вычисляет TVDSS, такова: отрицательные значения - это шаги глубины ниже уровня моря, а положительные значения - шаги глубины выше него.