服務熱線:010-64846440
您現在的位置:產品與服務
產品與服務
products
產品與服務
products
公司地址:北京市海淀區西小口路66號D-3-408室
郵 編:100192
電 話:010-64846440 / 010-64868477
E-mail:jiahewuxian@163.com
化學驅數值模擬軟件CDNS使用說明

化學驅數值模擬軟件

 CDNS(4.0)

使用說明

戴  家  林

2008年1月

目  錄

一、組分設計 1

二、化學離子設計 2

三、CDNS軟件運行流程 3

四、數據流的填寫方法 4

五、關鍵字詳細描述 7

1、TTT 7

2、DATE 7

3、TMAX 7

4、UNIT 8

5、INJECT 8

6、GRIDMD 9

7、GRIDVL 15

8、WELLVL 19

9、PARA 21

10、A-DATA 28

11、S-DATA 30

12、P-DATA 35

13、T-DATA 41

14、SLUG 44

15、DT 46

16、ENDFW 47

17、PERFOR 48

18、TIME 48

19、STOP 49

六、數據流的例子 49

七、結果輸出文件 57

1、數據檢查信息文件 57

2、指標分析信息文件 58

3、曲線繪制數據文件 58

4、等值圖和三維可視圖繪制數據文件 59

5、流線圖繪制數據文件 59


化學驅數值模擬軟件CDNS使用說明


一、組分設計

CDNS能夠模擬堿、表活劑、聚合物、示蹤劑等化學劑加入注入液的驅油過程。軟件設計了下列物質組分:

1、水

2、油

3、表活劑

4、聚合物

5、陰離子

6、二階陽離子

7、鈉

8、氫

9、碳酸根(擬組分)

10、石油酸(擬組分)

11、第1種示蹤劑

12、第2種示蹤劑

13、第3種示蹤劑

根據注入化學劑的情況和原油的性質來確定組分數。

1、如果注入液中沒有任何化學劑,即為純水驅,組分數為2(1-2)。

2、如果注入液中沒有堿,也沒有示蹤劑,組分數為6(1-6)。

3、如果注入液中沒有堿,有示蹤劑,組分數為6加示蹤劑的個數。

4、如果注入堿,原油中沒有石油酸,組分數為9加示蹤劑的個數。

5、如果注入堿且原油中有石油酸,組分數為10加示蹤劑個數。

二、化學離子設計

當注入液中加入堿以后,油藏流體中將會有一系列的化學反應發生,這些反應包括水相中的化學反應,固體的沉淀溶解反應、巖石表面的離子交換反應和表活劑上的離子交換反應。為了描述這些化學反應,CDNS設計了下列化學離子成份。

1、水相中的化學離子

⑴ H+

⑵ Na+

⑶ Ca2+ (Mg2+)

⑷ CO32-

⑸ HAo(油里的石油酸)

⑹ H2O

(以上6種是獨立的離子)

⑺ HCO3-

⑻ OH-

⑼ A-

⑽ HAw(水里的石油酸)  

2、固體化學離子

⑴ CaCO3

⑵ Ca (OH)2  

3、在巖石表面上吸附的離子

⑵()

三、CDNS軟件運行流程

1、CDNS共有兩個程序:

①CDNS1

②CDNS2

CDNS1是初始化程序,CDNS2是主模型模擬計算程序。

CDNS1需要一個用戶準備的輸入數據文件(數據流)。CDNS2需要的輸入數據從CDNS1運行后的輸出文件得到。

2、CDNS運行流程

四、數據流的填寫方法

1、關鍵字

CDNS的全部輸入數據部由關鍵字進行引導。共有18個一級關鍵字,它們是:

(1)TTT

(2)DATE

(3)TMAX

(4)UNIT

(5)INJECT

(6)GRIDMD

(7)GRIDVL

(8)WELLVL

(9)PARA

(10)A-DATA

(11)S-DATA

(12)P-DATA

(13)T-DATA

(14)SLUG

(15)DT

(16)ENDFW

(17)PERFOR

(18)TIME

(19)STOP

2、關鍵字說明

(1)關鍵字引導輸入的各類數據均采用自由F格式。二數之間至少相距一個空格。

(2)若同一關鍵字在數據流文件出現多次,新的將代替老的。

(3)當某些關鍵字在數據流文件中缺失時,程序自動對相應的數據置缺省值。

3、關鍵字的次序

CDNS對關鍵字之間的次序沒有限制,但從模擬計算的需要出發,一般次序應當是:

TTT(*)

DATE

TMAX

UNIT(*)

INJECT

GRIDMD

GRIDVL或WELLVL

PARA

A-DATA(*)

S-DATA(*)

P-DATA(*)

T-DATA(*)

SLUG

DT

ENDFW   (*)

PERFOR(*)

TIME

STOP

其中有(*)的關鍵字是可選擇的。

五、關鍵字詳細描述

1、TTT

格式:

TTT

ttt

說明:

ttt是標題內容,它是一個不超過50個字符的字符串。

例:

TTT

Gudong   52+3

2、DATE

格式:

DATE nn yy  dd  

說明:  

開始模擬的日期或油田投產日期

例:

DATE  1998  07  11

3、TMAX

格式:

TMAX   tmax   ooip

說明:

tmax-此次模擬的最大天數

ooip-原始地質儲量

單位:104m3

例:

TMAX  5000   1538

4、UNIT

當采用英制單位時,才使用此關鍵字,缺省單位是公制單位

格式:

UNIT   k

說明:

k――k=0 公制單位

k=1 英制單位

例:

UNIT 0

5、INJECT

格式:

INJECT ka ks kp kgel  kfoam  kt kgas kacid

說明

ka——ka=0  不注堿

ka=1  注堿,

ks——ks=0  不注表活劑

ks=1  注表活劑

kp——kp=0  不注聚合物

kp=1  注聚合物

kgel——kgel=0  不注GEL組分

kgel=1  注GEL組分

kfoam——kfoam=0  不注FOAM組分

kfoam=1  注FOAM組分

kt——kt=0  不注示蹤劑

kt=1  注示蹤劑

kgas——kgas=0  無氣相

kgas=1  有氣相

kacid——kacid=0  無石油酸

kacid=1  有石油酸

例:

INJECT  1  1  1  0   0   0   0  1

6、GRIDMD

格式:

GRIDMD  k1  k2

說明:

k1――坐標系統標識

k1=1 直角坐標系統

k1=2 柱坐標系統

k1=3 變寬度剖面坐標系統

k2――劃分網格方法標識

k2=1 人工劃分網格

k2=2 程序自動劃分網格

例:

GRIDMD  1  1

6.1如果k2=1,在GRIDMD行后面有下列數據行:

(1)數據行1

格式:

nx  ny  nz  nw   nwa

說明:

nx――x方向的網格塊數

ny――y方向的網格塊數

nx――z方向的網格塊數

nw――生產井數(包括注入井和采出井)

nwa――資料井數。如果nwa=0則全部生產井也是資料井。

例:

10    10     3     2    0

(2)數據行2

格式:

dxi  (i=1,nx)

說明:

dxi-x方向第i個網格塊的步長,

單位:m或feet。

例:(nx=10)

10*50.0

(3)數據行3

格式

dyi   (i=1,ny)

說明:

dyi-y方向第i個網格塊的步長,

單位:m或feet。

例:(ny=10)

3*50   20   30   5*40.0

(4)生產井信息行,每口井一行,共nw行。

格式:

wname iw  jw  flag  rw   swell

idir  ifirst  ilast  pwfmin  pwfmax   qmin   qmax

說明:

wname-井的名字,不超過10個字符

iw-井位在x方向上的網格號

jw-井位在y方向上的網格號

iflag-井別標識

iflag=1,定產注入井

iflag=2,定壓生產井

iflag=3,定壓注入井

iflag=4,定產生產井

rw-井筒的半徑

單位:m或feet。

swell-井的表皮因子

idir-井筒方向標識

idir=1,沿x方向的水平井

idir=2,沿y方向的水平井

idir=3,垂直井

ifirst-井射開層段的第一個網格塊號

filast-井射開層段的最后一個網格塊號

pwfmin-最小井底壓力限制

單位:Kpa 或pai

pwfmax-最大井底壓力限制

單位:Kpa 或pai

qwin-最小速度限制

單位:m3/day或ft3/day

qmax-最大速度限制

單位:m3/day或ft3/day


例:

A-15 5 3 1 0.15 0 3 1  3 5000 15000 0 500

(5)資料井信息行。如果nwa>0,則有nwa個資料井信息行,每口資料井一行,共nwa行。

格式:

wname   iw  jw

說明

wname-井名,不超過10個字符

iw-資料井在x方向的網格塊號

jw-資料井在y方向的網格塊號

例:B-301   8   17

6.2如果k2=2,在GRIDMD行后面有下列數據行:

(1)數據行1

格式

xoo  yoo  nz   nw   nwa   nbou

說明:

xoo  yoo-坐標系原點的坐標

單位:m或feet。

nz-Z方向上的網格數(層數)

nw-生產井數(包括注入井和采出井)

nwa-資料井數。若nwa為0,則全部生產井是資料井

nbou-邊界線的點數

例:

0  0  5  11  0  20

(2)生產井信息行,每口井一行,共NW行。

格式:

wname  xw   yw  iflag  rw  swell  idir  dfirst  dlast  pwfmin   pwfmax   qmin  qmax

說明:

wname-井名,不超過10個字符

xw,yw-井在x-y平面上的坐標

單位:m或feet。

dfirst-井射開段頂部的深度

單位:m或feet。

dlast-井射開段底部的深度

單位:m或feet。

其它與5.1中生產井信息行的說明相同。

例:

A-487  70  25  4  0.15  0  3  1250  1380  8000  16000  0  400

(3)資料井信息行

如果nwa>0,每口資料井一行,共nwa行。

格式:

wname  xw  yw

說明:

wname-井名

xw,yw-井的坐標

單位:m或feet。

例:

GD-33   84   35

(4)邊界點信息行

若nbou>0,則有nbou個邊界點信息行,共nbou行

格式:

xb  yb

說明

xb,yb-邊界點的坐標

單位:m或feet。

邊界點按逆時針方向排序,第一個邊界點與第nbou個邊界點重合。

例:

15  3

7、GRIDVL

這個關鍵字用于輸入地質屬性在網格塊上的分布值。

GRIDVL用雙關鍵字方法輸入數據。雙關鍵字的第一個關鍵字是地質屬性的名字,它可以是下列名字中的任何一個:

DEPTH(深度)

DZ (網格厚度)

DZNET(凈厚度)

POR(孔隙度)

PERMX(x方向上的滲透率)

PERMY(y方向上的滲透率)

PERMZ(z方向上的滲透率)

PW(初始水相壓力)

SW(初始水飽和度)

雙關鍵字的第二個關鍵字是賦值方法,它可以是下列關鍵字之一:

CC   XX   YY   ZZ   VV   L1   DIP   MLT   MOD

說明:

(1)CC

表明前面的數組網格值是一個常數。

例:

POR   CC  0.33

(2)xx

表明數組值只在x方向上變化,應輸入nx個值。

例:

nx=10

PERMX  xx

1500 1400 1300 1260 1250 1200 1000 950 900 850

(3)yy

表明數組值只在y方向上變化,應輸入ny個值。

例:(ny=5)

PERMX  YY

1400   1800   1900   2000  2100

(4)ZZ

表明數組值只在Z方向上變化,應輸入nz個值。

例:(nz=3)

DZNET   ZZ

4.0    5.5   3.2

(5)VV

表明輸入全部網格塊上的值

例:(nx=10,ny=5,nz=2)

DEPTH  VV

1500  1510  1520  1530  1540  1550  1560  1570  1580  1590

1510  1520  1530  1540  1550  1560  1570  1580  1590  1600

1520  1530  1540  1550  1560  1570  1580  1590  1600  1610

1530  1540  1550  1560  1570  1580  1590  1600  1610  1620

1540  1550  1560  1570  1580  1590  1600  1610  1620  1630

1505  1515  1525  1535  1545  1555  1565  1575  1585  1595

1515  1525  1535  1545  1555  1565  1575  1585  1595  1605

1525  1535  1545  1555  1565  1575  1585  1595  1605  1615

1535  1545  1555  1565  1575  1585  1595  1605  1615  1625

1545  1555  1565  1575  1585  1595  1605  1615  1625  1635

(6)L1

表明只輸入第一層的網格值,其它層由程序自動計算。L1只適用于第一個關鍵字是DEPTH(深度)。

例:(nx=10  ny=5  nz=3)

DEPTH    L1

50*1230

(7)DIP

DIP只用于DEPTH(深度),表明只輸第一個網格塊的值,其它網格塊的深度由程序自動計算。

格式:

DEPTH    DIP   d111     xdd    ydd

說明:

d111-第一個網格塊(1,1,1)的深度

單位:m或feet

xdd-x方向上的傾角

單位:度

ydd-y方向上的傾角

單位:度

例:(nx=10,ny=5,nz=3)

DEPTH     DIP     1230    10     7.5

(8)MLT

表明該數組的值是另一個數組的倍數

格式:

數組名1   MLT    fact    數組名2

說明:

fact-倍數因子

例:

PERMY    MLT      1.0    PERMX

(9)MOD

MOD用于對數組一部分網格上的值進行修改,其它值不動。

格式:

數組名   MOD  i1  i2  j1  j2  k1  k2  #   bb

說明:

i1-x方向起始網格號

i2-x方向結束網格號

j1-y方向起始網格號

j2-y方向結束網格號

k1-z方向起始網格號

k2-z方向結束網格號

#-算術運算符,它是“+”、“-”、“*”、“/”、“=”之一。

bb-數值

例:PERMX  MOD 2 3 3 4 1   1 = 1500

8、WELLVL

當地質屬性參數只在井點有值時,需要用WELLVL來輸入井點的地質參數,然后程序自動插值計算網格塊上的值。

格式:

WELLVL  intpl

說明:

intpl-插值方法標識

intpl=1 改進的距離加權插值法

intpl=2 克里金插值法

intpl=3  分形插值法

例:

WELLVL  1

8.1資料井數據信息行

關鍵字WELLVL行的后面跟有nwa口資料井的數據信息,每口井有nz行,共有nwa*nz行。

格式:

jh  kz  wdepth  wdz  wpor  wkx   wky   wkz  wpw   wsw

說明:

jh-資料井號,不超過10個字符

kz-層號,1≤kz≤nz

wdepth-該井層頂深

單位:m或feet

wdz-該井層的厚度

單位:m或feet

wpor-該井層的孔隙度

單位:分數

wkx-該井層x方向上的滲透率

單位:10-3um2

wky-該井層y方向上的滲透率

單位:10-3um2

wkz-該井層z方向上的滲透率

單位:10-3um2

wpw-該井層的初始水相壓力

單位:kpa或ps;

wsw-該井層的初始水飽和度

單位:分數

例:

A8 1 1200 3.5   0.35   1300  1300   130  12000  0.45

9、PARA

關鍵字PARA引導用戶輸入模擬需要的一些關于油藏的常數,包括巖石、水、油的物理性質、地層水和注入水的化學成份、相對滲透率曲線和毛管壓力曲線等。

常數用下列關鍵字輸入,這些關鍵字是可選擇的,如果沒有輸入,相應的常數將取缺省值。

ROCK

WATER

OIL

KRHI

KRLO

PC

FWCHEM

IWCHE1

IWCHE2

9.1  ROCK

用于輸入巖石的壓縮系數和油藏的標準條件

格式:

ROCK  compcr   pstand   tstand

說明:

    compcr-巖石的壓縮系數

單位:1/kpa 或1/psi

pstand-標準條件下的參照壓力

單位:Kpa 或psi;

tstand-標準條件下的參照溫度

單位:℃或℉

例:

ROCK  0.00001  12000   40.0

9.2 WATER

輸入水的物理性質

格式:

WATER   denw   visw  conpcw  xiftw

說明:

denw-水的密度

單位:g/cm3或psi/ft

visw-水的粘度

單位:mPa.s

compcw-水的壓縮系數

單位:1/kpa 或1/psi

xiftw――油水界面張力的對數值

例:

WATER   0.99  0.45  0.00005   1.3

9.3 OIL

輸入油的物理性質

格式:

OIL     deno   viso  conpco   acidv  

說明:

deno-油的密度

單位:g/cm3或psi/ft

viso-水的粘度

單位:mPa.s

compco-油的壓縮系數

單位:1/kpa 或1/psi

acidv-石油酸的酸值

單位:無量綱

例:

OIL  0.92  70.0  0.0001   1.7

9.4 KRHI

輸入高界面張力時的相對滲透率曲線

格式:

KRHI  S1r   S2r   P1r   P2r   e1   e2

說明:

S1r-水相的束縛水飽和度

單位:分數

S2r-油相的殘余油飽和度

單位:分數

P1-水相相對滲透率的最大油

單位:分數

P2-油相相對滲透率的最大值

單位:分數

e1-水相相對滲透率曲線的指數

單位:無量綱

e2-油相相對滲透率曲線的指數

單位:無量綱

例:

KRHI  0.30  0.20  0.2  0.95  1.8  2.3

9.5  KRLO

輸入低界面張力時的相對滲透率曲線。

格式:

KRLO  s1rc  s2rc  p1rc  p2rc  e1c  e2c

說明:

   s1rc——水相的束縛水飽和度

          單位:分數

   s2rc——油相的殘余油飽和度

           單位:分數

   p1rc——水相相對滲透率的最大值

           單位:分數

   p2rc——油相相對滲透率的最大值

           單位:分數

   e1c——水相相對滲透率曲線的指數,無量綱

   e2c——油相相對滲透率曲線的指數,無量綱

例:

KRLO  0.  0.  0.8  1.0  1.1  0.9

9.6 PC

輸入毛管壓力曲線數據

格式:

PC iow  cpc  epc

說明:

iow-巖石的潤濕性標識

iow=0 親水

iow=1 親油

cpc-毛壓曲線端點值

單位:或

epc-毛壓曲線指數

單位:無量綱

例:

PC  0  62  2.1

9.7 FWCHEM

輸入地層水的化學組成

格式:

FWCHEM  C-CL  C-H  C-Na  C-Ca  C-Mg  C-CO3   C-HCO3

C-H2CO3   C-OH

說明

C-CL-地層水cl-的濃度

C-H-地層水H+的濃度

C-Na-地層水Na+的濃度

C-Ca-地層水Ca2+的濃度

C-Mg-地層水Mg2+的濃度

C-CO3-地層水CO32-的濃度

C-HCO3-地層水HCO3-的濃度

C-H2CO3-地層水H2CO3的濃度

C-OH-地層水OH-的濃度

單位:PPm (即mg/L)

例:

FWCHEM 5600 0.000006 400 80 23 37 220 110 13

9.8 IWCHE1

輸入第一種注入水的化學組成

格式:

IWCHE1   C-CL  C-H  C-Na  C-Ca  C-Mg  C-CO3   C-HCO3

               C-H2CO3   C-OH

說明:

C-CL-注入水cl-的濃度

C-H-注入水H+的濃度

C-Na-注入水Na+的濃度

C-Ca-注入水Ca2+的濃度

C-Mg-注入水Mg2+的濃度

C-CO3-注入水CO32-的濃度

C-HCO3-注入水HCO3-的濃度

C-H2CO3-注入水H2CO3的濃度

C-OH-注入水OH-的濃度

單位:ppm(即mg/L)

例:

IWCHE1 4300 0.0000001 230 120 55 41 60 95 18


9.9 IWCHE2

輸入第一種注入水的化學組成

格式:

IWCHE2   C-CL  C-H  C-Na  C-Ca  C-Mg  C-CO3   C-HCO3

            C-H2CO3   C-OH

說明:

C-CL-注入水cl-的濃度

C-H-注入水H+的濃度

C-Na-注入水Na+的濃度

C-Ca-注入水Ca2+的濃度

C-Mg-注入水Mg2+的濃度

C-CO3-注入水CO32-的濃度

C-HCO3-注入水HCO3-的濃度

C-H2CO3-注入水H2CO3的濃度

C-OH-注入水OH-的濃度

單位:ppm(即mg/L)

例:

IWCHE2 4300 0.0000001 230 120 55 41 60 95 18


10、A-DATA

關鍵字A-DATA引導用戶輸入關于堿的參數。

堿的數據用下列關鍵字輸入:

   A-TYPE

   A-PH

   A-ACID

10.1 A-TYPE

  格式:

 A-TYPE itype

說明:

 itype——堿的類型標志

itype=1    NaOH

itype=2    NaCO3

10.2 A-PH

格式:

 A-PH phc pht

說明:

 phc——PH的臨界值,大于該值后表活劑吸附依賴于PH值

 pht——PH的臨界值,大于該值后表活劑不發生吸附

例:

A-PH     7.0   11.2

10.3  A-ACID

格式:

   A-ACID  cselp  cseup  cac2i  acidis  eqwps

說明:

   cselp——能夠產生表活劑的含鹽量低限

            單位:meq/ml

   cseup——能夠產生表活劑的含鹽量高限

            單位:meq/ml

   cac2i——油里的初始石油酸濃度

            單位:moles/l 油

acidis——石油酸的溶解系數

eqwps——石油酸的當量

例:

A-ACID  0.02  0.25  0.003  0.006  500


11、S-DATA

關鍵字S-DATA引導用戶輸入關于表活劑的參數。

表活劑的數據用下列關鍵字輸入:

  IFTTAB

  S-DEN

  S-ADSO

  S-EXCH

  S-CSEL

11.1 IFTTAB

格式:

    IFTTAB  nift  cseift  

說明:

     nift——實驗室界面張力曲線的數據點數

     cseift——測試時含鹽量

           單位:mg/l

     如果nift>0,則后面跟有nift個數據行。

格式:

xift   yift   zift

說明:

xift――堿的濃度

單位:wt%

yift――表活劑濃度

單位:wt%

zift――界面張力

單位:mN/m

例:

IFTTAB         6

0                0              20

0                0.11            5

0                0.1             0.5

0                0.2             0.01

0                0.3             0.002

0                0.4             0.001


11.2 S-DEN

格式:

     S-DEN    den3    compc3

說明:

     den3——表活劑的密度

           單位:g/m3或psi/ft

     compc3——表活劑的壓縮系數

           單位:1/kPa或1/psi

例:

S-DEN   0.99   0.0

11.3  S-ADSO

輸入吸附參數

此關鍵字有兩種格式任選其一 。

格式1:

    S-ADSO   1   ephi3   nsads   csesad    

說明:

   ephi3-表面活性劑的可及體積分數

nsads-吸附曲線的數據點數目

csesad-測吸附曲線的含鹽量    單位:mg/L

如果nsads>0則后面跟有nsads個數據行,

格式:

xsads      ysads

說明:

xsads-表面活性劑濃度     單位mg/L

ysads-表面活性劑的吸附量      單位mg/g

例:S-ADSO   1   1.0   15   5000  

0.0  0.0

500 0.15

10000.21

15000.25

2000           0.31

3000            0.4

4000               0.52

5000               0.64

60000.75

7000               0.86

8000               0.98

9000               1.02

10000              1.02

11000              1.02

12000              1.02

格式2:

    S-ADSO  2  ephi3   ad31   ad32   b3d    

說明:

     ephi3 -表活劑的可及體積分數

     ad31——表活劑的吸附參數,a31

           單位:體積比

     ad32——表活劑的吸附參數,a32

           單位:ml/meq

     b3d——表活劑吸附參數,b3

          單位:1/表活劑體積

     例:

   S-ADSO  2  1.0   0.6  0.25  1010  

11.4  S-EXCH

輸入離子交換參數

格式:

    S-EXCH   qv   xkc   xks   eqw

說明:

     qv——粘土的離子交換能力

         單位:(meq/ml)/pv

     xkc——粘土的陰離子交換常數,c

          單位:(meq/ml)-1

     xks——表活劑的陰離子交換常數,s

          單位:(meq/ml)-1

     eqw——表活劑的當量

例:

     S-EXCH   0.05   0.25   0.2   419

11.5  S-CSEL

輸入含鹽量的參數

格式:

   S-CSEL   csel   cseu   6

說明:

   csel——進入Ⅲ相的最低的含鹽量

          單位:meq/ml

   cseu——進入Ⅲ相的最高的含鹽量

          單位:meq/ml

   6——關于鈣的含鹽量的斜率參數,無量綱

例:

   S-CSEL   0.012   0.54   0.8



12、P-DATA

關鍵字P-DATA引導用戶輸入有關聚合物的參數。在P-DATA行后面用下列關鍵字輸入聚合物的具體參數:

VISTAB

P-GAMM

P-REDU

P-ADSO

12.1  VISTAB

關鍵字P-VIST用于輸入一種聚合物的粘濃曲線表。

格式:

VISTAB  nvis      csevis    

說明:

nvis-聚合物粘濃曲線的數據點數

csevis-測試時的含鹽量

單位:mg/l

如果nvis>0,則后面跟有nvis個數據行。

格式:

xvis         yvis         zvis

說明:

xvis-堿的濃度

單位:重量百分數

yvis-聚合物的濃度

單位:ppm(即mg/l)

zvis-聚合物溶液的粘度

單位:mpa.s

例:

VISTAB  6         5000

0    0             0.45

0    100              3.0

0    500              8.4

0    1000             12.5

0    1500             16.8

0    2000             25.1

12.2  P-GAMM

關鍵字 P-GAMM用于輸入聚合物的剪切參數。

格式:

P-GAMM    gammac     gamhf     pown

說明:

gammc-剪切方程中的系數

單位:

gamhf-零剪切速率時的粘度半值時的剪切速率

單位:sec-1

Pown-聚合物粘度隨剪切速率變化的指數

單位:無量綱

例:

P-GAMM     3.2      20     1.1

12.3 P-REDU

關鍵字P-REDU用于輸入滲透率降低的參數。

此關鍵字有兩種格式任選其一。

格式1:P-RERU  1  npredu   csered

說明:

npredu-滲透率下降曲線數據點數目

csered-測試滲透率下降曲線的含鹽量    單位:mg/L

如果npredu>0,后面跟有npredu個數據行。

格式:

   xpredu ,    rk

說明:

    xpredu-聚合物濃度     單位:wt%

    rk-滲透率下降系數     單位:無量綱

例:

  P-REDU  1   6    5000

0.01.0

0.05       3.0

0.10       5.0

0.15       6.1

0.20       7.0

0.25       7.9

格式2:

P-REDU   2   brk   crk

說明:

brk-計算滲透率降低因子RK的參數

單位:體積分數/重量百分數

crk-計算Rk的參數

單位:

例:

P-REDU   2   1000  0.0186

12.4  P-ADSO

關鍵字P-ADSO用于輸入聚合物的可及體積和吸附方面的參數。

此關鍵字有兩種格式任選其一。

格式1:

P-ADSO    1    ephi4   npads    csepad    

說明:

   ephi4 -聚合物的可及體積分數

npads-吸附曲線的數據點數

   csepad-側吸附曲線的含鹽量   單位:mg/L

   如果npads>0,則后面跟有npads個數據行,

格式:

xpads      ypads

說明:

   xpads-聚合物濃度     單位:mg/L

   ypads-聚合物吸附量     單位:mg/g

例:

  P-ADSO    1    0.8    8    5000  

0.0 0.0

2500.02

500          0.04

7500.065

10000.07

12500.07

15000.071

1750        0.071

格式2:

P-ADSO  2    ephi4   ad41    ad42   b4d

說明:

   ephi4 -聚合物的可及體積分數

單位:分數

ad41-聚合物的吸附參數

單位:無量綱

ad42-聚合物的吸附參數

單位:ml/meq

b4d-聚合物的吸附參數

單位:體積分數/重量百分數

例:

P-ADSO    0.85    0.7   0   100   




13、T-DATA

關鍵字T-DATA引導用戶輸入關于示蹤劑的參數。在T-DATA行后面用下列關鍵字輸入具休數據:

   T-NT

T-NAME

T-TK

T-TKS

T-RDC

T-RET

13.1 T-NT

   用于輸入示蹤劑的

  T—NT 用于輸入示蹤劑的數目.

  格式:

  T-NT   nt

   說明:

          nt ----- 示蹤劑數目

               nt = 1,2,3

13.2 T-NAME

T-NAME用于輸入示蹤劑的名字

格式:

T-NAME    tname i (i=1,NT)

說明:

tanme i-第i個示蹤劑的名字不超過10個字符。

例:(NT=3)

T-NAME   TRACER1   TRACER2   TRACER3  

13.3 T-TK

T-TK用于輸入示蹤劑在油水中的分配系數

格式:

T-TK  tki  (i=1,NT)

說明:

tki-第i種示蹤劑在初始條件下在水中的分配系數,

如果tki=0,表示是水示蹤劑

如果tki=-1,表示是油示蹤劑

例:(NT=1)

T-TK    0

13.4 T-TKS

T-TKS用于輸入計算示蹤劑分配系數的參數,此參數作為含鹽量變化的系數來計算分配系數的變化。

格式:

T-TKS   tksi   (i=1,NT)

說明:

tksi-第i個示蹤劑的參數

單位:(meq/ml)-1


例:(NT=1)

T-TKS  0.003

13.5  T-RDC

T-RDC用于輸入示蹤劑的放射性衰減系數。

格式:

T-RDC   rdci (i=1,NT)

說明:

rdci-第i個示蹤劑的放射性衰減系數

單位:1/day

例:(NT=1)

T-RDC   0.03

13.6 T-RET

T-RET用于輸入示蹤劑的滯后因子。

格式:

T-RET  reti   (i=1,NT)

說明:

reti-第i個示蹤劑的滯后因子

單位:無量綱

例:(NT=1)

T-RET   0.061

14、SLUG

關鍵字SLUG引導用戶輸入動態井數據(工作制度)。當工作制度改變或注入濃度改變時就要用SLUG輸入一次。因此SLUG行將在數據流中多次出現。

SLUG行后緊跟nw個數據行,每口井占一行。數據行對于不同的井別有不同的格式,共有4種格式:

格式1:(用于定產注入井,iflag=1)

iwd  q  cca  ccs  ccp   iiw

GELCC   cgel1  cgel2  cgel3  cgel4(*)

FOAMCC  cfoam1  cfoam2  cfoam3 (*)

TRACERCC  cti (I=1,NT) (*)

說明:

iwd-井的序號

q-井的注入速度

單位:m3/day或ft3/day

cca-堿的注入濃度

單位:wt%

ccs-表活劑的注入濃度

單位:wt%

ccp-聚合物的注入濃度

單位:wt%

iiw-注入水號

iiw=1 或 2

       cgell1,cgel2,cgel3,cgel4 是GEL各組分的濃度,如沒有GEL注入,則沒有GELCC這一行.

       Cfoam1,cfoam2,cfoam3 是FOAM各組分的濃度, 如沒有FOAM注入, 則沒有FOAMCC這一行.

cti-第i種示蹤劑的注入濃度

單位:wt%

       如沒有示蹤劑,即NT=0, 則沒有TRACERCC這一行.

例:

1  200  1.2  0.3  0.15

格式2:(用于定壓注入井,iflag=3)

iwd   pwf   cca  ccs  ccp   iiw

GELCC   cgel1  cgel2  cgel3  cgel4(*)

FOAMCC  cfoam1  cfoam2  cfoam3 (*)

TRACERCC  cti (I=1,NT) (*)

說明:

iwd-井的序號

pwf-井底壓力

單位:kpa 或psi

cca  ccs  ccp  等,同格式1

例:

2   12000   1.2   0.3  0.18

格式3:(用于定產采出井,iflag=4)

iwd  q

說明:

iwd-井的序號

q-井的產液速度

單位:m3/day或ft3/day

例:

5     106

格式4:(用于定壓采出井,iflag=2)

iwd  Pwf

說明:

iwd-井的序號

pwf-井底壓力

單位:kpa或psi

例:

8  10100

15、DT

關鍵字DT用于輸入時間步長控制,它在數據流可出現多次。

格式:

DT  kdt   dto   dclim   dtmin   dtmax

說明:

kdt-時間步長控制方式標識

   kdt=1  定步長

kdt=2    變步長

dto-當kdt=1時,表示定長長值

當kdt=2時,表示初始步長值

單位:day

delim-每時間步濃度變化限制

單位:分數

dtmin-最小的時間步長限制

(kdt=2)

單位:day

dtmax-最大的時間步長限制

(kdt=2)

單位:day

例:

   DT  1   0.5

DT  2   0.05   0.005   0.01   5

16、ENDFW

關鍵字ENDFW輸入最大含水值,用以控制模擬結束時間.。

格式:

ENDFW  endfw

說明:

endfw-含水值.

例:

ENDFW     0.98

17、PERFOR

關鍵字PERFOR輸入井的射孔信息,它可以多次出現。如果沒有PERFOR行,表明井對所有層射孔。

一個PERFOR行輸入一口井的射孔信息。這口井改變射孔時再用PERFOR行。

格式:

PERFOR  wname   Ki (i=L1,L2)

說明:

wname-井的名字

Ki-第i段的射孔標志

Ki=0不射孔

Ki=1射孔

L1-井的起始段

L2-井的結束段

例:(NZ=4)

PERFOR   A-71   1  1  0  1

18、TIME

關鍵字TIME輸入用戶需要輸出模擬計算結果的時間,它應當出現許多次。

格式:

TIME   t

說明:

t-輸出模擬計算結果的時間

單位:day

例:

TIME  30

TIME  1500

19、STOP

STOP是數據流結束的關鍵字

格式:

STOP

例:

STOP



六、數據流的例子


C --- TITLE

TTT

            GUDAO  A-S-P FLOOD TEST

C  DATE   YERS  MONTH  DAYS

DATE    1997    10   1

C

TMAX   15000    52.5

C

C --- METRIC

METRIC

C

C --- INJEC  k-a  k-s  k-p  k-t  k-g  k-acid

INJECT  1  1  1  0  0  0  0  1

C

C --- GRIDMD  kcoord  kgrid

GRIDMD  1   1

C

C *** nx ny nz nw

     10   10   2   2   0    0

C *** dx(i) i=1 nx

     10*50.

C *** dy(i) i=1 ny

     10*50.

C *** wname(m) iw(m) type(m) rw(m) swell(m) dir(m)

C     ifirst(m) ilast(m) prf(m) pwfmin(m) pwfmax(m)

C     qtmin(m) qtmax(m)         m=1 nw

C

AA1   2   2   1   .15  0.  3   1  2    0.0 125000. .0 500.

AA2   9   9   4   .15  0.  3   1  2    0.0 125000. .0 500.

C --- GRIDVL

GRIDVL

C

C

POR    VV

     200*0.3

C

DZNET   VV

    200*5.

C

DEPTH   L1

    100*1200.

C

PERMX    VV

   100*5000.     100*500.

C

C

PERMY  MLT  1.00  PERMX

C

C

PERMZ  MLT  0.01  PERMX

C

C

PW    VV

     100*12000.    100*12050.

C

C

SW     VV

     200*0.45

C

C

PARA

 ROCK   0.000001  12000.    0.

 OIL    0.96   70.0   0.000001   1.7

 WATER  1.00    0.46  0.000001   1.3

 KRHI   .30  .25   .25   .95  1.45  2.7

 KRLO   .0    .0    0.7   1.    1.95   0.9

 PC      0.      62.0689     2.0

C

C ---  1    2    3   4     5     6      7      8      9

C --- Cl-  H+  Na+  Ca++  Mg++  CO3--  HCO3-  H2CO3  OH-

C

 FWCHEM  6000.  0.000002   100.   99.   20.   20.  700.   0.   0.

C

C --WCI1  CL-      H+            Na+    Ca++  Mg++  CO3--HCO3-H2CO3 OH-(ppm)

 IWCHE1   6000.  0.000002   100.   99.   20.   20.  700.   0.   0.

C --WCI2  CL-      H+            Na+    Ca++  Mg++  CO3--HCO3-H2CO3 OH-(ppm)

 IWCHE2   6000.  0.000002   100.   99.   20.   20.  700.   0.   0.

C

C

C

C

C   --- SURFAC

S-DATA

IFTTAB  36

  0.       0.       20.

  0.       0.1      0.5

  0.       0.2      0.08

  0.       0.3      0.006

  0.       0.4      0.0009

  0.       0.5      0.0008

  0.5      0.       5.

  0.5      0.1      0.2

  0.5      0.2      0.05

  0.5      0.3      0.004

  0.5      0.4      0.0007

  0.5      0.5      0.0006

  1.       0.       3.

  1.       0.1      0.015

  1.       0.2      0.02

  1.       0.3      0.002

  1.       0.4      0.0005

  1.       0.5      0.0004

  1.5      0.       1.

  1.5      0.1      0.10

  1.5      0.2      0.009

  1.5      0.3      0.0009

  1.5      0.4      0.0003

  1.5      0.5      0.0002

  2.0      0.       0.5

  2.0      0.1      0.05

  2.0      0.2      0.006

  2.0      0.3      0.0006

  2.0      0.4      0.0001

  2.0      0.5      0.00009

  3.0      0.       0.1

  3.0      0.1      0.01

  3.0      0.2      0.003

  3.0      0.3      0.0003

  3.0      0.4      0.0001

  3.0      0.5      0.0001

C

C

C   S-CMC  0.    1.    0.00008

C   S-IFT  13.  -16.5   .007  13.   -16.1  .010

C   S-TRAP      1865.      59074      364.2

S-DEN    0.98     0.

C  S-ALPH    0.5     0.5     0.     0.9     0.7

C  S-PC      0.      62.0689     2.0

C

S-ADSO     1     1.0    15    5000.

          0.         0.

        500.         0.15

       1000.         0.21

       1500.         0.25

       2000.         0.31

       3000.         0.4

       4000.         0.52

       5000.         0.64

       6000.         0.75

       7000.         0.86

       8000.         0.98

       9000.         1.02

      10000.         1.02

      11000.         1.02

      12000.         1.02

C

S-EXCH    0.05  .25   .2   419.

S-CSEL    0.012    0.544   0.8

C

C

C

C

C --- POLYM

P-DATA

C

VISTAB       18

  1000.     0.         0.5

  1000.   200.         4.5

  1000.   600.        18.

  1000.  1000.        35.

  1000.  1500.        44.

  1000.  2000.        50.

  4000.     0.         0.5

  4000.   200.         4.5

  4000.   600.        17.

  4000.  1000.        34.

  4000.  1500.        43.

  4000.  2000.        49.

  7000.     0.         0.5

  7000.   200.         4.4

  7000.   600.        16.

  7000.  1000.        33.

  7000.  1500.        42.

  7000.  2000.        48.

C

C

C   P-AP       60.  1600.   0.    2.   .01   -0.175

P-GAMM     3.219      20.    1.1

C   P-REDU     1000.  0.0186

P-REDU   1    5

    0.0       1.0

 1000.        1.2

 1500.        1.5

 2000.        2.0

 3000.        2.5

C

P-ADSO    1     0.8     8    5000.

   0.0        0.0

 250.         0.02

 500.         0.04

 750.         0.065

1000.         0.07

1250.         0.07

1500.         0.071

1750.         0.071

C

C

C

C --- ALKAL

A-DATA

A-TYPE 2

C   A-PH    7.22  11.25  10.06  0.

A-PH    7.22  11.25

A-ACID  0.30  0.41    1.7    0.00000001      500.

C

C

C

C --- DT  kdt dt dclim dtmin dtmax

DT     2   0.05  0.005   0.005  5.0

C

C --- ENDFW    endfw

ENDFW     0.98

C

SLUG

C *** id(m) qt(m) c(1 m) c(2 m) .... c(n m)    for flag(m)=1

C *** id(m) pwf(m)                             for flag(m)=2

C *** id(m) pwf(m) c(1 m) c(2 m) .... c(n m)   for flag(m)=3

C *** id(m) qt(m)                              for flag(m)=4

C ***        m=1 nw

1  100.   0.  0.   0.   1    

2     100.

C

C

C --- TIME   time (days)

TIME   1.

TIME   2.

TIME   5.

TIME   30.

C

C

C --- SLUG

SLUG

C *** id(m) qt(m) c(1 m) c(2 m) .... c(n m)    for flag(m)=1

C *** id(m) pwf(m)                             for flag(m)=2

C *** id(m) pwf(m) c(1 m) c(2 m) .... c(n m)   for flag(m)=3

C *** id(m) qt(m)                              for flag(m)=4

C ***        m=1 nw

1  100.   0.   0.  .05      1

2     100.

C

C

C --- TIME   time (days)

TIME   60.

TIME   90.

TIME   120.

TIME   150.

TIME   200.

TIME   250.

TIME   300.

C

C

C --- SLUG

SLUG

C *** id(m) qt(m) c(1 m) c(2 m) .... c(n m)    for flag(m)=1

C *** id(m) pwf(m)                             for flag(m)=2

C *** id(m) pwf(m) c(1 m) c(2 m) .... c(n m)   for flag(m)=3

C *** id(m) qt(m)                              for flag(m)=4

C ***        m=1 nw

1  100.  1.2    0.3    0.15      1

2     100.

C

C

C

TIME   350.

TIME   400.

TIME   450.

TIME   500.

TIME   550.

TIME   600.

TIME   650.

TIME   700.

TIME   750.

TIME   800.

TIME   850.

TIME   900.

TIME  950.

TIME  1000.

TIME  1050.

TIME  1100.

TIME  1200.

C

C --- SLUG

SLUG

C *** id(m) qt(m) c(1 m) c(2 m) .... c(n m)    for flag(m)=1

C *** id(m) pwf(m)                             for flag(m)=2

C *** id(m) pwf(m) c(1 m) c(2 m) .... c(n m)   for flag(m)=3

C *** id(m) qt(m)                              for flag(m)=4

C ***        m=1 nw

1  100.    0.   0.   0.     1

2     100.

C

C

TIME  1300.

TIME  1400.

TIME  1500.

TIME  1600.

TIME  1700.

TIME  1800.

TIME  1900.

TIME  2000.

C  STOP

C

C   DT     2   0.05  0.006   0.005  5.

C

TIME  2200.

TIME  2600.

TIME  3000.

TIME  3600.

ENDFW   0.98

TIME  4000.

TIME  4500.

TIME  5000.

TIME  5500.

TIME  6000.

TIME  7000.

TIME  8000.

TIME  9000.

TIME  10000.

TIME  11000.

TIME  12000.

TIME  13000.

TIME  14000.

TIME  15000.

C

C

STOP

END








七、結果輸出文件

用戶在用CDNS軟件模擬計算一個方案后,可以得到5類結果數據文件..

這些數據文件有下列用途:

1.檢查運行過程的正確性和錯誤;

2.油藏工程技術人員分析指標;

3.繪制各種曲線;

4.繪制各種等值圖和三維可視圖;

5.繪制流線圖.

1、數據檢查信息文件

CDNS輸出2個數據檢查信息文件:

    $1.check1

    $1.check2

其中, $1 表示模擬計算方案數據流文件的名字.

    $1.check1是初始化數據檢查信息文件.

    $1.check2是模型計算數據檢查信息文件.


2、指標分析信息文件

CDNS輸出3個指標分析信息文件:

    $1.report0

    $1.report1

    $1.report2

其中, $1 表示模擬計算方案數據流文件的名字.

    $1.report1是油田指標數據文件.

    $1.report2是井指標數據文件.

    $1.report0 是從$1.report1中提出的14項指標.

3、曲線繪制數據文件

CDNS輸出2個曲線繪制數據文件:

    $1.field

    $1.wells

其中, $1 表示模擬計算方案數據流文件的名字.

    $1. field是針對油田指標繪制曲線的數據文件.

    $1. wells是針對井指標繪制曲線的數據文件.


4、等值圖和三維可視圖繪制數據文件

CDNS輸出4個等值圖和三維圖繪制的數據文件:

    $1.stat

    $1.psvi

    $1.conc

    $1.chem

其中, $1 表示模擬計算方案數據流文件的名字.

    $1.stat是靜態數據場,包括深度,厚度,孔隙度,滲透率,初始水相壓力和初始水飽和度.

    $1. psvi是動態數據場,包括相壓力,相飽和度,相粘度和界面張力.

$1.conc 是動態數據場,包括各組分的濃度.

$1.chem. 是動態數據場,包括各化學離子的濃度.


5、流線圖繪制數據文件

CDNS輸出1個流線圖繪制的數據文件:

    $1.lx000              

其中, $1 表示模擬計算方案數據流文件的名字.

    $1.lx000記載了各時間階段各層流線圖的軌跡,可直接用于繪制流線圖.  

您可以通過以下方式聯系我們010-64846440
友情鏈接 LINK 中國石油 中國石化 中國海油 延長石油
內容:
姓名:
電話:
郵箱:
掃碼關注
公司地址:北京市海淀區西小口路66號D-3-408室
郵 編:100192
電 話:010-64846440 / 010-64868477
E-mail:jiahewuxian@163.com
版權所有@北京嘉和無限科技有限公司 京ICP備18027440號-1
成年免费A级毛片-日本欧美41禁毛片大片