# Corva Welleng
This project is a fork of the [welleng](https://github.com/jonnymaserati/welleng) package originally developed by Jonathan Corcutt.
---
[](https://github.com/pro-well-plan/pwptemp/blob/master/LICENSE.md)
[](https://badge.fury.io/py/welleng)
[](https://pepy.tech/project/welleng)
[](https://opensource.org/licenses/Apache-2.0)
[](https://github.com/jonnymaserati/welleng/actions)
[welleng] aspires to be a collection of useful tools for Wells/Drilling Engineers, kicking off with a range of well trajectory tools.
- Generate survey listings and interpolation with minimum curvature
- Calculate well bore uncertainty data (utilizing either the [ISCWSA] MWD Rev5 models) - the coded error models are within 0.001% accuracy of the ISCWSA test data.
- Calculate well bore clearance and Separation Factors (SF)
- standard [ISCWSA] method within 0.5% accuracy of the ISCWSA test data.
- new mesh based method using the [Flexible Collision Library].
## Support welleng
welleng is fuelled by copious amounts of coffee, so if you wish to supercharge development please donate generously:
<a href="https://www.buymeacoffee.com/jonnymaserati" target="_blank"><img src="https://cdn.buymeacoffee.com/buttons/v2/arial-yellow.png" alt="Buy Me A Coffee" width="217px" ></a>
## New Features!
- **Maximum Curvature Method:** added an alternative `Survey` method for calculating a well trajectory from survey stations that add a more realistic (in terms of tortuosity) versus the traditional *minimum curvature method*. See [this post](https://jonnymaserati.github.io/2022/06/19/modified-tortuosity-index-survey-frequency.html) for more details.
- **Modified Tortuosity Index:** added a `Survey` method for calculating a *modified tortuosity index* as described [here](https://jonnymaserati.github.io/2022/05/26/a-modified-tortuosity-index.html).
- **Panel Plot:** added the `type='panel'` to the `Survey.figure()` method to return plan and section plots.
- **Torque and Drag:** added a simple `torque_drag` module and an `architecture` module for creating scenarios (well bore and simple strings) - see this [post](https://jonnymaserati.github.io/2022/05/22/an-example-of-welleng-torque-drag.html) for instructions.
- **Vertical Section:** this should have been included a long time ago, but finally a vertical section will be calculated if the `vertical_section_azimuth` parameter is included in the `SurveyHeader` when initiating a `Survey` instance. Otherwise, to return the vertical section for a given azimuth (e.g. 45 degrees), or to set the vertical section azimuth and add the vertical section data to the `Survey` then:
```python
survey.get_vertical_section(45)
survey.set_vertical_section(45)
```
- **Hello version 0.4:** major version update to reflect all of the changes happening in the back end. If you have code that's built on previous versions of [welleng] then please lock that version in your env since likely it will require modifying to run with version 0.4 and higher.
- **Project Ahead:** you can now project a survey from the last station to the bit or project to a target to see how to get back on track:
```terminal
>>> node_bit = survey.project_to_bit(delta_md=9.0)
>>> survey_to_target = survey.project_to_target(node_target, dls_design=3.0, delta_md=9.0)
```
- **Interpolate Survey on TVD Depth:** new `survey` function for interpolating fixed TVD intervals along a [welleng] `Survey` instance, e.g. to interpolate `survey` every 10mTVD and return the interpolated survey as `s_interp_tvd`:
```terminal
>>> s_interp_tvd = survey.interpolate_survey_tvd(step=10)
```
- **OWSG Tool Error Models:** the ISCWSA curated Rev 4 and Rev 5 tool models have been coded up and continue to honor the ISCWSA diagnostic data. The OWSG tool errors are ***experimental*** with the following status:
- **Working**: MWD, SRGM, _Fl, SAG, IFR1, IFR2, EMS
- **Not Currently Working Correctly**: AX, GYRO
The available error models can be listed with the following command:
```terminal
>>> welleng.errors.ERROR_MODELS
```
- **World Magnetic Model Calculator:** calculates magnetic field data from the [World Magnetic Model](http://www.geomag.bgs.ac.uk/research/modelling/WorldMagneticModel.html) if magnetic field strength is not provided with the survey data.
- **Import from Landmark .wbp files:** using the `exchange.wbp` module it's now possible to import .wbp files exported from Landmark's COMPASS or DecisionSpace software.
```python
import welleng as we
wp = we.exchange.wbp.load("demo.wbp") # import file
survey = we.exchange.wbp.wbp_to_survey(wp, step=30) # convert to survey
mesh = we.mesh.WellMesh(survey, method='circle') # convert to mesh
we.visual.plot(mesh.mesh) # plot the mesh
```
- **Export to .wbp files *(experimental)*:** using the `exchange.wbp` module, it's possible to convert a planned survey file into a list of turn points that can be exported to a .wbp file.
```python
import welleng as we
wp = we.exchange.wbp.WellPlan(survey) # convert Survey to WellPlan object
doc = we.exchange.wbp.export(wp) # create a .wbp document
we.exchange.wbp.save_to_file(doc, "demo.wbp") # save the document to file
```
- **Well Path Creation:** the addition of the `connector` module enables drilling well paths to be created simply by providing start and end locations (with some vector data like inclination and azimuth). No need to indicate *how* to connect the points, the module will figure that out itself.
- **Fast visualization of well trajectory meshes:** addition of the `visual` module for quick and simple viewing and QAQC of well meshes.
- **Mesh Based Collision Detection:** the current method for determining the Separation Factor between wells is constrained by the frequency and location of survey stations or necessitates interpolation of survey stations in order to determine if Anti-Collision Rules have been violated. Meshing the well bore interpolates between survey stations and as such is a more reliable method for identifying potential well bore collisions, especially wth more sparse data sets.
- More coming soon!
## Tech
[welleng] uses a number of open source projects to work properly:
* [trimesh] - awesome library for loading and using triangular meshes.
* [Flexible Collision Library] - for fast collision detection.
* [numpy] - the fundamental package for scientific computing with Python.
* [scipy] - a Python-based ecosystem of open-source software for mathematics, science, and engineering.
* [vedo] - a python module for scientific visualization, analysis of 3D objects and point clouds based on VTK.
* [magnetic-field-calculator] - a Python API for the British Geological Survey magnetic field calculator.
## Simple Installation
A default, minimal [welleng] installation requires [numpy] and [scipy] which is sufficient for importing or generating trajectories with error models. Other libraries are optional depending on usage - most of [welleng]'s functionality can be unlocked with the `easy` install tag, but if you wish to use mesh collision functionality, then an advanced install is required using the `all` install tag to get [python-fcl], after first installing the compiled dependencies as described below.
You'll receive some `ImportError` messages and a suggested install tag if you try to use functions for which the required dependencies are missing.
### Default install with minimal dependencies:
```
pip install welleng
```
### Easy install with most of the dependencies and no compiled dependencies:
```
pip install welleng[easy]
```
## Advanced Installation
If you want to use the mesh collision detection method, then the compiled dependencies are required prior to installing all of the [welleng] dependencies.
### Ubuntu
Here's how to get the trickier dependencies manually installed on Ubuntu (further instructions can be found [here](https://github.com/flexible-collision-library/fcl/blob/master/INSTALL)):
```terminal
sudo apt-get update
sudo apt-get install libeigen3-dev libccd-dev octomap-tools
```
On a Mac you should be able to install the above with brew and on a Windows machine you'll probably have to build these libraries following the instruction in the link, but it's not too tricky. Once the above are installed, then it should be a simple:
```terminal
pip install welleng[all]
```
For developers, the repository can be cloned and locally installed in your GitHub directory via your preferred Python env (the `dev` branch is usuall a version or two ahead of the `main`).
```terminal
git clone https://github.com/jonnymaserati/welleng.git
cd welleng
pip install -e .[all]
```
Make sure you include that `.` in the final line (it's not a typo) as this ensures that any changes to your development version are immediately implemented on save.
### Windows
Detailed instructions for installing [welleng] in a Windows OS can be found in this [post](https://jonnymaserati.github.io/2021/05/11/install-welleng-windows.html).
### Colaboratory
Perhaps the simplest way of getting up and running with [welleng] is to with a [colab notebook](https://colab.research.google.com/notebooks/intro.ipynb). The required dependencies can be installed with the following cell:
```python
!apt-get install -y xvfb x11-utils libeigen3-dev libccd-dev octomap-tools
!pip install welleng[all]
```
Unfortunately the visualization doesn't work with colab (or rather I've not been able to embed a VTK object) so some further work is needed to view the results. However, the [welleng] engine can be used to generate data in the notebook. Test it out with the following code:
```python
!pip install plotly jupyter-dash pint
!pip install -U git+https://github.com/Kitware/ipyvtk-simple.git
import welleng as we
import plotly.graph_objects as go
from jupyter_dash import JupyterDash
# create a survey
s = we.survey.Survey(
md=[0., 500., 2000., 5000.],
inc=[0., 0., 30., 90],
azi=[0., 0., 30., 90.]
)
# interpolate survey - generate points every 30 meters
s_interp = s.interpolate_survey(step=30)
# plot the results
fig = go.Figure()
fig.add_trace(
go.Scatter3d(
x=s_interp.x,
y=s_interp.y,
z=s_interp.z,
mode='lines',
line=dict(
color='blue'
),
name='survey_interpolated'
),
)
fig.add_trace(
go.Scatter3d(
x=s.x,
y=s.y,
z=s.z,
mode='markers',
marker=dict(
color='red'
),
name='survey'
)
)
fig.update_scenes(zaxis_autorange="reversed")
fig.show()
```
## Quick Start
Here's an example using `welleng` to construct a couple of simple well trajectories with the `connector` module, creating survey listings for the wells with well bore uncertainty data, using these surveys to create well bore meshes and finally printing the results and plotting the meshes with the closest lines and SF data.
```python
import welleng as we
from tabulate import tabulate
# construct simple well paths
print("Constructing wells...")
connector_reference = we.survey.from_connections(
we.connector.Connector(
pos1=[0., 0., 0.],
inc1=0.,
azi1=0.,
pos2=[-100., 0., 2000.],
inc2=90,
azi2=60,
),
step=50
)
connector_offset = we.survey.from_connections(
we.connector.Connector(
pos1=[0., 0., 0.],
inc1=0.,
azi1=225.,
pos2=[-280., -600., 2000.],
inc2=90.,
azi2=270.,
),
step=50
)
# make survey objects and calculate the uncertainty covariances
print("Making surveys...")
sh_reference = we.survey.SurveyHeader(
name="reference",
azi_reference="grid"
)
survey_reference = we.survey.Survey(
md=connector_reference.md,
inc=connector_reference.inc_deg,
azi=connector_reference.azi_grid_deg,
header=sh_reference,
error_model='ISCWSA MWD Rev4'
)
sh_offset = we.survey.SurveyHeader(
name="offset",
azi_reference="grid"
)
survey_offset = we.survey.Survey(
md=connector_offset.md,
inc=connector_offset.inc_deg,
azi=connector_offset.azi_grid_deg,
start_nev=[100., 200., 0.],
header=sh_offset,
error_model='ISCWSA MWD Rev4'
)
# generate mesh objects of the well paths
print("Generating well meshes...")
mesh_reference = we.mesh.WellMesh(
survey_reference
)
mesh_offset = we.mesh.WellMesh(
survey_offset
)
# determine clearances
print("Setting up clearance models...")
c = we.clearance.Clearance(
survey_reference,
survey_offset
)
print("Calculating ISCWSA clearance...")
clearance_ISCWSA = we.clearance.ISCWSA(c)
print("Calculating mesh clearance...")
clearance_mesh = we.clearance.MeshClearance(c, sigma=2.445)
# tabulate the Separation Factor results and print them
results = [
[md, sf0, sf1]
for md, sf0, sf1
in zip(c.reference.md, clearance_ISCWSA.SF, clearance_mesh.SF)
]
print("RESULTS\n-------")
print(tabulate(results, headers=['md', 'SF_ISCWSA', 'SF_MESH']))
# get closest lines between wells
lines = we.visual.get_lines(clearance_mesh)
# plot the result
we.visual.plot(
[mesh_reference.mesh, mesh_offset.mesh], # list of meshes
names=['reference', 'offset'], # list of names
colors=['red', 'blue'], # list of colors
lines=lines
)
print("Done!")
```
This results in a quick, interactive visualization of the well meshes that's great for QAQC. What's interesting about these results is that the ISCWSA method does not explicitly detect a collision in this scenario wheras the mesh method does.
For more examples, including how to build a well trajectory by joining up a series of sections created with the `welleng.connector` module (see pic below), check out the [examples] and follow the [jonnymaserati] blog.
Well trajectory generated by [build_a_well_from_sections.py]
## Todos
- Add a `Target` class to see what you're aiming for - **in progress**
- Documentation
- Generate a scene of offset wells to enable fast screening of collision risks (e.g. hundreds of wells in seconds)
- WebApp for those that just want answers
- Add a `units` module to handle any units system - **in progress**
It's possible to generate data for visualizing well trajectories with [welleng], as can be seen with the rendered scenes below.
ISCWSA Standard Set of Well Paths
The ISCWSA standard set of well paths for evaluating clearance scenarios have been rendered in [blender] above. See the [examples] for the code used to generate a [volve] scene, extracting the data from the [volve] EDM.xml file.
## License
[Apache 2.0](LICENSE)
Please note the terms of the license. Although this software endeavors to be accurate, it should not be used as is for real wells. If you want a production version or wish to develop this software for a particular application, then please get in touch with [jonnycorcutt], but the intent of this library is to assist development.
[//]: # (These are reference links used in the body of this note and get stripped out when the markdown processor does its job. There is no need to format nicely because it shouldn't be seen. Thanks SO - http://stackoverflow.com/questions/4823468/store-comments-in-markdown-syntax)
[jonnycorcutt]: <mailto:jonnycorcutt@gmail.com>
[welleng]: <https://github.com/jonnymaserati/welleng>
[Flexible Collision Library]: <https://github.com/flexible-collision-library/fcl>
[trimesh]: <https://github.com/mikedh/trimesh>
[python-fcl]: <https://github.com/BerkeleyAutomation/python-fcl>
[vedo]: <https://github.com/marcomusy/vedo>
[numpy]: <https://numpy.org/>
[scipy]: <https://www.scipy.org/>
[examples]: <https://github.com/jonnymaserati/welleng/tree/main/examples>
[blender]: <https://www.blender.org/>
[volve]: <https://www.equinor.com/en/how-and-why/digitalisation-in-our-dna/volve-field-data-village-download.html>
[ISCWSA]: <https://www.iscwsa.net/>
[build_a_well_from_sections.py]: <https://github.com/jonnymaserati/welleng/tree/main/examples/build_a_well_from_sections.py>
[magnetic-field-calculator]: <https://pypi.org/project/magnetic-field-calculator/>
[jonnymaserati]: <https://jonnymaserati.github.io/>
Raw data
{
"_id": null,
"home_page": "https://github.com/corva-ai/welleng",
"name": "corva-welleng",
"maintainer": "",
"docs_url": null,
"requires_python": ">=3.8.*",
"maintainer_email": "",
"keywords": "well,trajectory,wellpath,wellbore,drilling,drill,error,separation,minimum curvature,iscwsa,owsg,well engineering,wells,drilling engineering,directional drilling,mwd,survey,covariance,digitalization,automation,volve,witsml",
"author": "Mo Kamyab",
"author_email": "m.kamyab@corva.ai",
"download_url": "https://files.pythonhosted.org/packages/3c/14/011576296fce1c77f7e1fca62fd012269c072a517ad9a79e13ab5bd5373d/corva-welleng-0.5.4.tar.gz",
"platform": null,
"description": "# Corva Welleng\n\nThis project is a fork of the [welleng](https://github.com/jonnymaserati/welleng) package originally developed by Jonathan Corcutt.\n\n---\n\n[](https://github.com/pro-well-plan/pwptemp/blob/master/LICENSE.md)\n[](https://badge.fury.io/py/welleng)\n[](https://pepy.tech/project/welleng)\n[](https://opensource.org/licenses/Apache-2.0)\n[](https://github.com/jonnymaserati/welleng/actions)\n\n[welleng] aspires to be a collection of useful tools for Wells/Drilling Engineers, kicking off with a range of well trajectory tools.\n\n - Generate survey listings and interpolation with minimum curvature\n - Calculate well bore uncertainty data (utilizing either the [ISCWSA] MWD Rev5 models) - the coded error models are within 0.001% accuracy of the ISCWSA test data.\n - Calculate well bore clearance and Separation Factors (SF)\n - standard [ISCWSA] method within 0.5% accuracy of the ISCWSA test data.\n - new mesh based method using the [Flexible Collision Library].\n\n## Support welleng\nwelleng is fuelled by copious amounts of coffee, so if you wish to supercharge development please donate generously: \n\n<a href=\"https://www.buymeacoffee.com/jonnymaserati\" target=\"_blank\"><img src=\"https://cdn.buymeacoffee.com/buttons/v2/arial-yellow.png\" alt=\"Buy Me A Coffee\" width=\"217px\" ></a>\n\n## New Features!\n - **Maximum Curvature Method:** added an alternative `Survey` method for calculating a well trajectory from survey stations that add a more realistic (in terms of tortuosity) versus the traditional *minimum curvature method*. See [this post](https://jonnymaserati.github.io/2022/06/19/modified-tortuosity-index-survey-frequency.html) for more details.\n - **Modified Tortuosity Index:** added a `Survey` method for calculating a *modified tortuosity index* as described [here](https://jonnymaserati.github.io/2022/05/26/a-modified-tortuosity-index.html).\n - **Panel Plot:** added the `type='panel'` to the `Survey.figure()` method to return plan and section plots.\n - **Torque and Drag:** added a simple `torque_drag` module and an `architecture` module for creating scenarios (well bore and simple strings) - see this [post](https://jonnymaserati.github.io/2022/05/22/an-example-of-welleng-torque-drag.html) for instructions.\n - **Vertical Section:** this should have been included a long time ago, but finally a vertical section will be calculated if the `vertical_section_azimuth` parameter is included in the `SurveyHeader` when initiating a `Survey` instance. Otherwise, to return the vertical section for a given azimuth (e.g. 45 degrees), or to set the vertical section azimuth and add the vertical section data to the `Survey` then:\n ```python\n survey.get_vertical_section(45)\n survey.set_vertical_section(45)\n ```\n - **Hello version 0.4:** major version update to reflect all of the changes happening in the back end. If you have code that's built on previous versions of [welleng] then please lock that version in your env since likely it will require modifying to run with version 0.4 and higher.\n - **Project Ahead:** you can now project a survey from the last station to the bit or project to a target to see how to get back on track:\n ```terminal\n >>> node_bit = survey.project_to_bit(delta_md=9.0)\n >>> survey_to_target = survey.project_to_target(node_target, dls_design=3.0, delta_md=9.0)\n ```\n - **Interpolate Survey on TVD Depth:** new `survey` function for interpolating fixed TVD intervals along a [welleng] `Survey` instance, e.g. to interpolate `survey` every 10mTVD and return the interpolated survey as `s_interp_tvd`:\n ```terminal\n >>> s_interp_tvd = survey.interpolate_survey_tvd(step=10)\n ```\n - **OWSG Tool Error Models:** the ISCWSA curated Rev 4 and Rev 5 tool models have been coded up and continue to honor the ISCWSA diagnostic data. The OWSG tool errors are ***experimental*** with the following status:\n - **Working**: MWD, SRGM, _Fl, SAG, IFR1, IFR2, EMS\n - **Not Currently Working Correctly**: AX, GYRO\n \n The available error models can be listed with the following command:\n ```terminal\n >>> welleng.errors.ERROR_MODELS\n ```\n - **World Magnetic Model Calculator:** calculates magnetic field data from the [World Magnetic Model](http://www.geomag.bgs.ac.uk/research/modelling/WorldMagneticModel.html) if magnetic field strength is not provided with the survey data.\n - **Import from Landmark .wbp files:** using the `exchange.wbp` module it's now possible to import .wbp files exported from Landmark's COMPASS or DecisionSpace software.\n ```python\n import welleng as we\n\n wp = we.exchange.wbp.load(\"demo.wbp\") # import file\n survey = we.exchange.wbp.wbp_to_survey(wp, step=30) # convert to survey\n mesh = we.mesh.WellMesh(survey, method='circle') # convert to mesh\n we.visual.plot(mesh.mesh) # plot the mesh\n ```\n \n - **Export to .wbp files *(experimental)*:** using the `exchange.wbp` module, it's possible to convert a planned survey file into a list of turn points that can be exported to a .wbp file.\n\n ```python\n import welleng as we\n\n wp = we.exchange.wbp.WellPlan(survey) # convert Survey to WellPlan object\n doc = we.exchange.wbp.export(wp) # create a .wbp document\n we.exchange.wbp.save_to_file(doc, \"demo.wbp\") # save the document to file\n ```\n \n - **Well Path Creation:** the addition of the `connector` module enables drilling well paths to be created simply by providing start and end locations (with some vector data like inclination and azimuth). No need to indicate *how* to connect the points, the module will figure that out itself.\n - **Fast visualization of well trajectory meshes:** addition of the `visual` module for quick and simple viewing and QAQC of well meshes.\n - **Mesh Based Collision Detection:** the current method for determining the Separation Factor between wells is constrained by the frequency and location of survey stations or necessitates interpolation of survey stations in order to determine if Anti-Collision Rules have been violated. Meshing the well bore interpolates between survey stations and as such is a more reliable method for identifying potential well bore collisions, especially wth more sparse data sets.\n - More coming soon!\n\n## Tech\n\n[welleng] uses a number of open source projects to work properly:\n\n* [trimesh] - awesome library for loading and using triangular meshes.\n* [Flexible Collision Library] - for fast collision detection.\n* [numpy] - the fundamental package for scientific computing with Python.\n* [scipy] - a Python-based ecosystem of open-source software for mathematics, science, and engineering.\n* [vedo] - a python module for scientific visualization, analysis of 3D objects and point clouds based on VTK.\n* [magnetic-field-calculator] - a Python API for the British Geological Survey magnetic field calculator.\n\n## Simple Installation\n\nA default, minimal [welleng] installation requires [numpy] and [scipy] which is sufficient for importing or generating trajectories with error models. Other libraries are optional depending on usage - most of [welleng]'s functionality can be unlocked with the `easy` install tag, but if you wish to use mesh collision functionality, then an advanced install is required using the `all` install tag to get [python-fcl], after first installing the compiled dependencies as described below.\n\nYou'll receive some `ImportError` messages and a suggested install tag if you try to use functions for which the required dependencies are missing.\n\n### Default install with minimal dependencies:\n```\npip install welleng\n```\n### Easy install with most of the dependencies and no compiled dependencies:\n```\npip install welleng[easy]\n```\n## Advanced Installation\nIf you want to use the mesh collision detection method, then the compiled dependencies are required prior to installing all of the [welleng] dependencies.\n### Ubuntu\nHere's how to get the trickier dependencies manually installed on Ubuntu (further instructions can be found [here](https://github.com/flexible-collision-library/fcl/blob/master/INSTALL)):\n\n```terminal\nsudo apt-get update\nsudo apt-get install libeigen3-dev libccd-dev octomap-tools\n```\nOn a Mac you should be able to install the above with brew and on a Windows machine you'll probably have to build these libraries following the instruction in the link, but it's not too tricky. Once the above are installed, then it should be a simple:\n\n```terminal\npip install welleng[all]\n```\n\nFor developers, the repository can be cloned and locally installed in your GitHub directory via your preferred Python env (the `dev` branch is usuall a version or two ahead of the `main`).\n\n```terminal\ngit clone https://github.com/jonnymaserati/welleng.git\ncd welleng\npip install -e .[all]\n```\n\nMake sure you include that `.` in the final line (it's not a typo) as this ensures that any changes to your development version are immediately implemented on save.\n\n### Windows\nDetailed instructions for installing [welleng] in a Windows OS can be found in this [post](https://jonnymaserati.github.io/2021/05/11/install-welleng-windows.html).\n\n### Colaboratory\nPerhaps the simplest way of getting up and running with [welleng] is to with a [colab notebook](https://colab.research.google.com/notebooks/intro.ipynb). The required dependencies can be installed with the following cell:\n\n```python\n!apt-get install -y xvfb x11-utils libeigen3-dev libccd-dev octomap-tools\n!pip install welleng[all]\n```\nUnfortunately the visualization doesn't work with colab (or rather I've not been able to embed a VTK object) so some further work is needed to view the results. However, the [welleng] engine can be used to generate data in the notebook. Test it out with the following code:\n\n```python\n!pip install plotly jupyter-dash pint\n!pip install -U git+https://github.com/Kitware/ipyvtk-simple.git\n\nimport welleng as we\nimport plotly.graph_objects as go\nfrom jupyter_dash import JupyterDash\n\n\n# create a survey\ns = we.survey.Survey(\n md=[0., 500., 2000., 5000.],\n inc=[0., 0., 30., 90],\n azi=[0., 0., 30., 90.]\n)\n\n# interpolate survey - generate points every 30 meters\ns_interp = s.interpolate_survey(step=30)\n\n# plot the results\nfig = go.Figure()\nfig.add_trace(\n go.Scatter3d(\n x=s_interp.x,\n y=s_interp.y,\n z=s_interp.z,\n mode='lines',\n line=dict(\n color='blue'\n ),\n name='survey_interpolated'\n ),\n)\n\nfig.add_trace(\n go.Scatter3d(\n x=s.x,\n y=s.y,\n z=s.z,\n mode='markers',\n marker=dict(\n color='red'\n ),\n name='survey'\n )\n)\nfig.update_scenes(zaxis_autorange=\"reversed\")\nfig.show()\n```\n\n## Quick Start\n\nHere's an example using `welleng` to construct a couple of simple well trajectories with the `connector` module, creating survey listings for the wells with well bore uncertainty data, using these surveys to create well bore meshes and finally printing the results and plotting the meshes with the closest lines and SF data.\n\n```python\nimport welleng as we\nfrom tabulate import tabulate\n\n# construct simple well paths\nprint(\"Constructing wells...\")\nconnector_reference = we.survey.from_connections(\n we.connector.Connector(\n pos1=[0., 0., 0.],\n inc1=0.,\n azi1=0.,\n pos2=[-100., 0., 2000.],\n inc2=90,\n azi2=60,\n ),\n step=50\n)\n\nconnector_offset = we.survey.from_connections(\n we.connector.Connector(\n pos1=[0., 0., 0.],\n inc1=0.,\n azi1=225.,\n pos2=[-280., -600., 2000.],\n inc2=90.,\n azi2=270.,\n ),\n step=50\n)\n\n# make survey objects and calculate the uncertainty covariances\nprint(\"Making surveys...\")\nsh_reference = we.survey.SurveyHeader(\n name=\"reference\",\n azi_reference=\"grid\"\n)\nsurvey_reference = we.survey.Survey(\n md=connector_reference.md,\n inc=connector_reference.inc_deg,\n azi=connector_reference.azi_grid_deg,\n header=sh_reference,\n error_model='ISCWSA MWD Rev4'\n)\nsh_offset = we.survey.SurveyHeader(\n name=\"offset\",\n azi_reference=\"grid\"\n)\nsurvey_offset = we.survey.Survey(\n md=connector_offset.md,\n inc=connector_offset.inc_deg,\n azi=connector_offset.azi_grid_deg,\n start_nev=[100., 200., 0.],\n header=sh_offset,\n error_model='ISCWSA MWD Rev4'\n)\n\n# generate mesh objects of the well paths\nprint(\"Generating well meshes...\")\nmesh_reference = we.mesh.WellMesh(\n survey_reference\n)\nmesh_offset = we.mesh.WellMesh(\n survey_offset\n)\n\n# determine clearances\nprint(\"Setting up clearance models...\")\nc = we.clearance.Clearance(\n survey_reference,\n survey_offset\n)\n\nprint(\"Calculating ISCWSA clearance...\")\nclearance_ISCWSA = we.clearance.ISCWSA(c)\n\nprint(\"Calculating mesh clearance...\")\nclearance_mesh = we.clearance.MeshClearance(c, sigma=2.445)\n\n# tabulate the Separation Factor results and print them\nresults = [\n [md, sf0, sf1]\n for md, sf0, sf1\n in zip(c.reference.md, clearance_ISCWSA.SF, clearance_mesh.SF)\n]\n\nprint(\"RESULTS\\n-------\")\nprint(tabulate(results, headers=['md', 'SF_ISCWSA', 'SF_MESH']))\n\n# get closest lines between wells\nlines = we.visual.get_lines(clearance_mesh)\n\n# plot the result\nwe.visual.plot(\n [mesh_reference.mesh, mesh_offset.mesh], # list of meshes\n names=['reference', 'offset'], # list of names\n colors=['red', 'blue'], # list of colors\n lines=lines\n)\n\nprint(\"Done!\")\n\n```\n\nThis results in a quick, interactive visualization of the well meshes that's great for QAQC. What's interesting about these results is that the ISCWSA method does not explicitly detect a collision in this scenario wheras the mesh method does.\n\n\n\nFor more examples, including how to build a well trajectory by joining up a series of sections created with the `welleng.connector` module (see pic below), check out the [examples] and follow the [jonnymaserati] blog.\n\n\n\nWell trajectory generated by [build_a_well_from_sections.py]\n\n## Todos\n\n - Add a `Target` class to see what you're aiming for - **in progress**\n - Documentation\n - Generate a scene of offset wells to enable fast screening of collision risks (e.g. hundreds of wells in seconds)\n - WebApp for those that just want answers\n - Add a `units` module to handle any units system - **in progress**\n\nIt's possible to generate data for visualizing well trajectories with [welleng], as can be seen with the rendered scenes below.\n\nISCWSA Standard Set of Well Paths\n\nThe ISCWSA standard set of well paths for evaluating clearance scenarios have been rendered in [blender] above. See the [examples] for the code used to generate a [volve] scene, extracting the data from the [volve] EDM.xml file.\n\n## License\n\n[Apache 2.0](LICENSE)\n\nPlease note the terms of the license. Although this software endeavors to be accurate, it should not be used as is for real wells. If you want a production version or wish to develop this software for a particular application, then please get in touch with [jonnycorcutt], but the intent of this library is to assist development.\n\n[//]: # (These are reference links used in the body of this note and get stripped out when the markdown processor does its job. There is no need to format nicely because it shouldn't be seen. Thanks SO - http://stackoverflow.com/questions/4823468/store-comments-in-markdown-syntax)\n\n [jonnycorcutt]: <mailto:jonnycorcutt@gmail.com>\n [welleng]: <https://github.com/jonnymaserati/welleng>\n [Flexible Collision Library]: <https://github.com/flexible-collision-library/fcl>\n [trimesh]: <https://github.com/mikedh/trimesh>\n [python-fcl]: <https://github.com/BerkeleyAutomation/python-fcl>\n [vedo]: <https://github.com/marcomusy/vedo>\n [numpy]: <https://numpy.org/>\n [scipy]: <https://www.scipy.org/>\n [examples]: <https://github.com/jonnymaserati/welleng/tree/main/examples>\n [blender]: <https://www.blender.org/>\n [volve]: <https://www.equinor.com/en/how-and-why/digitalisation-in-our-dna/volve-field-data-village-download.html>\n [ISCWSA]: <https://www.iscwsa.net/>\n [build_a_well_from_sections.py]: <https://github.com/jonnymaserati/welleng/tree/main/examples/build_a_well_from_sections.py>\n [magnetic-field-calculator]: <https://pypi.org/project/magnetic-field-calculator/>\n [jonnymaserati]: <https://jonnymaserati.github.io/>",
"bugtrack_url": null,
"license": "Apache 2.0",
"summary": "A collection of Well Engineering tools",
"version": "0.5.4",
"split_keywords": [
"well",
"trajectory",
"wellpath",
"wellbore",
"drilling",
"drill",
"error",
"separation",
"minimum curvature",
"iscwsa",
"owsg",
"well engineering",
"wells",
"drilling engineering",
"directional drilling",
"mwd",
"survey",
"covariance",
"digitalization",
"automation",
"volve",
"witsml"
],
"urls": [
{
"comment_text": "",
"digests": {
"blake2b_256": "3c14011576296fce1c77f7e1fca62fd012269c072a517ad9a79e13ab5bd5373d",
"md5": "9d12020f69249f161b32dee4e4f96234",
"sha256": "297180336cfb112a694b01c780e23444a2d469fe44b6f1eabccbb076f519c312"
},
"downloads": -1,
"filename": "corva-welleng-0.5.4.tar.gz",
"has_sig": false,
"md5_digest": "9d12020f69249f161b32dee4e4f96234",
"packagetype": "sdist",
"python_version": "source",
"requires_python": ">=3.8.*",
"size": 122085,
"upload_time": "2023-01-27T19:06:00",
"upload_time_iso_8601": "2023-01-27T19:06:00.390500Z",
"url": "https://files.pythonhosted.org/packages/3c/14/011576296fce1c77f7e1fca62fd012269c072a517ad9a79e13ab5bd5373d/corva-welleng-0.5.4.tar.gz",
"yanked": false,
"yanked_reason": null
}
],
"upload_time": "2023-01-27 19:06:00",
"github": true,
"gitlab": false,
"bitbucket": false,
"github_user": "corva-ai",
"github_project": "welleng",
"travis_ci": false,
"coveralls": false,
"github_actions": true,
"requirements": [],
"lcname": "corva-welleng"
}
JSON
