Scribd
Upload
16 views5 pages

Lec1

The document outlines the syllabus for a Reservoir Simulation course at the University of Basra's Petroleum Engineering Department, focusing on numerical methods such as the tridiagonal matrix algorithm and finite difference methods for flow equations. It includes topics like irregular gridding, averaging rock and fluid properties, and multi-phase flow through porous media. Additionally, it assigns homework to create and solve a tridiagonal matrix using Matlab.

Uploaded by

rya336437
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
16 views5 pages

Lec1

The document outlines the syllabus for a Reservoir Simulation course at the University of Basra's Petroleum Engineering Department, focusing on numerical methods such as the tridiagonal matrix algorithm and finite difference methods for flow equations. It includes topics like irregular gridding, averaging rock and fluid properties, and multi-phase flow through porous media. Additionally, it assigns homework to create and solve a tridiagonal matrix using Matlab.

Uploaded by

rya336437
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

University of Basra-College of

Engineering
Petroleum Engineering Department

Subject: Reservoir Simulation Class: 4th Year

Lecturer: Dr. Amani J. Majeed Al-Husseini

Syllabus

 solution of system of difference equations tridiagonal


algorithms, use of irregular Gridding, transmissibility,
 the finite difference from of the flow equation in terms of
transmissibility,
 Averaging of rock and fluid properties, solution of radial from
of the flow equation, two dimensional flow,
 setting up the finite difference from, ordering schemes,
standard row ordering, standard column ordering,
 resulting matrix structure,
 Introduction to multi‐phase flow through porous media.

AJM
Tridiagonal matrix algorithm

In numerical linear algebra, the tridiagonal matrix


algorithm, also known as the Thomas algorithm (named
after Llewellyn Thomas), is a simplified form of Gaussian
elimination that can be used to solve tridiagonal systems of
equations.

or

𝑎 𝑐 𝑥 𝑑
⎡ ⎤
𝑏 𝑎 𝑐 𝑥 𝑑
⎢ ⎥
⎢ 𝑏 𝑎 𝑐 ⎥ 𝑥 𝑑
⎣ 𝑏 𝑎 ⎦ 𝑥 𝑑

For Example
Let A is a Matrix ; X, B are Vectors and
A * X=B, then:

AJM
This is a tri‐diagonal Matrix.

Solution of The tri‐diagonal Matrix


Example:‐

𝟐𝒙𝟏 𝒙𝟐 𝟏
𝟐𝒙𝟏 𝟑𝒙𝟐 𝒙𝟑 𝟐
𝒙𝟐 𝟒𝒙𝟑 𝟐𝒙𝟒 𝟑
𝒙𝟑 𝟑𝒙𝟒 𝟒

For:
A X =B

Since:
𝒂 𝒄𝟏 𝒙𝟏 𝒅𝟏
⎡ 𝟏 ⎤
𝒃 𝒂𝟐 𝒄𝟐 𝒙𝟐 𝒅𝟐
⎢ 𝟐 ⎥
⎢ 𝒃𝟑 𝒂𝟑 𝒄𝟑 ⎥ 𝒙𝟑 𝒅𝟑
⎣ 𝒃𝟒 𝒂𝟒 ⎦ 𝒙𝟒 𝒅𝟒
Let;

∝𝒌 𝒂𝒌 𝒄𝒌 𝟏 𝜷𝒌

𝒃𝒌
𝜷𝒌
𝜶𝒌 𝟏

AJM
𝒂𝟏 𝟐 𝒂𝟐 𝟑 𝒂𝟑 𝟒 𝒂𝟒 𝟑
𝒃𝟐 𝟐 𝒃𝟑 𝟏 𝒃𝟒 𝟏
𝒄𝟏 𝟏 𝒄𝟐 𝟏 𝒄𝟑 𝟐
𝒅𝟏 𝟏 𝒅𝟐 𝟐 𝒅𝟑 𝟑 𝒅𝟒 𝟒

Then:

∝𝟏 𝒂𝟏 𝟐

𝒃𝟐 𝟐
𝜷𝟐 𝟏
𝜶𝟏 𝟐

∝𝟐 𝒂𝟐 𝑪𝟏 𝜷𝟐 𝟑 𝟏∗𝟏 𝟐

𝒃𝟑 𝟏
𝜷𝟑
𝜶𝟐 𝟐
𝟏 𝟕
∝𝟑 𝒂𝟑 𝑪𝟐 𝜷𝟑 𝟒 𝟏∗ .
𝟐 𝟐

𝒃𝟒 𝟏 𝟐
𝜷𝟒
𝜶𝟑 𝟑. 𝟓 𝟕
𝟐 𝟏𝟕
∝𝟒 𝒂𝟒 𝑪𝟑 𝜷𝟒 𝟑 𝟐∗ .
𝟕 𝟕

𝑽𝒌 𝒅𝒌 𝜷𝒌 𝑽𝒌 𝟏

𝑽𝟏 𝒅𝟏 𝜷𝟏 𝑽𝟏 𝟏 𝒅𝟏 𝟏

𝑽𝟐 𝒅𝟐 𝜷𝟐 𝑽𝟏 𝟐 𝟏∗𝟏 𝟏

𝟏 𝟓
𝑽𝟑 𝒅𝟑 𝜷𝟑 𝑽𝟐 𝟑 ∗𝟏
𝟐 𝟐

AJM
𝟐 𝟓 𝟐𝟑
𝑽𝟒 𝒅𝟒 𝜷𝟒 𝑽𝟑 𝟒 ∗
𝟕 𝟐 𝟕

Now we will use the following formula to find the value of (x)

𝑽𝒌 𝒄𝒌 𝒙𝒌 𝟏
𝒙𝒌
∝𝒌

𝑽𝟒 𝑪𝟒 𝒙 𝟒 𝟏 𝟐𝟑
→ 𝒙𝟒
∝𝟒 𝟏
𝑽𝟑 𝑪𝟑 𝒙 𝟒 𝟏
→ 𝒙𝟑
∝𝟑 𝟕
𝑽𝟐 𝑪𝟐 𝒙 𝟑 𝟗
→ 𝒙𝟐
∝𝟐 𝟏𝟕
𝑽𝟏 𝑪𝟏 𝒙 𝟐 𝟒
→ 𝒙𝟏
∝𝟏 𝟏𝟕

HW: By using Matlab, write files for the following;

1- How to create Tridiagonal matrix

2- The steps to solve your own example

AJM

You might also like