{ "cells": [ { "cell_type": "markdown", "id": "cfba9926", "metadata": {}, "source": [ "# 65 Displaying Seismic Horizons and Faults\n", "\n", "This notebook illustrates how to display seismic interpretations created in Petrel and exported as Shape Files in Python. The seismic data was acquired in 2021 and was obtained from the Geological Survey of NRW. \n", "\n", "" ] }, { "cell_type": "markdown", "id": "448e8438", "metadata": {}, "source": [ "## Set File Paths and download Tutorial Data\n", "\n", "If you downloaded the latest `GemGIS` version from the Github repository, append the path so that the package can be imported successfully. Otherwise, it is recommended to install `GemGIS` via `pip install gemgis` and import `GemGIS` using `import gemgis as gg`. In addition, the file path to the folder where the data is being stored is set. The tutorial data is downloaded using Pooch (https://www.fatiando.org/pooch/latest/index.html) and stored in the specified folder. Use `pip install pooch` if Pooch is not installed on your system yet." ] }, { "cell_type": "code", "execution_count": 1, "id": "00408871", "metadata": {}, "outputs": [], "source": [ "import warnings\n", "warnings.filterwarnings(\"ignore\")\n", "\n", "import gemgis as gg\n", "import pandas as pd\n", "import geopandas as gpd\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import shapely\n", "from typing import Union\n", "from segysak.segy import segy_loader" ] }, { "cell_type": "code", "execution_count": 2, "id": "9657b36d", "metadata": {}, "outputs": [], "source": [ "file_path ='data/65_displaying_seismic_horizons_and_faults/'\n", "gg.download_gemgis_data.download_tutorial_data(filename=\"65_displaying_seismic_horizons_and_faults.zip\", dirpath=file_path)" ] }, { "cell_type": "markdown", "id": "2120486e", "metadata": {}, "source": [ "## Opening traces of seismic lines\n", "\n", "For this tutorial, data of the Landesseismik Münsterland (see https://www.gd.nrw.de/zip/gd_report_2301s.pdf as reference) is used. " ] }, { "cell_type": "code", "execution_count": 3, "id": "60ab47b6", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", " | Profile | \n", "length | \n", "geometry | \n", "
---|---|---|---|
0 | \n", "GD_NRW_2101 | \n", "25436.00 | \n", "LINESTRING (389796.000 5749009.000, 389810.000... | \n", "
1 | \n", "GD_NRW_2102 | \n", "48120.00 | \n", "LINESTRING (379724.000 5765301.500, 379734.000... | \n", "
<xarray.Dataset>\n", "Dimensions: (cdp: 4806, twt: 1751)\n", "Coordinates:\n", " * cdp (cdp) uint16 10203 10204 10205 10206 ... 15005 15006 15007 15008\n", " * twt (twt) float64 0.0 4.0 8.0 12.0 ... 6.992e+03 6.996e+03 7e+03\n", " cdp_x (cdp) float32 4.215e+05 4.215e+05 4.215e+05 ... 3.797e+05 3.797e+05\n", " cdp_y (cdp) float32 5.743e+06 5.743e+06 5.743e+06 ... 5.765e+06 5.765e+06\n", "Data variables:\n", " data (cdp, twt) float32 0.03581 -0.02099 -0.0793 ... -0.01154 0.04135\n", "Attributes: (12/13)\n", " ns: None\n", " sample_rate: 4.0\n", " text: C 1 Client Geologischer Dienst NRW\\nC 2 Contra...\n", " measurement_system: m\n", " d3_domain: None\n", " epsg: None\n", " ... ...\n", " corner_points_xy: None\n", " source_file: Muenster_2D_PreSTM_Stack_2102_AGC.sgy\n", " srd: None\n", " datatype: None\n", " percentiles: [-3.0195783868438255, -2.780390889497134, -1.2279652...\n", " coord_scalar: -100.0
\n", " | \n", " | data | \n", "cdp_x | \n", "cdp_y | \n", "
---|---|---|---|---|
cdp | \n", "twt | \n", "\n", " | \n", " | \n", " |
10203 | \n", "0.00 | \n", "0.04 | \n", "421512.00 | \n", "5743166.00 | \n", "
4.00 | \n", "-0.02 | \n", "421512.00 | \n", "5743166.00 | \n", "|
8.00 | \n", "-0.08 | \n", "421512.00 | \n", "5743166.00 | \n", "|
12.00 | \n", "-0.12 | \n", "421512.00 | \n", "5743166.00 | \n", "|
16.00 | \n", "-0.12 | \n", "421512.00 | \n", "5743166.00 | \n", "|
... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
15008 | \n", "6984.00 | \n", "-0.11 | \n", "379724.00 | \n", "5765301.50 | \n", "
6988.00 | \n", "-0.03 | \n", "379724.00 | \n", "5765301.50 | \n", "|
6992.00 | \n", "-0.03 | \n", "379724.00 | \n", "5765301.50 | \n", "|
6996.00 | \n", "-0.01 | \n", "379724.00 | \n", "5765301.50 | \n", "|
7000.00 | \n", "0.04 | \n", "379724.00 | \n", "5765301.50 | \n", "
8415306 rows × 3 columns
\n", "<xarray.DataArray 'data' (twt: 1751, cdp: 4806)>\n", "array([[ 0.03581395, 0.03755527, 0.04379281, ..., -0.19376808,\n", " -0.16737616, -0.15387064],\n", " [-0.02099377, -0.02159597, -0.01925519, ..., -0.18929183,\n", " -0.15299934, -0.13157344],\n", " [-0.07930481, -0.08407128, -0.09289104, ..., -0.09974575,\n", " -0.07530349, -0.06080835],\n", " ...,\n", " [ 0.02854965, 0.03610878, 0.03544764, ..., -0.03994363,\n", " -0.0334743 , -0.02788334],\n", " [-0.3293916 , -0.32478482, -0.32344526, ..., -0.02570425,\n", " -0.01836406, -0.01154189],\n", " [-0.3131938 , -0.3141619 , -0.31639624, ..., 0.03036585,\n", " 0.0359443 , 0.04135381]], dtype=float32)\n", "Coordinates:\n", " * cdp (cdp) uint16 10203 10204 10205 10206 ... 15005 15006 15007 15008\n", " * twt (twt) float64 0.0 4.0 8.0 12.0 ... 6.992e+03 6.996e+03 7e+03\n", " cdp_x (cdp) float32 4.215e+05 4.215e+05 4.215e+05 ... 3.797e+05 3.797e+05\n", " cdp_y (cdp) float32 5.743e+06 5.743e+06 5.743e+06 ... 5.765e+06 5.765e+06
\n", " | Type | \n", "Domain | \n", "Droid | \n", "Comment | \n", "ShapeName | \n", "Project | \n", "geometry | \n", "X | \n", "Y | \n", "Z | \n", "Name | \n", "
---|---|---|---|---|---|---|---|---|---|---|---|
0 | \n", "Seismic horizon | \n", "Unknown | \n", "://1d9a2dd1-dd1d-4676-a92e-6057e64d33c2/cb463e... | \n", "NaN | \n", "Base Cretaceous | \n", "C:\\Users\\Nicklas.Ackermann\\Desktop\\Muensterlan... | \n", "POINT Z (421497.819 5743179.962 -708.911) | \n", "421497.82 | \n", "5743179.96 | \n", "558.91 | \n", "2102 | \n", "
1 | \n", "Seismic horizon | \n", "Unknown | \n", "://1d9a2dd1-dd1d-4676-a92e-6057e64d33c2/cb463e... | \n", "NaN | \n", "Base Cretaceous | \n", "C:\\Users\\Nicklas.Ackermann\\Desktop\\Muensterlan... | \n", "POINT Z (421490.638 5743186.924 -708.911) | \n", "421490.64 | \n", "5743186.92 | \n", "558.91 | \n", "2102 | \n", "
2 | \n", "Seismic horizon | \n", "Unknown | \n", "://1d9a2dd1-dd1d-4676-a92e-6057e64d33c2/cb463e... | \n", "NaN | \n", "Base Cretaceous | \n", "C:\\Users\\Nicklas.Ackermann\\Desktop\\Muensterlan... | \n", "POINT Z (421483.457 5743193.886 -710.030) | \n", "421483.46 | \n", "5743193.89 | \n", "560.03 | \n", "2102 | \n", "
3 | \n", "Seismic horizon | \n", "Unknown | \n", "://1d9a2dd1-dd1d-4676-a92e-6057e64d33c2/cb463e... | \n", "NaN | \n", "Base Cretaceous | \n", "C:\\Users\\Nicklas.Ackermann\\Desktop\\Muensterlan... | \n", "POINT Z (421476.276 5743200.848 -711.149) | \n", "421476.28 | \n", "5743200.85 | \n", "561.15 | \n", "2102 | \n", "
4 | \n", "Seismic horizon | \n", "Unknown | \n", "://1d9a2dd1-dd1d-4676-a92e-6057e64d33c2/cb463e... | \n", "NaN | \n", "Base Cretaceous | \n", "C:\\Users\\Nicklas.Ackermann\\Desktop\\Muensterlan... | \n", "POINT Z (421469.095 5743207.810 -711.894) | \n", "421469.10 | \n", "5743207.81 | \n", "561.89 | \n", "2102 | \n", "
\n", " | Name | \n", "geometry | \n", "
---|---|---|
0 | \n", "2102 | \n", "LINESTRING (10204.968 558.911, 10205.975 558.9... | \n", "
\n", " | Name | \n", "geometry | \n", "X | \n", "Y | \n", "distance | \n", "
---|---|---|---|---|---|
0 | \n", "2102 | \n", "POINT (10204.968 558.911) | \n", "10204.97 | \n", "558.91 | \n", "NaN | \n", "
1 | \n", "2102 | \n", "POINT (10205.975 558.911) | \n", "10205.97 | \n", "558.91 | \n", "1.01 | \n", "
2 | \n", "2102 | \n", "POINT (10206.981 560.030) | \n", "10206.98 | \n", "560.03 | \n", "1.50 | \n", "
3 | \n", "2102 | \n", "POINT (10207.985 561.149) | \n", "10207.98 | \n", "561.15 | \n", "1.50 | \n", "
4 | \n", "2102 | \n", "POINT (10208.985 561.894) | \n", "10208.99 | \n", "561.89 | \n", "1.25 | \n", "
... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
4759 | \n", "2102 | \n", "POINT (15000.002 909.671) | \n", "15000.00 | \n", "909.67 | \n", "0.99 | \n", "
4760 | \n", "2102 | \n", "POINT (15000.980 909.075) | \n", "15000.98 | \n", "909.07 | \n", "1.15 | \n", "
4761 | \n", "2102 | \n", "POINT (15002.008 950.317) | \n", "15002.01 | \n", "950.32 | \n", "41.26 | \n", "
4762 | \n", "2102 | \n", "POINT (15003.002 950.094) | \n", "15003.00 | \n", "950.09 | \n", "1.02 | \n", "
4763 | \n", "2102 | \n", "POINT (15004.007 949.935) | \n", "15004.01 | \n", "949.93 | \n", "1.02 | \n", "
4764 rows × 5 columns
\n", "\n", " | Name | \n", "geometry | \n", "
---|---|---|
0 | \n", "Fault interpretation 164 | \n", "LINESTRING (10410.600 309.512, 10417.264 362.367) | \n", "
1 | \n", "Fault interpretation 165 | \n", "LINESTRING (10348.015 117.902, 10359.164 262.483) | \n", "
2 | \n", "Fault interpretation 166 | \n", "LINESTRING (10825.354 111.007, 10827.121 201.4... | \n", "
3 | \n", "Fault interpretation 167 | \n", "LINESTRING (10806.579 103.932, 10812.443 187.4... | \n", "
4 | \n", "Fault interpretation 168 | \n", "LINESTRING (10919.610 20.529, 10929.574 120.41... | \n", "
... | \n", "... | \n", "... | \n", "
140 | \n", "Fault interpretation 52 | \n", "LINESTRING (10378.146 176.226, 10363.864 266.9... | \n", "
141 | \n", "Fault interpretation 53 | \n", "LINESTRING (10438.399 532.834, 10438.399 430.5... | \n", "
142 | \n", "Fault interpretation 56 | \n", "LINESTRING (10661.375 452.519, 10656.012 384.4... | \n", "
143 | \n", "Fault interpretation 57 | \n", "LINESTRING (10783.751 231.979, 10763.524 400.9... | \n", "
144 | \n", "Fault interpretation 58 | \n", "LINESTRING (10737.597 262.077, 10739.699 399.9... | \n", "
145 rows × 2 columns
\n", "