<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Drilling and Completion Material Management Improvement by Stochastic Modelling</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Carlos Felipe Acevedo-Gutierrez</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Gabriel Villalobos-Camargo</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Jorge Ivan Romero-Gelvez</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Universidad de Bogotá Jorge Tadeo Lozano</institution>
          ,
          <addr-line>Bogotá</addr-line>
          ,
          <country country="CO">Colombia</country>
        </aff>
      </contrib-group>
      <fpage>66</fpage>
      <lpage>76</lpage>
      <abstract>
        <p>The understanding of logistic phenomena entails a broad field of variables that relate to a wider set of considerations. There is a trade-of between general logistic models and mixed models. The former tends to underfit the real case scenarios, while the latter have to be carefully designed in order to have a clear managerial process. In this paper, it is considered a situation observed in a remote production field in Colombian East flatlands, Rubiales. In this case, the supply of materials had challenges in delivery opportunity, inventory rotation and forecast. Within this paper, we expose the data analysis process needed to construct a stochastic model to resolve those issues. Using the historical information available from the first stages of development in this Field some similarities could be noted and lead to a full evaluation which ended in a full forecast model that reduced material costs in 8,65% (approx.1.75 MUSD in first semester of 2019).</p>
      </abstract>
      <kwd-group>
        <kwd>eol&gt;Logistics</kwd>
        <kwd>Forecast</kwd>
        <kwd>stochastic model</kwd>
        <kwd>material management</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>where  ( ) ∈ ℕ represents the population size at time</p>
      <p>( ) ∶=  1 −  1 2 &gt; 0,  ( ) ∶=  2 +  2 2 &gt; 0, ∀ ∈ (0,  )
and</p>
      <p>∶=  1/ 1 ∈ ℕ,  1 &gt; 0,  2 &gt; 0,  1 &gt; 0,  2 ≥ 0
We propose a variation of this model for deterministic stochastic modelling. See section 3.
1.2. Type well
The considered Field, from now Rubiales, is a remote field in Colombian East Flatlands. The
nearest urban area (7705 inhabitants approx.) is located a hundred kilometers away and does
not count with industrial facilities that could provide the amount of material required by the
drilling and completion development plans. Rubiales has a development plan that covers nine
years with near a thousand wells distributed almost symmetrically in time, additionally, the
geostructure holding the production formation is quite isotropic and regular through the area
which allows considering a unique type well for almost all considerations related to
engineering, design and supply plans. This consideration was accepted during the first two years of
operation until some surplus materials began to grow into the inventory in the second semester
of 2017.</p>
      <p>The main constrain in designing the supply chain of Rubiales was its remoteness.
Furthermore in this paper, we focus primarily on the tubular material required by the drilling and
completion plans, primarily because in ECOPETROL, the operator company, there is diferent
acquisition models for construction and consumable materials. For consumables the model
is consignation and only the used material gets paid, but construction material was delivered
in the warehouse and got paid once confirmed as received, thus, the cost of the construction
materials that would not be used during the campaign will indeed be charged to the project.
Although a well uses several construction materials diferent from tubular material, the
quantity is unitary, it means that a unit of each material is used per well, so in that case, it was not
necessary a forecast to adjust the quantity in order to avoid surplus and therefore related costs.</p>
      <p>As seen in Figure 1 there are 6 types of tubular materials used in the well construction, those
materials vary in length and weight therefore their quantities would be considered in tons
instead of the usual length field units. For the specific quantities refer to Table 1.</p>
    </sec>
    <sec id="sec-2">
      <title>2. Initial material quantity definition</title>
      <p>We start by discussing the schedule of operations documented in the well construction
chronograms, as they are the main guides for the drilling and completion operations. For the
development plan in Rubiales, there was a main assumption that guided the organization of the
activities of the project: each well has the same drilling and completion time no matter which
rig is assigned to. Figure 2 shows a typical well construction chronogram for Rubiales in a
6-month period. It could be noted that there are several rigs working at the same time and that
are short time periods between wells called mobilizations that correspond to the time that the
rig uses to move from well to well and from well location to well location, this particularity
will cause that the full analysis could not be a time-dependent one, this is resolved by using a
stochastic approximation based in the same assumption.</p>
      <p>Considering the former plan documents (Type well design and Well construction
chronogram) it is possible to make a crude forecast by simply operating quantity of materials by the
number of wells (seventy-one) in the considered time period (Table 2). This approximation
could lead to a good estimate of the total of required material, but it does not consider any
of the conditions required by the supply chain. At this point we need to evaluate
commercial conditions for the acquisition of these materials, as per usual any stored material must
be checked and maintained. More importantly, a warranty period covers the first six months
after warehouse material reception, and once this period is over any of the possible
fabrication issues could not be resolved through contractual means. Additionally, the material could
sufer a degradation caused by the storage conditions. Therefore, a high inventory rotation is
preferred, thus the possibility of receiving the whole quantity exposed in table 2 at once is not
acceptable, the deliveries must be thoroughly planned.</p>
      <p>It is possible to construct a planned consumption curve and note stabilized slopes that could
be understood as consumption behavior (consumption rate), two ways of presenting this
information could be used, time-dependent or time-independent, both have important applications
for the model (Figure 4). On the time-dependent curve, we could observe a kind of regular
slope which becomes clearer on the time-independent curve (see section 3.1) proving that the
amount of material considered is constant for each well since the slope shows no variation.
This behavior is, of course, an ideal one, the real behavior difers significantly and that’s why
a further revision of forecasting is necessary.</p>
      <p>In Figure 5 is shown the real delivery of 9 5/8” tubulars, the delivery does not seem to follow
a regular pattern. Overlapping Figure 3 and Figure 5 it is evident that there is little
correspondence between them (Figure 6). An overlap of Figure 5 with the time-dependent curve shown
in Figure 4 leads to the initial system state graphic (Figure 7). In it we could the horizontal
distance between the real delivery curve and the consumption curve corresponds to the material
storage time before use, in other words, it represents the inventory rotation; the largest the
distance between the curves, the longest the storage period, and equivalently, risks associated
to material miss storage and wrongful handling grows.</p>
    </sec>
    <sec id="sec-3">
      <title>3. Consumption rate and proposed model</title>
      <p>The consumption rate could be defined as the speed at which the material is consumed. Since
the operation time is not constant due to the mobilization between wells or well locations,
it is necessary to define an arbitrary independent variable to observe this behavior, so the
consumption rate was established as ton/well. Furthermore, since wells were deployed
sequentially the well number can be used as an arbitrary time variable, thus allowing to create a
comparison graph between planned consumption and real consumption (Figure 8) that shows
the divergence between planned and real consumption rates. Hereby those rates seem to have
a constant slope, an expected trend for the planned consumption rate, since the quantities
follow the initial estimation. zooming on the real consumption curve (Figure 9) the slope shows
a variation from well to well showing that the real consumption is not constant.</p>
      <p>Since the expected consumption rate is defined from plan documents the plot of planned
consumption slope will be a constant equal to the average consumption per well. If the real
consumption well to well is operated to seek the historic average a behavior emerges showing
an asymptotic tendency to the real consumption average.</p>
      <p>
        Next, based in the theorem in [
        <xref ref-type="bibr" rid="ref11">11</xref>
        ] we propose for our deterministic model:
well .
      </p>
      <p>stands for the consumption average in the well  and   stands for the consumption in the

  = ∑  /
 =1
 stands for the derivative in consumption average in instant (well)  ,
 stands for the
consumption average in the well  ,</p>
      <p>− 1 stands for the consumption average in the well  − 1.</p>
      <p>Applying equations 1 and 2 a plot for average planned consumption and real average
consumption could be constructed (Figure 10) and finally noted the real divergence in the
consumption and therefore the surplus material that had been purchased for the whole project
which leads to a reevaluation in the plan documents to reduce the initial quantities to more
realistic quantities. Figure 11 shows the final state of the system once this model was applied
and shows the reduction in the acquisition, the use of the previous storage material and the
adjusted inventory rotation.</p>
    </sec>
    <sec id="sec-4">
      <title>4. Results</title>
      <p>Once the model was applied it was verified in the next 66 wells. Figure 12 shows the final
consumption behavior versus the initially planned consumption and the stabilization in the
slope for the model reveals a significant divergence between planned and real consumption.
Table 3 summarizes the total diferences in each material category which could lead to a buy
surplus of 414.56 tons of material. In Figure 12, 7” tubulars show abnormal behavior due to
operational troubles occurred in wells 58,59 and 64 where consumption rises above any previous
consideration, but the model absorbs the change and return to a regular value once the normal
operation is reached again.</p>
      <p>For the next wells the new stabilization average will be used plus some contingency value
based on normal assurance company policies. This will have a total potential efect over approx.
750 wells to develop until 2027.</p>
      <p>Figure 13 shows the final system state. Three situations are observed:
• For 9 5/8”, 7” and 4 1/2” slotted materials, the system reaches a stabilization with an
inventory rotation around 3 months which match with the general commercial
considerations.
• For 4 1/2” and 3 1/2” materials it was necessary to suspend any acquisition be-cause the
surplus obtained during the information capture phase and previous campaigns fulfilled
the actual campaign, in this case, the added value was not to continue storing material
with a possible lack of future use.
• For 4 1/2” tubing, a stabilization in inventory rotation was reached but a change in the
acquisition policies required to buy additional material in order to avoid stoking out
during the campaign execution, nonetheless the consumption model reached its stabilized
phase.</p>
      <p>In conclusion, the model application for 66 wells reported a substantial reduction in the
planned requirements of material in about 8.65%, representing a total cost reduction of 1.35
million USD dollars, and a reduction in risks and costs associated to the storing of 414 tons of
surplus material that could be estimated in 0,4 million US Dollars. Considering that about 750
wells are still in process of construction in Rubiales the impact of the model could represent a
huge improvement in material management for Ecopetrol in the next seven years.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>E.</given-names>
            <surname>Ahmadi</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D. T.</given-names>
            <surname>Masel</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A. Y.</given-names>
            <surname>Metcalf</surname>
          </string-name>
          ,
          <string-name>
            <given-names>K.</given-names>
            <surname>Schuller</surname>
          </string-name>
          ,
          <article-title>Inventory management of surgical supplies and sterile instruments in hospitals: a literature review, Health Systems 8 (</article-title>
          <year>2019</year>
          )
          <fpage>134</fpage>
          -
          <lpage>151</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>L.</given-names>
            <surname>Galli</surname>
          </string-name>
          ,
          <string-name>
            <given-names>T.</given-names>
            <surname>Levato</surname>
          </string-name>
          ,
          <string-name>
            <given-names>F.</given-names>
            <surname>Schoen</surname>
          </string-name>
          , L. Tigli,
          <article-title>Prescriptive analytics for inventory management in health care</article-title>
          ,
          <source>Journal of the Operational Research Society</source>
          (
          <year>2020</year>
          )
          <fpage>1</fpage>
          -
          <lpage>14</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>F. I. JOSEPH</given-names>
            ,
            <surname>C. O. OMODERO</surname>
          </string-name>
          , U. C.
          <article-title>Okezie, Inventory control management and revenue generating capabilities of oil and gas drilling firms in nigeria</article-title>
          ,
          <source>Annals of Spiru</source>
          Haret University.
          <source>Economic Series</source>
          <volume>19</volume>
          (
          <year>2019</year>
          )
          <fpage>75</fpage>
          -
          <lpage>95</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>S.</given-names>
            <surname>Rajendran</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A. R.</given-names>
            <surname>Ravindran</surname>
          </string-name>
          ,
          <article-title>Inventory management of platelets along blood supply chain to minimize wastage and shortage</article-title>
          ,
          <source>Computers &amp; Industrial Engineering</source>
          <volume>130</volume>
          (
          <year>2019</year>
          )
          <fpage>714</fpage>
          -
          <lpage>730</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>C.</given-names>
            <surname>Zhang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>W.</given-names>
            <surname>Gao</surname>
          </string-name>
          ,
          <string-name>
            <given-names>T.</given-names>
            <surname>Yang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S.</given-names>
            <surname>Guo</surname>
          </string-name>
          ,
          <article-title>Opportunistic maintenance strategy for wind turbines considering weather conditions and spare parts inventory management</article-title>
          ,
          <source>Renewable Energy</source>
          <volume>133</volume>
          (
          <year>2019</year>
          )
          <fpage>703</fpage>
          -
          <lpage>711</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>O.</given-names>
            <surname>Berman</surname>
          </string-name>
          , E. Kim,
          <article-title>Stochastic models for inventory management at service facilities</article-title>
          ,
          <source>Stochastic Models</source>
          <volume>15</volume>
          (
          <year>1999</year>
          )
          <fpage>695</fpage>
          -
          <lpage>718</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>K. K.</given-names>
            <surname>Yin</surname>
          </string-name>
          ,
          <string-name>
            <given-names>G. G.</given-names>
            <surname>Yin</surname>
          </string-name>
          , H. Liu,
          <article-title>Stochastic modeling for inventory and production planning in the paper industry</article-title>
          ,
          <source>AIChE journal 50</source>
          (
          <year>2004</year>
          )
          <fpage>2877</fpage>
          -
          <lpage>2890</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [8]
          <string-name>
            <given-names>A.</given-names>
            <surname>Singh</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J. P.</given-names>
            <surname>Hespanha</surname>
          </string-name>
          ,
          <article-title>A derivative matching approach to moment closure for the stochastic logistic model</article-title>
          ,
          <source>Bulletin of mathematical biology</source>
          <volume>69</volume>
          (
          <year>2007</year>
          )
          <fpage>1909</fpage>
          -
          <lpage>1925</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          [9]
          <string-name>
            <given-names>J. H.</given-names>
            <surname>Matis</surname>
          </string-name>
          ,
          <string-name>
            <given-names>T. R.</given-names>
            <surname>Kife</surname>
          </string-name>
          ,
          <string-name>
            <given-names>P.</given-names>
            <surname>Parthasarathy</surname>
          </string-name>
          ,
          <article-title>On the cumulants of population size for the stochastic power law logistic model</article-title>
          ,
          <source>Theoretical Population Biology</source>
          <volume>53</volume>
          (
          <year>1998</year>
          )
          <fpage>16</fpage>
          -
          <lpage>29</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [10]
          <string-name>
            <given-names>J. H.</given-names>
            <surname>Matis</surname>
          </string-name>
          ,
          <string-name>
            <given-names>T. R.</given-names>
            <surname>Kife</surname>
          </string-name>
          ,
          <article-title>On interacting bee/mite populations: a stochastic model with analysis using cumulant truncation</article-title>
          ,
          <source>Environmental and Ecological Statistics</source>
          <volume>9</volume>
          (
          <year>2002</year>
          )
          <fpage>237</fpage>
          -
          <lpage>258</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          [11]
          <string-name>
            <given-names>A.</given-names>
            <surname>Singh</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J. P.</given-names>
            <surname>Hespanha</surname>
          </string-name>
          ,
          <article-title>Moment closure techniques for stochastic models in population biology</article-title>
          , in: 2006 American Control Conference, IEEE,
          <year>2006</year>
          , pp.
          <fpage>6</fpage>
          -pp.
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>