Building on from my previous post, I wanted to follow up with an example. But first this.
Permitting in GIS
To create a land parcel I can digitise a satellite image, directly record the boundary on the ground, adapt an existing dataset or work with a scanned paper map. I could create a series of buffers around a point, for example create a perimeter around any well pad. The polygon can have land owners’ details, status of the permit and a date the land parcel was last visited.
In the example below I have 10 source lines, the line interval is 400m and the source interval is every 12.5m. My land parcels are highlighted in purple when no permission is granted for access and green for when permission has been granted. I could click on each land parcel and extract more information.
A spatial query provides a way to quickly communicate information about the permitted or unpermitted areas.
Above are 4 fields that are impacted, but permission has been obtained. The image shows a count of the number of shots in that area and whether permission has been granted. It is also possible to show the fields where permission has not been granted and the number of shots, below. In this example 1407 shots do not have permission. Information like this can act as a proactive way to ensure the crew is not stood waiting for permission. If I could have included receiver lines (running at 90 degrees) to show the impact of land permission has on the spread.
If I used each one of these land parcels as a geofence I could record every time I crossed into a permitted or unpermitted area.
To model this, I used QGIS to plot a seismic shot every 10 seconds. In reality the sweep is almost certainly longer than that. The coordinates are in a random location in Sussex, UK (not oilfield) and are generated to plot south to north every 12.5m. The points are then mapped on the screen, created in memory, in QGIS. Every time a point is plotted a call is made to check whether it is in a geofence.
Here is a video of me talking through the process:
The code is available here:
The functions I used are adapted from these sources:
Stopwatch method adapted from http://codereview.stackexchange.com/questions/26534/is-there-a-better-way-to-count-seconds-in-python
Createpoint method adapted from http://gis.stackexchange.com/questions/86812/how-to-draw-polygons-from-the-python-console
point_in_poly method adapted from http://geospatialpython.com/2011/01/point-in-polygon.html
I used QGIS for this example but it could be easily implemented into ArcGIS or any other GIS software.
Want to learn more about GIS and EO for Oil and Gas? Then my page contains all my blogs, plus case studies and links http://gis.acgeospatial.co.uk/