ProdX
Modeling
Bottomhole Pressure
summary bottomhole pressure is computed through a pressure traverse methodology using standard multiphase flow correlations traverses start with the measured wellhead pressure from daily production docid\ jnogjwoaowkpzzaho3pkp pvt is established through the properties supplied in the wellbore formation docid\ z0gzmed5eninbjj5bbcrp table, and daily rates and pressures come from the daily production docid\ jnogjwoaowkpzzaho3pkp table when establishing rates, downtime correction is applied to scale rates to a full 24 hours of production details surrounding the algorithm, lift specific nuances, and settings are outlined below node depths bottomhole pressure is extracted from the traverse at two node depths for each day of the well's history the gauge node and the reservoir node in the wellperformanceanalysisresults section of the graphql api, there are 7 relevant fields to bottomhole pressure gauge node bottomholepressure gauge node bottomhole pressure, psia bottomholedepth gauge node depth, ft bottomholedepthsource source definition of the gauge node depth reservoir node bottomholepressurereservoir reservoir node bottomhole pressure, psia bottomholedepthreservoir reservoir node depth, ft bottomholedepthsourcereservoir source definition of the reservoir node depth general bottomholepressurecorrelation correlation used for traverse on that day depth sources are defined as follows in the graphql api pump intake end of tubing bottomhole gauge formation top end of casing gauge node the gauge node is a dynamic field the attempts to report a bottomhole pressure at the most common measurement point in the well this is the value plotted on production charts in the prodx application the purpose of this value is to give a comparison point to common measurements (gauge pressure, esp intake, srp intake) to evaluate accuracy and correlation selection the following logic applies based on the lift type and equipment in the well any well type if a downhole gauge is defined in the wellbore downhole equipment docid\ uavqolnoyqyq7fvadhspx table, this will be the gauge node depth while it was installed superceeding any lift specific definitions below if the gauge is defined at a depth below tubing, the system will assume the depth and node definition to be end of tubing if the well does not have tubing and the gauge is defined at a depth below casing, the system will assume the depth and node definition to be end of casing natural flow end of tubing (if tubing was ran), otherwise formation depth gas lift end of tubing esp esp intake depth srp srp intake depth reservoir node the reservoir node will be reported at the top wellbore formation docid\ z0gzmed5eninbjj5bbcrp depth for wells without time varing formations, this will be a fixed depth for the life of the well for wells which have time varying formations, this can, in rare cases, shift if a shallower zone gets introduced later in the life of the well important considerations in a typical configuration, the gauge node is shallower than the reservoir node, so both pressures are extracted from the traverse directly in atypical situations, where the reservoir node is shallower than the gauge node, a hydrostatic correction is applied based on the primary fluid in the well using fluid density as computed at specific depths and temperatures using pvt and pressure traverse methodology (gas or liquid) in all situations where a bottomhole pressure node depth is requested that is shallower than a downhole pump intake, the traverse will always compute pressures down to pump intake and then use a hydrostatic correction to translate that pressure upward natural flow bhp natural flow bhp is modeled using a standard traverse from whp to bottomhole node using the configured multiphase flow correlation tubing flow (or casing flow in the case that no tubing has been ran) is assumed to be the default flow path unless the isannularflow property is set to true for a given day in the daily production docid\ jnogjwoaowkpzzaho3pkp table gas lift bhp gas lift bhp is modeled using a standard traverse from whp to reservoir node using the configured multiphase flow correlation tubing flow/annular injection or annular flow/tubing injection is determined by comparing tubing and casing pressures injection pressure is inferred to be the higher of the two unless isannularflow is explicitly defined in the daily production docid\ jnogjwoaowkpzzaho3pkp table however, unlike natural flow, gas lift requires a depth of injection to be determined to apply appropriate gas rates to a traverse segment this is determined from a formation gas traverse down the annulus (or tubing in the case of annular flow/tubing injection) at the measured injection rate injection pressure at depth is compared to valve opening and closing pressure to determine which valve is open, establishing a depth of injection once depth of injection is established, segments above use formation gas + injection gas rate; segments below only use formation gas rate more details on gas lift methodology can be found in gas lift docid\ sgmcxvyojjtgmnnchfcr esp bhp esp bhp is modeled by modeling a stage by stage head calculation across the pump, when information is available in situations where no pump information is available, intake pressure will be used directly and traversed down to the reservoir node more details of esp methodology can be found in esp docid 5vt8sdbvmypcdxmhgiknp rod pump bhp documentation coming soon model settings settings can be modified at the tenant level (will apply to all new wells added to the system), or at the individual well level it is important to confirm an individual well's settings in the model qualty / model settings or well performance / bottomhole pressure / settings sections of the application before analyzing bottomhole pressure flow method multiphase flow correlation used for hydrostatic, frictional, and accelerational pressure drop calculations across each segment of the well the default correlation is hagedorn and brown petroleum production engineering standard correlations include single phase gas single phase liquid homogeneous liquid homogeneous no slip homeogeneous no slip no friction gray hagedorn and brown petroleum production engineering (default) https //doi org/10 1016/b978 0 7506 8270 1 x5000 2 https //doi org/10 1016/b978 0 7506 8270 1 x5000 2 hagedorn and brown maurer engineering https //www bsee gov/sites/bsee gov/files/tap technical assessment program/300ae pdf https //www bsee gov/sites/bsee gov/files/tap technical assessment program/300ae pdf duns and ros beggs and brill mukherjee and brill baxendall thomas fancher brown poettmann carpenter tuffp ansari tuffp ansari modified orkiszewski custom correlations can be designed using downhole gauge data to train various types of machine learning models to provide enhanced accuracy (<10% median absolute percentage error), provided enough gauge data is available in a given field contact a customer success representative for more information about these correlations xecta ml2 classification algorithm to pick the ideal multiphase flow correlation on any given day based on a variety of static and dynamic features for a given well xecta ml3 uses ml2 to pick the best correlation for the day and then applies a correction factor to the final bottomhole pressure bhp smoothing there are three settings related to bhp smoothing enabled, smoothing method, and window size (days) if enabled, will use a trailing window to smooth the final bottomhole pressures to remove noise introduced by allocated production data by default, this is disabled beggs and brill payne correction applies payne's correction for liquid holdup to the beggs and brill correlation by default this is enabled use gauge data if enabled, this will use gauge data (gaugepressure or espintakepressure from daily production docid\ jnogjwoaowkpzzaho3pkp ) as the source of truth for the output gauge node bottomhole pressure a traverse will still be used to move from the gauge node to the reservoir node bottomhole pressure, which is used for all internal modeling (ddv/reservoir pressure, nodal analysis, pi, pi based forecasting, and artificial lift timing and selection) for natural flow / gas lift wells with a permanent tubing gauge, the reference depth of that gauge (and by proxy the depth of the gauge bhp node) will be read from the wellbore downhole equipment docid\ uavqolnoyqyq7fvadhspx table if there is a gaugepressure defined in the daily production docid\ jnogjwoaowkpzzaho3pkp dataset without a corresponding gauge depth in wellbore downhole equipment docid\ uavqolnoyqyq7fvadhspx table, the gauge depth will be assumed to be end of tubing this setting can also be used in conjunction with the ml2 and ml3 machine learning based correlation methods, but the ml3 correction will not be applied in practice since there is zero implied error at the gauge node it is assumed that any bottomhole gauge data available in conjunction with esp and srp installations is at pump intake, regardless of whether it is imported into the daily production docid\ jnogjwoaowkpzzaho3pkp fields gaugepressure, espintakepressure, or srpintakepressure since the methodology of the use gauge pressure algorithm ignores pump delta pressure calculations, any data in the gaugepressure, espintakepressure, or srpintakepressure fields must be intake pressures espdischargepressure can still be imported into the daily production docid\ jnogjwoaowkpzzaho3pkp table, but a gauge depth must not be defined in the wellbore downhole equipment docid\ uavqolnoyqyq7fvadhspx table in order for this setting to function properly by default, this is disabled max pressure iterations maximum limit of solver iterations to compute pressure gradients for a given segment of the traverse typically multiphase flow correlations converge in less than 5 iterations use caution when increasing this limit as it could cause increases in compute durations for a well in situations where the solver may not ever converge by default, this is set to 100 base segment length to compute a traverse, the wellbore must be discritized into segments the base segment length is effectively the longest segment the algorithm will consider many segments will be automatically subdivided into smaller lengths to account for changing equipment and dimensions by default, this is set to 200 feet friction factor method selects the method used to compute friction factors for pressure gradient calculations by default, this is set to swamee jain available methods include colebrook jain chen fang swamee jain zigrang silvester moody delta pressure tolerance when solving for the total pressure gradient in a given segment, this tolerance is used to define convergence decreasing this tolerance could significantly extent the computation time to achieve convergence, while increasing could sacrifice accuracy by default, this is set to 0 001 psi critical rate method usually considered in gas wells, a liquid loading critical rate is computed for each segment of the traverse this setting defines the desired methodology for critical rates used in basic liquid loading anomalies note xecta also offers advanced, data driven algorithms to detect and predict liquid loading using scada data by default, this is set to coleman available options include turner coleman nagoo xecta (requires scada data) erosional velocity there are two settings related to erosional velocity method and c factor by default, these are set to api rp 14e and 150, respectively bhp scaling if enabled, this allows the user to specific a slope and intercept value to scale the output bottomhole pressure at the gauge and reservoir node using the following formula final bhp = slope traverse bhp + intercept by default, this is disabled