.. _AdvancedExampleDeviatedElasticWellbore: #################################################### Deviated Elastic Wellbore Problem #################################################### ------------------------------------------------------------------ Problem description ------------------------------------------------------------------ This example uses the solid mechanics solver to handle a deviated wellbore problem with open hole completion. This wellbore is subjected to a mud pressure at wellbore wall and undrained condition is assumed (no fluid flow in the rock formation). A segment of the wellbore with isotropic linear elastic deformation is simulated in this case. Far field stresses and gravity effect are excluded. The main goal of this example is to validate the internal wellbore mesh generator and mechanics solver for the case of an inclined wellbore. Analytical results of the radial and hoop stresses, :math:`\sigma_{rr}` and :math:`\sigma_{\theta\theta}`, around the wellbore are expressed as `(Detournay and Cheng, 1988) `__ : .. math:: \sigma_{rr} = p_0 \frac{a^2}{r^2} .. math:: \sigma_{\theta\theta} = -p_0 \frac{a^2}{r^2} where :math:`p_0` is the applied mud pressure at wellbore wall, :math:`a` is the wellbore radius and :math:`r` is the radial coordinate. **Input file** This benchmark example uses no external input files and everything required is contained within two GEOS xml files that are located at: .. code-block:: console inputFiles/wellbore/DeviatedElasticWellbore_base.xml and .. code-block:: console inputFiles/wellbore/DeviatedElasticWellbore_benchmark.xml The corresponding xml file for the integrated test is .. code-block:: console inputFiles/wellbore/DeviatedElasticWellbore_smoke.xml In this example, we would focus our attention on the ``Mesh`` tag. ----------------------------------------------------------- Solid mechanics solver ----------------------------------------------------------- As fluid flow is not considered, only the solid mechanics solver ``SolidMechanicsLagrangianFEM`` is required for solving this wellbore problem. .. literalinclude:: ../../../../../../../inputFiles/wellbore/DeviatedElasticWellbore_benchmark.xml :language: xml :start-after: :end-before: -------------------------------------------------------------------- Deviated wellbore mesh -------------------------------------------------------------------- The internal wellbore mesh generator ``InternalWellbore`` is employed to create the mesh of this wellbore problem. The radius of the wellbore and the size of the surrounding rock formation are defined by a vector ``radius``. In the tangent direction, ``theta`` angle is specified from 0 to 180 degree for a half of the domain regarding its symmetry. Note that the whole domain could be specified with a ``theta`` angle from 0 to 360 degree, if modeling complicated scenarios. The trajectory of the well is defined by ``trajectory``. In this example, the wellbore is inclined in the x-z plane by an angle of 45 degree. The ``autoSpaceRadialElems`` parameter allows optimally increasing the element size from local zone around the wellbore to the far-field zone, which is set to 1 to activate this option. The ``useCartesianOuterBoundary`` transforms the far-field boundary to a squared shape to enforce a Cartesian aligned outer boundary, which eases the loading of the far-field boundary conditions. In this example, this value is set to 0 for the single region along the radial direction. .. literalinclude:: ../../../../../../../inputFiles/wellbore/DeviatedElasticWellbore_benchmark.xml :language: xml :start-after: :end-before: .. figure:: mesh.png :align: center :width: 800 :figclass: align-center -------------------------------------------------------------------- Constitutive law -------------------------------------------------------------------- Isotropic linear elastic constitutive behavior is considered for the rock around the wellbore. Note that the default density is useless in this specific example, as gravity effect is neglected. .. literalinclude:: ../../../../../../../inputFiles/wellbore/DeviatedElasticWellbore_base.xml :language: xml :start-after: :end-before: -------------------------------------------------------------------- Boundary conditions -------------------------------------------------------------------- Far-field boundaries are subjected to roller constraints and in-situ stresses are not considered. The mud pressure on the wellbore wall is defined by ``Traction`` field specification. The nodeset generated by the internal wellbore generator for this face is named as ``rneg``. The traction type is ``normal`` to mimic a pressure that is applied normal to the wellbore wall. The negative sign of the scale is attributed to the negative sign convention for compressive stresses in GEOS. .. literalinclude:: ../../../../../../../inputFiles/wellbore/DeviatedElasticWellbore_base.xml :language: xml :start-after: :end-before: --------------------------------- Results and benchmark --------------------------------- A good agreement between the GEOS results and the corresponding analytical solutions is shown in the figure below: .. plot:: docs/sphinx/advancedExamples/validationStudies/wellboreProblems/deviatedElasticWellbore/deviatedElasticWellboreFigure.py ------------------------------------------------------------------ To go further ------------------------------------------------------------------ **Feedback on this example** This concludes the deviated elastic wellbore example. For any feedback on this example, please submit a `GitHub issue on the project's GitHub page `_.