石油工程開發(fā)課程設(shè)計(jì)報(bào)告

1、東 北 石 油 大 學(xué)課 程 設(shè) 計(jì)課 程 石油工程課程設(shè)計(jì) 題 目 瑞利(Rayleigh)模型的應(yīng)用 院 系 石油工程學(xué)院 專業(yè)班級(jí) 油工班 學(xué)生姓名 學(xué)生學(xué)號(hào) 指導(dǎo)教師 王立軍 2013年7月19日32 / 33文檔可自由編輯打印東北石油大學(xué)課程設(shè)計(jì)任務(wù)書課程 石油工程課程設(shè)計(jì)題目 瑞利(Rayleigh)模型的應(yīng)用專業(yè) 石油工程 姓名 學(xué)號(hào) 主要內(nèi)容、基本要求、主要參考資料等主要內(nèi)容：(1) 推導(dǎo)瑞利預(yù)測(cè)模型；(2) 根據(jù)油氣田實(shí)際生產(chǎn)數(shù)據(jù)，進(jìn)行線性回歸，求得直線的截距和斜率；(3) 根據(jù)相關(guān)公式，確定瑞利模型常數(shù)a和c；(4) 計(jì)算油田年產(chǎn)油量；(5) 計(jì)算油田累積產(chǎn)油量；(6) 計(jì)

2、算可采儲(chǔ)量；(7) 計(jì)算最高年產(chǎn)量發(fā)生的時(shí)間；(8) 計(jì)算最高年產(chǎn)量；(9) 計(jì)算油田達(dá)到最高產(chǎn)量時(shí)的累積產(chǎn)量；(10) 繪制油田實(shí)際年產(chǎn)量與預(yù)測(cè)產(chǎn)量對(duì)比曲線；(11) 繪制油田實(shí)際累積產(chǎn)量與預(yù)測(cè)累積產(chǎn)量對(duì)比曲線?；疽螅?) 基礎(chǔ)數(shù)據(jù):2) 某油田或區(qū)塊的開發(fā)數(shù)據(jù)。2)設(shè)計(jì)要求該專題設(shè)計(jì)最終要求是，學(xué)生通過(guò)自選上述基礎(chǔ)數(shù)據(jù)，利用所學(xué)知識(shí)完成規(guī)定設(shè)計(jì)內(nèi)容，編制相應(yīng)軟件，并提交規(guī)范設(shè)計(jì)報(bào)告。主要參考資料陳元千 油藏工程實(shí)踐M.石油工業(yè)出版社. 完成期限 2013年7月1日2013年7月19日 指導(dǎo)教師 王立軍 專業(yè)負(fù)責(zé)人 王立軍 2013年7 月 19日目 錄1 前 言31.1 設(shè)計(jì)的目的意義

3、31.2 設(shè)計(jì)的主要內(nèi)容32 基礎(chǔ)數(shù)據(jù)43 基礎(chǔ)理論53.1瑞利(Rayleigh)模型的應(yīng)用模型的建立53.2模型的求解方法74 設(shè)計(jì)結(jié)果8認(rèn)識(shí)與結(jié)論12參考文獻(xiàn)13附 錄141 前 言1.1 設(shè)計(jì)的目的意義本課題的目的是讓學(xué)生通過(guò)自選一組數(shù)據(jù)，利用所學(xué)專業(yè)知識(shí)在指導(dǎo)教師的指導(dǎo)下獨(dú)立完成對(duì)某一油田或區(qū)塊開發(fā)指標(biāo)的預(yù)測(cè)。本課題要求學(xué)生對(duì)瑞利(Rayleigh)模型的關(guān)系式進(jìn)行推導(dǎo)，其結(jié)果包括，對(duì)油氣田的年產(chǎn)量、累積產(chǎn)量、可采儲(chǔ)量、最高年產(chǎn)量、最高年產(chǎn)量發(fā)生的時(shí)間的預(yù)測(cè)方法。從而將理論知識(shí)和實(shí)際問(wèn)題相結(jié)合，通過(guò)該專題設(shè)計(jì)的訓(xùn)練，加強(qiáng)學(xué)生理論知識(shí)運(yùn)用能力，計(jì)算機(jī)技術(shù)應(yīng)用能力及解決實(shí)際問(wèn)題的工程應(yīng)用

4、能力。1.2 設(shè)計(jì)的主要內(nèi)容根據(jù)已有的基礎(chǔ)數(shù)據(jù)，利用所學(xué)的專業(yè)知識(shí)，在指導(dǎo)教師指導(dǎo)下獨(dú)立完成并提交一個(gè)油田或一個(gè)區(qū)塊油田開發(fā)指標(biāo)預(yù)測(cè)結(jié)果，設(shè)計(jì)主要內(nèi)容如下：(1) 推導(dǎo)瑞利預(yù)測(cè)模型；(2) 根據(jù)油氣田實(shí)際生產(chǎn)數(shù)據(jù)，進(jìn)行線性回歸，求得直線的截距和斜率；(3) 根據(jù)相關(guān)公式，確定瑞利模型常數(shù)a和c；(4) 計(jì)算油田年產(chǎn)油量；(5) 計(jì)算油田累積產(chǎn)油量；(6) 計(jì)算可采儲(chǔ)量；(7) 計(jì)算最高年產(chǎn)量發(fā)生的時(shí)間；(8) 計(jì)算最高年產(chǎn)量；(9) 計(jì)算油田達(dá)到最高產(chǎn)量時(shí)的累積產(chǎn)量；(10) 繪制油田實(shí)際年產(chǎn)量與預(yù)測(cè)產(chǎn)量對(duì)比曲線；(11) 繪制油田實(shí)際累積產(chǎn)量與預(yù)測(cè)累積產(chǎn)量對(duì)比曲線。2 基礎(chǔ)數(shù)據(jù)中國(guó)華北地區(qū)的

5、任丘霧迷山古潛山油田于1975年7月被發(fā)現(xiàn)，1976年投產(chǎn)。其油藏類型為塊狀底水潛山，孔、洞、縫發(fā)育，儲(chǔ)層位于震旦系，以白云巖為主。油藏埋藏深度為25883510m，含油面積為56.9km2，油層有效厚度為272.0m.，有效孔隙度為6%，空氣滲透率為106×10-3m2。原始地層壓力為32.6MPa，原始飽和壓力為1.3MPa，原始?xì)庥捅葹?.4m3/t，地面原油密度為0.889g/cm3，地層油粘度為8.20mPa·s。開始上報(bào)的原油地質(zhì)儲(chǔ)量為53038×104t，可采儲(chǔ)量為11760×104t。油藏采用邊緣底部的注水方式。1990年底綜合含水率為8

6、2%。該油田的開發(fā)數(shù)據(jù)列于下表。 表2-1 開發(fā)數(shù)據(jù)表年份開發(fā)時(shí)間(a)年產(chǎn)量(104t/a)累積產(chǎn)量(104t/a)1975112.2012.2019762606.10618.30197731127.301745.60197841191.502937.10197951240.204177.30198061235.205412.5019817977.956390.4519828901.827292.2719839807.628099.89198410739.388839.27198511675.399514.66198612581.9410096.60198713352.4510449.051

7、98814209.9210658.97198915153.6210812.59199016120.2410932.833 基礎(chǔ)理論瑞利(Rayleigh)模型能夠很好地預(yù)測(cè)對(duì)油氣田的年產(chǎn)量、累積產(chǎn)量、可采儲(chǔ)量、最高年產(chǎn)量、最高年產(chǎn)量發(fā)生的時(shí)間，以及最高年產(chǎn)量對(duì)應(yīng)的累計(jì)產(chǎn)量。3.1瑞利(Rayleigh)模型的應(yīng)用模型的建立在數(shù)理統(tǒng)計(jì)中瑞利分布的分布密度表示為： (3-1)式中 f(x)瑞利分布的分布密度函數(shù)；x分布變量，其值大于等于0；控制分布峰值的參數(shù)。由式(3-l)表示的瑞利分布，在x從0到區(qū)間內(nèi)的累積分布（或分布函數(shù)）等于1，因此可以得到關(guān)系式如下 （3-2）對(duì)于油氣田來(lái)說(shuō)，開發(fā)時(shí)間從0

8、到之間的累計(jì)產(chǎn)量可視為油氣田的可采儲(chǔ)量，可表示為 (3-3) 式中：NR為油氣田的可采儲(chǔ)量，104t（油），108m3（氣）；NP（t）為累積產(chǎn)量，104t（油），108m3（氣）；Q年產(chǎn)量，104t/a（油），108m3/a（氣）；t為開發(fā)時(shí)間，a。將式（3-3）等號(hào)兩端同除以NR得 （3-4）由式（3-2）與式（3-4）對(duì)比可以看出，因此，若將式（3-1）轉(zhuǎn)換為預(yù)測(cè)油氣田產(chǎn)量分布的模型時(shí)，需將該式的右端乘以可采儲(chǔ)量。該可采儲(chǔ)量可視為，將數(shù)理統(tǒng)計(jì)中的理論分布轉(zhuǎn)為預(yù)測(cè)模型的轉(zhuǎn)換常數(shù)。由此，式（3-1）可寫為 （3-5）設(shè)，則；又設(shè)和得 （3-6）再作如下設(shè)定 （3-7）則得到預(yù)測(cè)油氣田產(chǎn)量的模

9、型為 （3-8）油氣田的累積產(chǎn)量表示為 （3-9）將式（3-8）代入式（3-9）進(jìn)行積分得 （3-10）再講式（3-7）代入式（3-10），得到預(yù)測(cè)油氣田累積產(chǎn)量的模型為 （3-11）對(duì)式（3-8）兩邊取導(dǎo)數(shù)，整理得 （3-12）當(dāng)時(shí)，必然有 （3-13）由式（3-13）可以得到最高年產(chǎn)量發(fā)生的時(shí)間為 （3-14）式中：為最高年產(chǎn)量發(fā)生的時(shí)間，a將式（3-14）代入式（3-8）的得到預(yù)測(cè)油氣田最高年產(chǎn)量的公式為 （3-15）式中：為最高年產(chǎn)量，104t/a（油），108m3/a（氣）。再將式（3-14）代入式（3-11）得到預(yù)測(cè)最高年產(chǎn)量對(duì)應(yīng)的累積產(chǎn)量為 （3-16）式中：為最高年產(chǎn)量對(duì)應(yīng)的累

10、積產(chǎn)量，104t/a（油），108m3/a（氣）。將式（3-7）代入式（3-16）得 （3-17）已知，因此式（3-17）又可寫為 （3-18）由式（3-18）可以看出，對(duì)于瑞利模型來(lái)說(shuō)，當(dāng)可采儲(chǔ)量的40%時(shí)，油田進(jìn)入遞減期。 3.2模型的求解方法由上述推導(dǎo)的結(jié)果可以看出，要想利用瑞利模型進(jìn)行各項(xiàng)預(yù)測(cè)，就必須首先確定模型常數(shù)a和c的數(shù)值。為此，將式（3-18）等號(hào)兩端先除以t再取常用對(duì)數(shù)，得 （3-19）若設(shè)： （3-20） （3-21）則得： （3-22）由式（3-22）看出，Q/t與之間具有半對(duì)數(shù)直線關(guān)系。當(dāng)根據(jù)油氣田的實(shí)際生產(chǎn)數(shù)據(jù)，由式（3-22）進(jìn)行線性回歸之后，可以求得直線的截距和斜

11、率的數(shù)值，并由式（3-22）和式（21）確定模型常數(shù)a和c的數(shù)值。 4 設(shè)計(jì)結(jié)果處理表1中的數(shù)據(jù)得到和相應(yīng)的，按照（3-22）式的直線關(guān)系繪于圖4-1，得到了一條很好的直線。經(jīng)線性回歸求得直線的截距=2.5816；直線的斜率=0.007528；直線的相關(guān)系數(shù)r=0.9945。將和的數(shù)值分別代入式（3-20）和式（3-21）求得a為387.59，c為57.76。圖4-1 該油田與的半對(duì)數(shù)關(guān)系將a和c的數(shù)值代入式（3-7）得到該油田的可采儲(chǔ)量為11020×104t；將c的值代入式（3-14）得到最高年產(chǎn)量發(fā)生的時(shí)間為5.3692a；將a和c的數(shù)值代入式（3-15）得到該油田的最高年產(chǎn)量為

12、1242.8×104t；將該油田可采儲(chǔ)量的數(shù)值代入式（3-18）得到該當(dāng)油田年產(chǎn)量達(dá)到最高年產(chǎn)量時(shí)的累積產(chǎn)量為4401.14×104t。將a和c的數(shù)值代入式（3-8）得，預(yù)測(cè)該油田年產(chǎn)量的公式為 （3-23）將a和c的數(shù)值代入式（3-10）得到預(yù)測(cè)該油田累積產(chǎn)量的公式為 （3-24）由式（3-23）和式（3-24）預(yù)測(cè)該油田的產(chǎn)量和累積產(chǎn)量如表4-1表4-1 開發(fā)數(shù)據(jù)表年份開發(fā)時(shí)間(a)預(yù)測(cè)年產(chǎn)量(104t/a)預(yù)測(cè)累積產(chǎn)量(104t/a)19751352.6832177.712919762673.3874694.654619773934.91751504.67197841

13、118.6542538.444197951216.6253713.276198061231.564943.923198171175.1326152.55198281064.9517276.2519839921.08638271.242198410762.85629113.677198511606.43899797.615198612463.547410331.12198713341.14710731.59198814241.979411021.2198915165.56211223.16199016109.33911359.07根據(jù)表4-1中數(shù)據(jù)，在同一坐標(biāo)系中繪制實(shí)際產(chǎn)量與預(yù)測(cè)產(chǎn)量及其對(duì)應(yīng)時(shí)

14、間t之間的關(guān)系曲線，繪于圖4-2。 圖4-2 實(shí)際與預(yù)測(cè)產(chǎn)量隨時(shí)間變化曲線在同一坐標(biāo)系中繪制實(shí)際累計(jì)產(chǎn)量與預(yù)測(cè)累計(jì)產(chǎn)量及其對(duì)應(yīng)時(shí)間t之間的關(guān)系曲線，繪于圖4-3。圖4-3 實(shí)際與預(yù)測(cè)累計(jì)產(chǎn)量隨時(shí)間變化曲線 認(rèn)識(shí)與結(jié)論通過(guò)本次課程設(shè)計(jì)，得出如下認(rèn)識(shí)與結(jié)論：（1）通過(guò)此次課程設(shè)計(jì)，基本掌握了瑞利(Rayleigh)模型的應(yīng)用方法，并能初步對(duì)油田進(jìn)行預(yù)測(cè)分析。（2）能熟練掌握、運(yùn)用瑞利(Rayleigh)模型原理、推導(dǎo)及公式，可以對(duì)油氣田的年產(chǎn)量、累計(jì)產(chǎn)量、可采儲(chǔ)量、最高年產(chǎn)量、最高年產(chǎn)量發(fā)生的時(shí)間以及與最高年產(chǎn)量對(duì)應(yīng)的累積產(chǎn)量進(jìn)行預(yù)測(cè)。（3）通過(guò)計(jì)算機(jī)的編程學(xué)習(xí)與操作，使我對(duì)VB語(yǔ)言有了更進(jìn)一步的

15、認(rèn)識(shí)和了解。同時(shí)，提高了使用計(jì)算機(jī)技應(yīng)用的能力以及解決實(shí)際問(wèn)題的工程應(yīng)用能力，使理論與實(shí)踐相結(jié)合。（4）這次課程設(shè)計(jì)，使我認(rèn)識(shí)到學(xué)好計(jì)算機(jī)的重要性與實(shí)踐操作的必要性，在老師的指導(dǎo)下和與同學(xué)的討論中，不斷地發(fā)現(xiàn)問(wèn)題、不斷地解決問(wèn)題、不斷地領(lǐng)悟、不斷地獲取，使自己的程序不斷完善。（5）通過(guò)這次的課程設(shè)計(jì)，培養(yǎng)了我獨(dú)立思考和動(dòng)手操作的能力，也使我在日常的學(xué)習(xí)與討論當(dāng)中增長(zhǎng)了許多課外知識(shí)，認(rèn)識(shí)到數(shù)學(xué)模型對(duì)油氣田的產(chǎn)量預(yù)測(cè)的重要性。雖然在課程設(shè)計(jì)的過(guò)程中遇到了一些問(wèn)題，但經(jīng)過(guò)反復(fù)的思考與檢查找出原因所在并進(jìn)行改正。參考文獻(xiàn)1陳元千. 瑞利模型的完善推導(dǎo)與應(yīng)用J. 油氣地質(zhì)與采收率，2004，11（4）：

16、39-41.2陳元千 油藏工程實(shí)踐M.石油工業(yè)出版社,2003附 錄 程 序Private Sub Command2_Click()Dim a() As StringDim b() As StringDim c() As StringDim d() As StringDim i As Integeri = 1Open "d:aa.txt" For Input As #1On Error GoTo hDo While Not EOF(1)ReDim Preserve a(i)ReDim Preserve b(i)ReDim Preserve c(i)ReDim Preserv

17、e d(i)Input #1, a(i), b(i), c(i), d(i)i = i + 1Looph:Close #1With MSFlexGrid1.Rows = UBound(a)For i = 1 To UBound(a)MSFlexGrid1.TextMatrix(i - 1, 0) = a(i)MSFlexGrid1.TextMatrix(i - 1, 1) = b(i)MSFlexGrid1.TextMatrix(i - 1, 2) = c(i)MSFlexGrid1.TextMatrix(i - 1, 3) = d(i)Next iEnd WithEnd SubPrivate

18、 Sub Command1_Click()Open "d:aa.txt" For Output As #1With MSFlexGrid1k = .RowsFor i = 0 To k - 1strTemp = ""For j = 0 To .Cols - 1DoEventsstrTemp = strTemp & .TextMatrix(i, j) & ","Next jPrint #1, Left(strTemp, Len(strTemp) - 1)Next iEnd WithClose #1End SubPriva

19、te Sub Command3_Click()Dim a() As StringDim b() As StringDim c() As StringDim d() As StringDim e() As StringDim f() As StringDim i As Integeri = 1Open "d:aa.txt" For Input As #1On Error GoTo hDo While Not EOF(1)ReDim Preserve a(i)ReDim Preserve b(i)ReDim Preserve c(i)ReDim Preserve d(i)Inp

20、ut #1, a(i), b(i), c(i), d(i)i = i + 1Looph:Close #1n = UBound(a)ReDim e(1 To n) As StringReDim f(1 To n) As StringFor i = 2 To UBound(a)e(i) = Str(Log(Val(c(i) / Val(b(i) / Log(Val(10)f(i) = Str(Val(b(i) 2)Next iWith MSFlexGrid3MSFlexGrid3.TextMatrix(0, 0) = "年份"MSFlexGrid3.TextMatrix(0,

21、1) = "LogQ/t"MSFlexGrid3.TextMatrix(0, 2) = "t2".Rows = UBound(a)For i = 2 To UBound(a)MSFlexGrid3.TextMatrix(i - 1, 0) = a(i)MSFlexGrid3.TextMatrix(i - 1, 1) = e(i)MSFlexGrid3.TextMatrix(i - 1, 2) = f(i)Next iEnd WithEnd SubPrivate Sub Command4_Click()k = 10 Val(Text2.Text)l = 1

22、 / (2.303 * Val(Text3.Text)Dim a() As StringDim b() As StringDim c() As StringDim d() As StringDim g() As SingleDim i As Integeri = 1Open "d:aa.txt" For Input As #1On Error GoTo hDo While Not EOF(1)ReDim Preserve a(i)ReDim Preserve b(i)ReDim Preserve c(i)ReDim Preserve d(i)Input #1, a(i),

23、b(i), c(i), d(i)i = i + 1Looph:Close #1n = UBound(a)ReDim g(2 To n + 2) As SingleFor i = 2 To ng(i) = k * Val(b(i) * Exp(-(Val(b(i) 2 / l)Next iWith MSFlexGrid4.Rows = UBound(a)For i = 0 To n - 2MSFlexGrid4.TextMatrix(i + 1, 4) = g(i + 2)Next iEnd WithReDim MyData(n - 2, 3) As Double'-x軸坐標(biāo)值-Y軸坐標(biāo)

24、值-For i = 0 To 15MyData(i, 2) = b(i + 2): MyData(i, 3) = g(i + 2) '本句代表了：第一點(diǎn)數(shù)據(jù)的X軸坐標(biāo)，Y軸坐標(biāo)Next iFor i = 0 To 15MyData(i, 0) = b(i + 2): MyData(i, 1) = c(i + 2)Next i'波形圖外觀設(shè)置With MSChart3.TitleText = "實(shí)際與預(yù)測(cè)年產(chǎn)量隨時(shí)間變化曲線"' '設(shè)置圖線的外觀.Plot.SeriesCollection(1).Pen.Width = 30.Plot.Seri

25、esCollection(1).Pen.Style = VtPenStyleSolid' '設(shè)置XY軸.Plot.Axis(VtChAxisIdX).ValueScale.Auto = False.Plot.Axis(VtChAxisIdY).ValueScale.Auto = False'/ 設(shè)置最大值.Plot.Axis(VtChAxisIdX).ValueScale.Maximum = 18.Plot.Axis(VtChAxisIdY).ValueScale.Maximum = 2000'/ 設(shè)置最小值.Plot.Axis(VtChAxisIdY).Val

26、ueScale.Minimum = 0.Plot.Axis(VtChAxisIdX).ValueScale.Minimum = 0'/.Plot.Axis(VtChAxisIdX).ValueScale.MajorDivision = 12 'X軸主要網(wǎng)格數(shù)量.Plot.Axis(VtChAxisIdY).ValueScale.MajorDivision = 10 'Y軸主要網(wǎng)格數(shù)量.Plot.Axis(VtChAxisIdX).ValueScale.MinorDivision = 0 'X軸次要網(wǎng)格數(shù)量.Plot.Axis(VtChAxisIdY).Value

27、Scale.MinorDivision = 0 'Y軸次要網(wǎng)格數(shù)量.Plot.Axis(VtChAxisIdX).AxisGrid.MajorPen.Style = VtPenStyleDotted.Plot.Axis(VtChAxisIdY).AxisGrid.MajorPen.Style = VtPenStyleDottedMSChart3.Plot.AutoLayout = FalseMSChart3.Plot.UniformAxis = FalseMSChart3.chartType = VtChChartType2dXY '設(shè)置圖形為二維散點(diǎn)圖MSChart3.Cha

28、rtData = MyData '數(shù)據(jù)Dim Index As Integer Index = 1 '設(shè)置圖標(biāo)類型為非自動(dòng)選擇，這樣可以設(shè)置為自己需要的類型MSChart3.Plot.SeriesCollection(Index).ShowLine = False '去除連線MSChart3.Plot.SeriesCollection.Item(Index).SeriesMarker.Auto = False With MSChart3.Plot.SeriesCollection.Item(Index).DataPoints.Item(-1).Marker .Visib

29、le = True '散點(diǎn)的點(diǎn)是否可見 .Size = 80 '散點(diǎn)的大小 .Style = VtMarkerStyleFilledDiamond '選擇菱形的 .FillColor.Automatic = False '顏色不自動(dòng)填充 .FillColor.Set 400, 40, 400 '設(shè)置為黑色 .Pen.VtColor.Set 255, 40, 400 '同樣這個(gè)也必須設(shè)置，否則中間有默認(rèn)的顏色 End WithEnd WithEnd SubPrivate Sub Command5_Click()k = 10 Val(Text2.Te

30、xt)l = 1 / (2.303 * Val(Text3.Text)Dim a() As StringDim b() As StringDim c() As StringDim d() As StringDim h() As SingleDim i As Integeri = 1Open "d:aa.txt" For Input As #1On Error GoTo hDo While Not EOF(1)ReDim Preserve a(i)ReDim Preserve b(i)ReDim Preserve c(i)ReDim Preserve d(i)Input #1

31、, a(i), b(i), c(i), d(i)i = i + 1Looph:Close #1n = UBound(a)ReDim h(2 To n + 2) As SingleFor i = 2 To nh(i) = k * l / 2 * (1 - Exp(-(Val(b(i) 2 / l)Next iWith MSFlexGrid4.Rows = UBound(a)For i = 0 To n - 2MSFlexGrid4.TextMatrix(i + 1, 5) = h(i + 2)Next iEnd WithReDim MyData(n - 2, 3) As Double'-

32、x軸坐標(biāo)值-Y軸坐標(biāo)值-For i = 0 To 15MyData(i, 2) = b(i + 2): MyData(i, 3) = h(i + 2) '本句代表了：第一點(diǎn)數(shù)據(jù)的X軸坐標(biāo)，Y軸坐標(biāo)Next iFor i = 0 To 15MyData(i, 0) = b(i + 2): MyData(i, 1) = d(i + 2)Next i'波形圖外觀設(shè)置With MSChart2.TitleText = "實(shí)際與預(yù)測(cè)累計(jì)產(chǎn)量隨時(shí)間變化曲線"' '設(shè)置圖線的外觀.Plot.SeriesCollection(1).Pen.Width = 3

33、0.Plot.SeriesCollection(1).Pen.Style = VtPenStyleSolid' '設(shè)置XY軸.Plot.Axis(VtChAxisIdX).ValueScale.Auto = False.Plot.Axis(VtChAxisIdY).ValueScale.Auto = False'/ 設(shè)置最大值.Plot.Axis(VtChAxisIdX).ValueScale.Maximum = 18.Plot.Axis(VtChAxisIdY).ValueScale.Maximum = 12000'/ 設(shè)置最小值.Plot.Axis(VtCh

34、AxisIdY).ValueScale.Minimum = 0.Plot.Axis(VtChAxisIdX).ValueScale.Minimum = 0'/.Plot.Axis(VtChAxisIdX).ValueScale.MajorDivision = 12 'X軸主要網(wǎng)格數(shù)量.Plot.Axis(VtChAxisIdY).ValueScale.MajorDivision = 10 'Y軸主要網(wǎng)格數(shù)量.Plot.Axis(VtChAxisIdX).ValueScale.MinorDivision = 0 'X軸次要網(wǎng)格數(shù)量.Plot.Axis(VtChAx

35、isIdY).ValueScale.MinorDivision = 0 'Y軸次要網(wǎng)格數(shù)量.Plot.Axis(VtChAxisIdX).AxisGrid.MajorPen.Style = VtPenStyleDotted.Plot.Axis(VtChAxisIdY).AxisGrid.MajorPen.Style = VtPenStyleDottedMSChart2.Plot.AutoLayout = FalseMSChart2.Plot.UniformAxis = FalseMSChart2.chartType = VtChChartType2dXY '設(shè)置圖形為二維散點(diǎn)圖

36、MSChart2.ChartData = MyData '數(shù)據(jù)Dim Index As Integer Index = 1 '設(shè)置圖標(biāo)類型為非自動(dòng)選擇，這樣可以設(shè)置為自己需要的類型MSChart2.Plot.SeriesCollection(Index).ShowLine = False '去除連線MSChart2.Plot.SeriesCollection.Item(Index).SeriesMarker.Auto = False With MSChart2.Plot.SeriesCollection.Item(Index).DataPoints.Item(-1).M

37、arker .Visible = True '散點(diǎn)的點(diǎn)是否可見 .Size = 80 '散點(diǎn)的大小 .Style = VtMarkerStyleFilledDiamond '選擇菱形的 .FillColor.Automatic = False '顏色不自動(dòng)填充 .FillColor.Set 400, 40, 400 '設(shè)置為黑色 .Pen.VtColor.Set 55, 40, 400 '同樣這個(gè)也必須設(shè)置，否則中間有默認(rèn)的顏色 End WithEnd WithEnd SubPrivate Sub Command6_Click()Dim a()

38、As StringDim b() As StringDim c() As StringDim d() As StringDim e() As StringDim f() As StringDim i As Integeri = 1Open "d:aa.txt" For Input As #1On Error GoTo hDo While Not EOF(1)ReDim Preserve a(i)ReDim Preserve b(i)ReDim Preserve c(i)ReDim Preserve d(i)Input #1, a(i), b(i), c(i), d(i)i

39、= i + 1Looph:Close #1n = UBound(a)ReDim e(2 To n) As StringReDim f(2 To n) As StringFor i = 2 To UBound(a)e(i) = Str(Log(Val(c(i) / Val(b(i) / Log(Val(10)f(i) = Str(Val(b(i) 2)Next iReDim MyData(n - 2, 1) As Double '-x軸坐標(biāo)值-Y軸坐標(biāo)值- For i = 0 To UBound(a) - 2 MyData(i, 0) = f(i + 2): MyData(i, 1) =

40、 e(i + 2) '本句代表了：第一點(diǎn)數(shù)據(jù)的X軸坐標(biāo),Y軸坐標(biāo) Next i'波形圖外觀設(shè)置With MSChart1 .TitleText = "Rayleigh曲線"' '設(shè)置圖線的外觀 .Plot.SeriesCollection(1).Pen.Width = 30 .Plot.SeriesCollection(1).Pen.Style = VtPenStyleSolid' '設(shè)置XY軸 .Plot.Axis(VtChAxisIdX).ValueScale.Auto = False .Plot.Axis(VtChAxi

41、sIdY).ValueScale.Auto = False '/ 設(shè)置最大值 .Plot.Axis(VtChAxisIdX).ValueScale.Maximum = 300 .Plot.Axis(VtChAxisIdY).ValueScale.Maximum = 3 '/ 設(shè)置最小值 .Plot.Axis(VtChAxisIdY).ValueScale.Minimum = 0 .Plot.Axis(VtChAxisIdX).ValueScale.Minimum = 0 '/ .Plot.Axis(VtChAxisIdX).ValueScale.MajorDivisio

42、n = 20 'X軸主要網(wǎng)格數(shù)量 .Plot.Axis(VtChAxisIdY).ValueScale.MajorDivision = 20 'Y軸主要網(wǎng)格數(shù)量 .Plot.Axis(VtChAxisIdX).ValueScale.MinorDivision = 0 'X軸次要網(wǎng)格數(shù)量 .Plot.Axis(VtChAxisIdY).ValueScale.MinorDivision = 0 'Y軸次要網(wǎng)格數(shù)量 .Plot.Axis(VtChAxisIdX).AxisGrid.MajorPen.Style = VtPenStyleDotted .Plot.Axis

43、(VtChAxisIdY).AxisGrid.MajorPen.Style = VtPenStyleDotted MSChart1.Plot.AutoLayout = False MSChart1.Plot.UniformAxis = False MSChart1.chartType = VtChChartType2dXY '設(shè)置圖形為二維散點(diǎn)圖 MSChart1.ChartData = MyData '數(shù)據(jù) Dim Index As Integer Index = 1 '設(shè)置圖標(biāo)類型為非自動(dòng)選擇，這樣可以設(shè)置為自己需要的類型MSChart1.Plot.SeriesCo

44、llection(Index).ShowLine = False '去除連線MSChart1.Plot.SeriesCollection.Item(Index).SeriesMarker.Auto = False With MSChart1.Plot.SeriesCollection.Item(Index).DataPoints.Item(-1).Marker .Visible = True '散點(diǎn)的點(diǎn)是否可見 .Size = 80 '散點(diǎn)的大小 .Style = VtMarkerStyleFilledDiamond '選擇菱形的 .FillColor.Auto

45、matic = False '顏色不自動(dòng)填充 .FillColor.Set 400, 40, 400 '設(shè)置為黑色 .Pen.VtColor.Set 255, 40, 400 '同樣這個(gè)也必須設(shè)置，否則中間有默認(rèn)的顏色 End WithEnd WithEnd SubPrivate Sub Command7_Click()Dim a() As StringDim b() As StringDim c() As StringDim d() As StringDim e() As StringDim f() As StringDim i As Integeri = 1Open

46、 "d:aa.txt" For Input As #1On Error GoTo hDo While Not EOF(1)ReDim Preserve a(i)ReDim Preserve b(i)ReDim Preserve c(i)ReDim Preserve d(i)Input #1, a(i), b(i), c(i), d(i)i = i + 1Looph:Close #1n = UBound(a)ReDim e(2 To n) As StringReDim f(2 To n) As StringFor i = 2 To UBound(a)e(i) = Str(Lo

47、g(Val(c(i) / Val(b(i) / Log(Val(10)f(i) = Str(Val(b(i) 2)Next iReDim MyData(n - 2, 3) As Double '-x軸坐標(biāo)值-Y軸坐標(biāo)值- For i = 0 To UBound(a) - 2 MyData(i, 0) = f(i + 2): MyData(i, 1) = e(i + 2) '本句代表了：第一點(diǎn)數(shù)據(jù)的X軸坐標(biāo)，Y軸坐標(biāo) Next i MyData(0, 2) = Text4.Text: MyData(0, 3) = Text5.Text For i = 1 To UBound(a)

48、 - 2 MyData(i, 2) = Text6.Text: MyData(i, 3) = Text7.Text Next i '波形圖外觀設(shè)置With MSChart1 .TitleText = "Rayleigh曲線"' '設(shè)置圖線的外觀 .Plot.SeriesCollection(1).Pen.Width = 30 .Plot.SeriesCollection(1).Pen.Style = VtPenStyleSolid' '設(shè)置XY軸 .Plot.Axis(VtChAxisIdX).ValueScale.Auto = Fa

49、lse .Plot.Axis(VtChAxisIdY).ValueScale.Auto = False '/ 設(shè)置最大值 .Plot.Axis(VtChAxisIdX).ValueScale.Maximum = 300 .Plot.Axis(VtChAxisIdY).ValueScale.Maximum = 3 '/ 設(shè)置最小值 .Plot.Axis(VtChAxisIdY).ValueScale.Minimum = 0 .Plot.Axis(VtChAxisIdX).ValueScale.Minimum = 0 '/ .Plot.Axis(VtChAxisIdX).V

50、alueScale.MajorDivision = 20 'X軸主要網(wǎng)格數(shù)量 .Plot.Axis(VtChAxisIdY).ValueScale.MajorDivision = 20 'Y軸主要網(wǎng)格數(shù)量 .Plot.Axis(VtChAxisIdX).ValueScale.MinorDivision = 0 'X軸次要網(wǎng)格數(shù)量 .Plot.Axis(VtChAxisIdY).ValueScale.MinorDivision = 0 'Y軸次要網(wǎng)格數(shù)量 .Plot.Axis(VtChAxisIdX).AxisGrid.MajorPen.Style = VtPen

51、StyleDotted .Plot.Axis(VtChAxisIdY).AxisGrid.MajorPen.Style = VtPenStyleDotted MSChart1.Plot.AutoLayout = False MSChart1.Plot.UniformAxis = False MSChart1.chartType = VtChChartType2dXY '設(shè)置圖形為二維散點(diǎn)圖 MSChart1.ChartData = MyData '數(shù)據(jù) .Plot.Axis(VtChAxisIdX, 0).AxisTitle = "t2" .Plot.Axi

52、s(VtChAxisIdY, 0).AxisTitle = "LogQ/t" Dim Index As Integer Index = 1 '設(shè)置圖標(biāo)類型為非自動(dòng)選擇，這樣可以設(shè)置為自己需要的類型MSChart1.Plot.SeriesCollection(Index).ShowLine = False '去除連線MSChart1.Plot.SeriesCollection.Item(Index).SeriesMarker.Auto = False With MSChart1.Plot.SeriesCollection.Item(Index).DataPoints.Item(-1).Marker .Visible = True '散點(diǎn)的點(diǎn)是否可見 .Size = 80 '散點(diǎn)的大小 .Style = VtMarkerStyle