Uncertainty Modelling
Description of methodology related to uncertainty modelling
Last updated
Was this helpful?
Description of methodology related to uncertainty modelling
Last updated
Was this helpful?
In this section we describe the content and functionality available in Oliasoft’s Uncertainty modelling module. We begin with a description of the error modelling framework we use, based on the standard developed by The Industry Steering Committee on Wellbore Survey Accuracy (ISCWSA). Then we discuss minimum distance calculations between wellbores, both 3D closest approach, and the earlier horizontal- and perpendicular- scan methods. Finally, we present separation factor calculations based on framework developed by ISCWSA and associates.
In the following, a simple description is given. The framework is as follows. A single wellbore is under consideration, and a set of sensors are used to measure the position along the wellbore in NEV-coordinates. All sensors are assumed to be statistically independent, i.e. a measurement by one does not affect any of the others.
This makes the model linear as a function of sensors, and one can add/subtract sensors as needed.
Furthermore, the wellbore is assumed to consist of a set of survey legs, and each survey leg consists of a set of survey stations. To make it formal, denote the wellbore by W, the survey legs by and the survey stations by . Then,
where and are the cardinality of and , respectively. Also, given a survey station, , there exists a unique survey leg , such that , i.e. if . The wellbore as a continuous curve is described by a parametrisation between the survey stations, usually using the minimum curvature method.
For every sensor , at a given survey station , in survey leg , there is an associated so-called error vector, . This error vector is given by
where is given by the parametrization used of the wellbore between survey stations, is the measurement vector, is the magnitude of the error source over the survey leg, usually quoted at 1 , and finally is the weighting function for the error source. The error summation terminates at the station of interest, and the vector error at this station is given by
Usually, the wellbore is parametrized by the minimum curvature method between survey stations. However, there is no significant loss of accuracy in using the balanced tangential method [ref Williamson], and this is the preferred parametrisation in uncertainty modelling. From these error vectors, a covariance matrix, , is constructed, effectively the outer product of such. This covariance matrix can furthermore be used to calculate uncertainty in position of the wellbore in any direction.
Sensors come in four flavors, depending on how the associated error/uncertainty propagate down the wellbore. Assume we have survey legs, and the survey station of interest is in leg . A short description of the different types of sensors/errors follows.
These are randomly propagating errors, and the contribution to survey station uncertainty from such a source over leg , is
where is the number of survey stations in leg . The total contribution is
These are systematically propagating errors, and the contribution to survey station uncertainty from such a source over leg , is
where is the number of survey stations in leg . The total contribution is
These are error types systematic among all the stations in the well. Hence, we can construct a total vector error from beginning to station , i.e.
where is the number of survey stations in leg . The total contribution to the uncertainty at station is
Well-by-Well errors are (perfectly) correlated within a well, while Global errors are (perfectly) correlated between all wellbores.
The total position covariance at survey station is the sum over all error sources, i.e.
A. We observe that for Random , Well-by-Well and Global error sources, the contribution to the covariance matrix is independent of survey legs, and only dependent on the survey stations. This is not the case for Systematic error sources, which is seen from equation 6 and 7 and the fact that the contribution is not 'linear' in survey stations.
B. When comparing the positional uncertainty between wellbores, the contribution from global error sources should be excluded. If in addition, the wellbores are in the same well, also the well-by-well sources should be omitted.
C. Details about uploading of error/uncertainty- models to the app can be found here (link to decription of IPM/manually defined error models).
An important quantity in its own right, and also in separation factor calculations is the distance between two wellbores. Precisely, given two wellbores and , and denote as the reference wellbore, and as the offset wellbore. Then, given a survey station , calculate the minimum distance between and . Oliasoft offers three different methods for calculating/approximating this distance
3D Closest Approach
Horizontal Scan
Perpendicular Scan
The methods are implemented under the assumptions that the wellbore trajectories are described by the minimum curvature method between survey stations. The last two are included for completeness, and are approximations to the true minimum distance, found by 3D closest approach.
The 3D closest approach method gives the true minimum distance between a survey station and the wellbore , under the assumption that the wellbore trajectory is described by the minimum curvature method between survey stations. The algorithm is conceptually simple, although computationally relatively heavy
Given a survey station .
Parametrize the offset wellbore trajectory, , using the minimum curvature method between survey stations.
Minimize the distance function .
The horizontal scan method is an approximation to the true minimum distance between a survey station and the wellbore . It is a good approximation between vertical/near vertical well, and is useless between horizontal wells. The algorithm goes as follows
Given a survey station , at a vertical depth .
Find the point(s), on the offset wellbore trajectory at the same vertical depth, (such point(s) does not necessarily exist).
Calculate the distance between and (and minimize if more than one).
The perpendicular scan method is an approximation to the true minimum distance between a survey station and the wellbore . It works good between vertical/horizontal wells, but has limitations between deviated wells. The algorithm goes as follows
Given a survey station .
Calculate the normal vector to the wellbore at , and prolong it to it hits the offset well at (this point does not necessarily exist).
Calculate the distance between and .
The separation factor, calculated for every survey station in a reference well against an offset well, is essentially the ratio between the minimum centre-to-centre distance between the two wellbores and the relative positional uncertainty between the two. Currently, Oliasoft offers the Pedal Curve Method as separation factor model which is described in details here. Below we describe the essentials.
Separation factor calculations involve two wellbores, and , referred to as reference well and offset well respectively. For both the wellbores, NEV-covariance matrices, at every survey station have been calculated based on uncertainty models, not necessarily the same model applied to both the wellbores. In addition, for every survey station in the reference well, , the minimum centre-to-centre distance to the offset well is calculated. Then, the separation factor, at survey station is defined as
where is the minimum distance between the survey station and the offset wellbore, , and is the relative uncertainty at one standard deviation between the survey station and the point on the offset well, giving the minimum distance. The other quantities involved are
and are the open hole radius, at the point of interest, of the reference and offset wellbore respectively.
is a surface margin, which effectively increases the radius of the offset well. It is introduced to accommodate small, unidentified errors. Recommended value is m.
is a dimensional scaling factor that determines the probability of well crossing. Recommended value is .
is introduced to quantify, at one standard deviation, the uncertainty in the projection ahead of the current survey station. Recommended value is m.
Calculate the NEV-covariance matrices for the reference and offset wellbore. Remember to exclude any common global error source, and also common well error sources if the wellbores originate from the same well.
Calculate the minimum distance between the survey stations in the reference well and the offset well.
Project the NEV-covariance matrices along the direction defined by the minimum distance, to get the relative uncertainty.
Assemble the separation factor .
The minimum distance calculation depends on which scan method is used, hence, also the separation factor calculations depend on this. Not only through the minimum distance, , but also thorough the relative uncertainty . Best policy is to use the 3D closest approach which gives the true minimum distance.
[1] ISCWSA. Definition of the ISCWSA error model. ISCWSA, 4.3, September 2017.
[2] H. S. Williamson. Accuracy prediction for directional measurement while drilling. SPE Drilling and Completion, 15(04), 2000.
[3] S.J.Sawaryn et.al. Well collision avoidance - separation rule. SPE Annual Technical Conference and Exhibition, 9-11 October, San Antonio, Texas, USA, 2017.