Skip to article

 

Poro-Elastic-Plastic Well Model (2D)

Stress-strain state in the vicinity of a vertical well of radius Rw drilled to a depth of h is determined. The reservoir is considered to be isotropic and homogeneous. The problem is solved in a cylindrical coordinate system.

Geometry creation

1. Create the first circle with radius 10.

On the command bar, select the module for constructing surface geometry (Mode — Geometry, Entity— Surface, Action — Create).

Select Circle in the list of geometric elements.

Set block sizes:

Click Apply.

2. Create the second circle with radius 1.

On the command bar, select the module for constructing surface geometry (Mode - Geometry, Entity - Surface, Action - Create).

From the list of geometric primitives, select Circle.

Set block sizes:

Click Apply.

3. Subtract the first circle from the second one.

On the command bar, select the module for constructing surface geometry (Mode — Geometry, Entity — Surface, Action — Boolean).

Select Subtract in the list of operations.

Set the following parameters:

Click Apply.

As a result, only one body (Body 1) will remain in the model tree.

4. Leave a quarter of the volume (condition of symmetry).

On the command bar, select the module for constructing surface geometry (Mode — Geometry, Entity — Surface, Action — Webcut).

From the list of available section views select Coordinate Plane.

Set the following parameters:

Click Apply.

Do the same, but in the YZ plane:

Click Apply.

Then delete volumes 4 and 1.

5. On the command bar, select the module for constructing volume geometry (Mode — Geometry, Entity — Volume, Action — Delete).

Set the following parameters:

Click Apply.

Setting the material

1. In the command bar, select the module for specifying material properties (Mode — Material, Entity — Material Management).

Specify the name of the material Material 1. Expand the item Elasticity in the left column and drag the Isotropic Material to the Material Properties column. Set the following parameters:

Expand Plasticity in the left column and drag the Second Drucker-Prager Strength Criterion into the Material Properties column. Set the following parameters:

Expand Geomechanics in the left column and drag Bio Isotropic Model to the Material Properties column. Set the following parameters:

Click Apply.

Close the Materials management window.

Meshing

1. On the command bar, select the mesh module (Mode — Mesh, Entity — Curve, Action — Mesh).

Specify the parameters of mesh refinement:

Click Apply.

2. On the command bar, select the mesh module (Mode — Mesh, Entity — Curve, Action — Mesh).

Specify the parameters of mesh refinement:

Click Apply.

3. On the command bar, select the mesh module (Mode — Mesh, Entity — Curve, Action — Mesh).

Specify the parameters of mesh refinement:

Click Apply.

4. Building the mesh.

On the command bar, select the surface mesh module (Mode — Mesh, Entity — Surface, Action — Mesh).

Specify the degree of mesh refinement:

Click Mesh.

Setting boundary conditions

1. Attach curves 8 and 12 in the direction Y and X respectively.

On the command panel, select the boundary conditions module (Mode - Boundary Conditions, Entity - Displacement, Action - Create).

Set the following parameters:

Click Apply.

Set the following parameters:

Click Apply.

2. Set the pore pressure.

On the command panel, select the boundary conditions module (Mode - Boundary Conditions, Entity – Pore Pressure, Action - Create).

Set the following parameters:

Click Apply.

3. Set the pressure on curves 13 and 14.

On the command panel, select the boundary conditions module (Mode - Boundary Conditions, Entity – Pressure, Action - Create).

Set the following parameters:

Click Apply.

Set the following parameters:

Click Apply.

Setting the block properties

1. Create a block of one material type.

On the command bar, select the block management module (Mode — Blocks, Entity — Block, Action — Add).

Set the following parameters:

Click Apply.

2. Set the block parameters.

On the command bar, select the block management module (Mode - Blocks, Entity - Block, Action - Block properties/parameters).

Set the following parameters:

Click Apply .

Starting calculation

1. Set the analysis type.

On the command bar, select the calculation settings module (Mode — Calculation Settings, Calculation Settings — Static, Static — General).

Select:

Click Apply.

Click Start Calculation.

2. In the window that appears, select the directory in which the result will be saved, and enter the file name.

3. In the case of a successful calculation, the console displays the message: Calculation finished successfully at "date" "time".

Results analysis

1. Open the file with the results. There are three ways to do that.

The FidesysViewer window will appear, in which you can view the calculation results.

2. On the toolbar, select FiltersAlphabeticalCoordinate System Conversions.

In the Properties window that opens, set:

Click Apply.

3. Display the σθθ component of the stress field (cylinder) on the model.

On the toolbar, set the following parameters:

4. On the toolbar, select FiltersAlphabeticalPlot Over Line. In the Properties window that opens, set:

 

Using Console Interface

Geometry creating, meshing, setting boundary conditions and materials can be performed using the console interface. Below is a link to the program code that allows you to perform the steps described above manual, you only need to specify the full path and name.

poroelastoplasticity