1、東 北 石 油 大 學課 程 設 計課 程 石油工程課程設計 題 目 瑞利(Rayleigh)模型的應用 院 系 石油工程學院 專業(yè)班級 油工12-5班 學生姓名 學生學號 120201140111 指導教師 王立軍 2015年7月24日東北石油大學課程設計任務書課程 石油工程課程設計 題目 瑞利(Rayleigh)模型的應用 專業(yè) 石油工程 姓名 XXX 學號 120201140111 主要內(nèi)容、基本要求、主要參考資料等主要內(nèi)容：(1) 推導瑞利預測模型；(2) 根據(jù)油氣田實際生產(chǎn)數(shù)據(jù)，進行線性回歸，求得直線的截距和斜率；(3) 根據(jù)相關公式，確定瑞利模型常數(shù)a和c；(4) 計算油田年產(chǎn)油量

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

3、目 錄1 前言·················································

4、;·······31.1設計的目的意義·········································

5、31.2設計的主要內(nèi)容·············································32 基礎數(shù)據(jù)··

6、··················································

7、43 基礎理論·················································

8、;···53.1 瑞利（Rayleigh）模型的應用模型的建立·······················53.2模型的求解方法·················

9、83;···························74 設計結果·····················&#

10、183;······························8認識與結論··················&#

11、183;·································12參考文獻···············&#

12、183;······································13附錄··········

13、3;···············································141.前 言1.1 設計的目的

14、意義本課題的目的是讓學生通過自選一組數(shù)據(jù)，利用所學專業(yè)知識在指導教師的指導下獨立完成對某一油田或區(qū)塊開發(fā)指標的預測。本課題要求學生對瑞利(Rayleigh)模型的關系式進行推導，其結果包括，對油氣田的年產(chǎn)量、累積產(chǎn)量、可采儲量、最高年產(chǎn)量、最高年產(chǎn)量發(fā)生的時間的預測方法。從而將理論知識和實際問題相結合，通過該專題設計的訓練，加強學生理論知識運用能力，計算機技術應用能力及解決實際問題的工程應用能力。1.2 設計的主要內(nèi)容根據(jù)已有的基礎數(shù)據(jù)，利用所學的專業(yè)知識，在指導教師指導下獨立完成并提交一個油田或一個區(qū)塊油田開發(fā)指標預測結果，設計主要內(nèi)容如下：(1) 推導瑞利預測模型；(2) 根據(jù)油氣田實際生

15、產(chǎn)數(shù)據(jù)，進行線性回歸，求得直線的截距和斜率；(3) 根據(jù)相關公式，確定瑞利模型常數(shù)a和c；(4) 計算油田年產(chǎn)油量；(5) 計算油田累積產(chǎn)油量；(6) 計算可采儲量；(7) 計算最高年產(chǎn)量發(fā)生的時間；(8) 計算最高年產(chǎn)量；(9) 計算油田達到最高產(chǎn)量時的累積產(chǎn)量；(10) 繪制油田實際年產(chǎn)量與預測產(chǎn)量對比曲線；(11) 繪制油田實際累積產(chǎn)量與預測累積產(chǎn)量對比曲線。2.基礎數(shù)據(jù)中國華北地區(qū)的任丘霧迷山古潛山油田于1975年7月被發(fā)現(xiàn)，1976年投產(chǎn)。其油藏類型為塊狀底水潛山，孔、洞、縫發(fā)育，儲層位于震旦系，以白云巖為主。油藏埋藏深度為25883510m，含油面積為56.9km2，油層有效厚度為

16、272.0m.，有效孔隙度為6%，空氣滲透率為106×10-3m2。原始地層壓力為32.6MPa，原始飽和壓力為1.3MPa，原始氣油比為4.4m3/t，地面原油密度為0.889g/cm3，地層油粘度為8.20mPa·s。開始上報的原油地質儲量為53038×104t，可采儲量為11760×104t。油藏采用邊緣底部的注水方式。1990年底綜合含水率為82%。該油田的開發(fā)數(shù)據(jù)列于下表。 表2-1 開發(fā)數(shù)據(jù)表年份開發(fā)時間(a)年產(chǎn)量(104t/a)累積產(chǎn)量(104t/a)1975112.2012.2019762606.10618.30197731127.30

17、1745.60197841191.502937.10197951240.204177.30198061235.205412.5019817977.956390.4519828901.827292.2719839807.628099.89198410739.388839.27198511675.399514.66198612581.9410096.60198713352.4510449.05198814209.9210658.97198915153.6210812.59199016120.2410932.833.基礎理論3.1瑞利(Rayleigh)模型的應用模型的建立在數(shù)理統(tǒng)計中瑞利分布的分布

18、密度表示為： (3-1)式中 f(x)瑞利分布的分布密度函數(shù)；x分布變量，其值大于等于0；控制分布峰值的參數(shù)。由式(3-l)表示的瑞利分布，在x從0到區(qū)間內(nèi)的累積分布（或分布函數(shù)）等于1，因此可以得到關系式如下 （3-2）對于油氣田來說，開發(fā)時間從0到之間的累計產(chǎn)量可視為油氣田的可采儲量，可表示為 (3-3) 式中：NR為油氣田的可采儲量，104t（油），108m3（氣）；NP（t）為累積產(chǎn)量，104t（油），108m3（氣）；Q年產(chǎn)量，104t/a（油），108m3/a（氣）；t為開發(fā)時間，a。將式（3-3）等號兩端同除以NR得 （3-4）由式（3-2）與式（3-4）對比可以看出，因此，若將

19、式（3-1）轉換為預測油氣田產(chǎn)量分布的模型時，需將該式的右端乘以可采儲量。該可采儲量可視為，將數(shù)理統(tǒng)計中的理論分布轉為預測模型的轉換常數(shù)。由此，式（3-1）可寫為 （3-5）設，則；又設和得 （3-6）再作如下設定 （3-7）則得到預測油氣田產(chǎn)量的模型為 （3-8）油氣田的累積產(chǎn)量表示為 （3-9）將式（3-8）代入式（3-9）進行積分得 （3-10）再講式（3-7）代入式（3-10），得到預測油氣田累積產(chǎn)量的模型為 （3-11）對式（3-8）兩邊取導數(shù)，整理得 （3-12）當時，必然有 （3-13）由式（3-13）可以得到最高年產(chǎn)量發(fā)生的時間為 （3-14）式中：為最高年產(chǎn)量發(fā)生的時間，a將

20、式（3-14）代入式（3-8）的得到預測油氣田最高年產(chǎn)量的公式為 （3-15）式中：為最高年產(chǎn)量，104t/a（油），108m3/a（氣）。再將式（3-14）代入式（3-11）得到預測最高年產(chǎn)量對應的累積產(chǎn)量為 （3-16）式中：為最高年產(chǎn)量對應的累積產(chǎn)量，104t/a（油），108m3/a（氣）。將式（3-7）代入式（3-16）得 （3-17）已知，因此式（3-17）又可寫為 （3-18）由式（3-18）可以看出，對于瑞利模型來說，當可采儲量的40%時，油田進入遞減期。 3.2模型的求解方法由上述推導的結果可以看出，要想利用瑞利模型進行各項預測，就必須首先確定模型常數(shù)a和c的數(shù)值。為此，將式

21、（3-18）等號兩端先除以t再取常用對數(shù)，得 （3-19）若設： （3-20） （3-21）則得： （3-22）由式（3-22）看出，Q/t與之間具有半對數(shù)直線關系。當根據(jù)油氣田的實際生產(chǎn)數(shù)據(jù)，由式（3-22）進行線性回歸之后，可以求得直線的截距和斜率的數(shù)值，并由式（3-22）和式（21）確定模型常數(shù)a和c的數(shù)值。4.設計結果本程序主要以下拉菜單控制各個過程為主體，實現(xiàn)各個不同的功能，簡潔明了的展示出各個過程。當啟動程序時，初始界面如圖（3-1）所示： 圖3-1初始運行界面單擊“導入數(shù)據(jù)”，程序將從存放數(shù)據(jù)的默認位置導入初始數(shù)據(jù)。單擊“導出數(shù)據(jù)”，原始數(shù)據(jù)將被保存在指定的默認位置。單擊“保存數(shù)

22、據(jù)”，預測結果將被保存在指定的默認位置。如果想修改數(shù)據(jù)可直接單擊想要修改的數(shù)據(jù)所對應的單元格進行修改，累計產(chǎn)量將在再次單擊參數(shù)計算時進行改正。當單擊“參數(shù)計算”時，程序將算出不同時間所對應的和lg（q/t）,并根據(jù)最小二乘法（詳情見附錄程序），計算出回歸直線的斜率=-0.0067；直線的截距=2.5495，將和的數(shù)值分別代入式（3-20）和式（3-21）求得a為354.4051，c為64.8929。單擊”擬合公式”，回歸直線的公式便可顯示在表格的下方，運行界面如圖（3-2）所示：圖3-2 參數(shù)計算及公式擬合運行界面單擊“瑞麗曲線”按鈕，可將（3-22）式的直線關系繪于圖3-3，得到了一條擬合很

23、好的直線，如圖（3-3）所示： 圖3-3 瑞利曲線運行界面將a和c的數(shù)值代入式（3-7）得到該油田的可采儲量為11499.1931×104t；將c的值代入式（3-14）得到最高年產(chǎn)量發(fā)生的時間為5.6953a；將a和c的數(shù)值代入式（3-15）得到該油田的最高年產(chǎn)量為1224.4885×104t；將該油田可采儲量的數(shù)值代入式（3-18）得到該當油田年產(chǎn)量達到最高年產(chǎn)量時的累積產(chǎn)量為4524.5606×104。由式（3-23）和式（3-24）預測該油田的產(chǎn)量和累積產(chǎn)量，單擊 “預測結果”，可將預測出的各個結果顯示在圖（3-4）的界面上。3-4 預測結果運行界面單擊“實

24、際與預測年產(chǎn)量與時間變化曲線”按鈕，可在同一坐標系中繪制實際產(chǎn)量與預測產(chǎn)量及其對應時間t之間的關系曲線，繪于圖（3-5）。圖3-5 實際與預測產(chǎn)量隨時間變化曲線運行界面單擊“實際與預測累計產(chǎn)量與時間變化曲線”按鈕，可在同一坐標系中繪制實際累計產(chǎn)量與預測累計產(chǎn)量及其對應時間t之間的關系曲線繪于圖3-6。圖3-6實際與預測累計產(chǎn)量隨時間變化曲線運行界面認識與結論經(jīng)過了三周從沒接觸過的課程設計，讓我獲得了很多收獲，和同學們每天準時相約機房，一起探討編程的團結協(xié)作，面對難題，絞盡腦汁，使難題獲得解決時內(nèi)心的滿足，都讓我感覺非常的充實。大家一起經(jīng)歷的三周，有探討時的說笑，有鉆研時的認真，也有請教他人時的

25、謙遜，當然還有老師指導時的孜孜不倦。這次課程設計，我所得到的收獲大致如下：(1)通過課程設計所設計的瑞利模型程序，讓我對vb編程有了更深入的了解，運用企業(yè)版vb，比上課時的基 礎版實現(xiàn)的功能更多樣，也更加強大，經(jīng)過查閱資料，學會了不少新控件的使用方法，和一些新的編程語言，相信在以后的工作中還會更好的得以運用。(2)課程設計雖然結束了，但是我所設計的程序并不完善，還有很多功能有待實現(xiàn)，這就要求自己還要不斷的學習編程知識，將程序加以完善，以便在以后的工作中，加以運用。(3)這次課程設計讓我們養(yǎng)成了遇到問題運用編程加以解決的習慣，這在以后的工作中會變得非常實用，通過編程，讓以后參加工作的我們工作起來

26、更加高效。(4)這次課程設計，感謝張老師對我們編程的指導，已經(jīng)王老師在中期對我的程序的評判，讓我在后面有了最終的方向，對我的課設的完成起了關鍵作用，最后感謝兩位老師在假期的陪伴。參考文獻1陳濤平，石油工程，石油工業(yè)出版社，20112吳雅娟，王莉利，Visual Basic 程序設計案例教程，中國石化出版社，20113蔣加伏、張林峰.Visual Basic程序設計教程M.北京郵電大學出版社,2011。4陳元千. 瑞利模型的完善推導與應用J. 油氣地質與采收率，2004，11（4）：39-41.5陳元千 油藏工程實踐M.石油工業(yè)出版社,2003附 錄一、程 序Dim t2(16) As

27、 Single, lg(16) As Single, j As Single, t(16) As Single, q(16) As Single, ql(16) As Single, X As Single, Z As Single, xx As Single, zx As Single Const ASC_ENTER = 13 'Enter Dim gRow As Integer Dim gCol As IntegerDim nf(16) As Single, kfsj(16) As Single, cl(16) As Single, ljcl(16) As Single '

28、年份，開發(fā)時間，產(chǎn)量，累計產(chǎn)量Private Sub Form_Resize()Me.PaintPicture Me.Picture, 0, 0, Me.ScaleWidth, Me.ScaleHeight '使圖片填充End Sub'初始加載Private Sub Form_Load()Me.AutoRedraw = TrueGrid1.ColWidth(0) = 800 '確定表格寬度Grid1.ColWidth(1) = 1100Grid1.ColWidth(2) = 1400Grid1.ColWidth(3) = 1600Grid1.ColWidth(4) =

29、1500Grid1.ColWidth(5) = 1500Grid1.TextMatrix(0, 4) = "t2"Grid1.TextMatrix(0, 5) = "lg(q/t)"Grid3.TextMatrix(0, 0) = "年份"Grid3.TextMatrix(0, 1) = "開發(fā)時間(a)"Grid3.TextMatrix(0, 2) = "年產(chǎn)量(104 t/a)"Grid3.TextMatrix(0, 3) = "預測年產(chǎn)油量(104 t/a)"Grid3

30、.TextMatrix(0, 4) = "累積產(chǎn)量(104 t/a)"Grid3.TextMatrix(0, 5) = "預測累積產(chǎn)油量(104 t/a)"Grid3.ColWidth(0) = 600 '表格大小Grid3.ColWidth(1) = 1200Grid3.ColWidth(2) = 1500Grid3.ColWidth(3) = 2200Grid3.ColWidth(4) = 1800Grid3.ColWidth(5) = 2500Text1.Visible = False '設置可視性Text2.Visible = F

31、alseText3.Visible = FalseText4.Visible = FalseText5.Visible = FalseText6.Visible = FalseText7.Visible = FalseText8.Visible = FalseText9.Visible = FalseText10.Visible = FalseGrid1.Visible = FalseGrid3.Visible = False Label1.Visible = FalseLabel1.Visible = FalseLabel2.Visible = FalseLabel3.Visible = F

32、alseLabel4.Visible = FalseLabel5.Visible = FalseLabel6.Visible = FalseLabel7.Visible = FalseLabel8.Visible = FalseLabel9.Visible = FalseMSChart1.Visible = FalseMSChart2.Visible = FalseMSChart3.Visible = FalseGrid1.Left = (Form1.Width - Grid1.Width) / 2 - 2000 '各圖標位置加載Grid1.Top = (Form1.Height -

33、Grid1.Height) / 2Label1.Left = (Form1.Width - Grid1.Width) / 2 + 6300Text2.Left = (Form1.Width - Grid1.Width) / 2 + 8000Label2.Left = (Form1.Width - Grid1.Width) / 2 + 6300Label3.Left = (Form1.Width - Grid1.Width) / 2 + 2000Label4.Left = (Form1.Width - Grid1.Width) / 2 + 6300Label5.Left = (Form1.Wid

34、th - Grid1.Width) / 2 + 6300Text3.Left = (Form1.Width - Grid1.Width) / 2 + 8000Text4.Left = (Form1.Width - Grid1.Width) / 2 + 3700Text5.Left = (Form1.Width - Grid1.Width) / 2 + 8720Text6.Left = (Form1.Width - Grid1.Width) / 2 + 9300Label1.Top = (Form1.Height - Grid1.Height) / 2Label2.Top = (Form1.He

35、ight - Grid1.Height) / 2 + 1000Label3.Top = (Form1.Height - Grid1.Height) / 2 + 5000Label4.Top = (Form1.Height - Grid1.Height) / 2 + 2000Label5.Top = (Form1.Height - Grid1.Height) / 2 + 3000Text2.Top = (Form1.Height - Grid1.Height) / 2Text3.Top = (Form1.Height - Grid1.Height) / 2 + 1000Text4.Top = (

36、Form1.Height - Grid1.Height) / 2 + 5000Text5.Top = (Form1.Height - Grid1.Height) / 2 + 2000Text6.Top = (Form1.Height - Grid1.Height) / 2 + 3000Grid3.Left = (Form1.Width - Grid1.Width) / 2 - 1500Grid3.Top = (Form1.Height - Grid1.Height) / 2 - 1500Label6.Top = 5500Label9.Top = 5500Text7.Top = 5500Text

37、10.Top = 5500Label7.Top = 6200Label8.Top = 6200Text8.Top = 6200Text9.Top = 6200Label6.Left = 600Label9.Left = 7080Label7.Left = 600Label8.Left = 7080Text7.Left = 3480Text10.Left = 9960Text8.Left = 3480Text9.Left = 9960MSChart1.Left = (Form1.Width - MSChart1.Width) / 2MSChart1.Top = (Form1.Width - MS

38、Chart1.Height) / 2 - 3000MSChart2.Left = (Form1.Width - MSChart2.Width) / 2MSChart2.Top = (Form1.Width - MSChart2.Height) / 2 - 3000MSChart3.Left = (Form1.Width - MSChart3.Width) / 2MSChart3.Top = (Form1.Width - MSChart3.Height) / 2 - 3000save.Enabled = False '按鍵是否可用csjs.Enabled = Falsenihegongs

39、hi.Enabled = Falseyucejieguo.Enabled = Falserayleigh.Enabled = Falsenianchanliang.Enabled = Falseleijichanliang.Enabled = Falsedaochushuju.Enabled = FalseEnd SubPrivate Sub daorushuju_Click()Label10.Visible = Falsedaochushuju.Enabled = Truecsjs.Enabled = TrueGrid1.Visible = TrueDim i As Single, k As

40、 SingleOpen "C:Documents and Settingsxx桌面基礎數(shù)據(jù).txt" For Input As #1i = 0k = Grid1.RowsWhile Not EOF(1)Input #1, nf(i), kfsj(i), cl(i), ljcl(i)Grid1.TextMatrix(i, 0) = nf(i)Grid1.TextMatrix(i, 1) = kfsj(i)Grid1.TextMatrix(i, 2) = cl(i)Grid1.TextMatrix(i, 3) = ljcl(i)i = i + 1 If i > Grid1

41、.Rows - 1 Then Grid1.Rows = i + 1WendIf k <> Grid1.Rows Then Grid1.Rows = Grid1.Rows - 1Grid1.TextMatrix(0, 0) = "年份"Grid1.TextMatrix(0, 1) = "開發(fā)時間(a)"Grid1.TextMatrix(0, 2) = "年產(chǎn)量104t/a"Grid1.TextMatrix(0, 3) = "累計產(chǎn)量104t/a"Close #1End SubPrivate Sub dao

42、chushuju_Click()If MsgBox("是否導出", vbYesNo) = vbYes ThenOpen "C:Documents and Settingsxx桌面導出數(shù)據(jù).txt" For Output As #1 For i = 1 To Grid1.Rows - 1 Print #1, Grid1.TextMatrix(i, 0), Grid1.TextMatrix(i, 1), Grid1.TextMatrix(i, 2), Grid1.TextMatrix(i, 3), Grid1.TextMatrix(i, 4) '將數(shù)

43、據(jù)保存到文本文檔中 Next iClose #1MsgBox ("導出成功")End IfEnd SubPrivate Sub save_Click()If MsgBox("是否保存", vbYesNo) = vbYes Then Open "C:Documents and Settingsxx桌面數(shù)據(jù)保存.txt" For Output As #1For i = 1 To Grid3.Rows - 1Print #1, Grid3.TextMatrix(i, 0), Grid3.TextMatrix(i, 1), Grid3.Tex

44、tMatrix(i, 2), Grid3.TextMatrix(i, 3), Grid3.TextMatrix(i, 4), Grid3.TextMatrix(i, 5)Next iClose #1MsgBox ("保存成功")End IfEnd SubPrivate Sub csjs_Click()nihegongshi.Enabled = Trueyucejieguo.Enabled = Truerayleigh.Enabled = TrueLabel1.Visible = TrueLabel2.Visible = TrueLabel5.Visible = TrueLa

45、bel4.Visible = TrueText2.Visible = TrueText3.Visible = TrueText5.Visible = TrueText6.Visible = TrueMSChart1.Visible = FalseMSChart2.Visible = FalseMSChart3.Visible = FalseDim n!n = 0For j = 1 To Grid1.Rows - 1 n = Grid1.TextMatrix(j, 2) + n Grid1.TextMatrix(j, 3) = nNext jFor j = 1 To Grid1.Rows - 1

46、 '最小二乘法 t(j) = Grid1.TextMatrix(j, 1) q(j) = Grid1.TextMatrix(j, 2) t2(j) = (t(j) 2 Grid1.TextMatrix(j, 4) = t2(j) lg(j) = Log(q(j) / t(j) / Log(10) Grid1.TextMatrix(j, 5) = lg(j) If j > Grid1.Rows - 1 Then Grid1.Rows = j + 1Next jFor j = 2 To Grid1.Rows - 1 X = X + t(j) 2 Z = Z + Log(q(j) /

47、t(j) / Log(10) zx = zx + Log(q(j) / t(j) / Log(10) * t(j) 2 xx = xx + (t(j) 2) 2Next jText2 = Int(Grid1.Rows - 2) * zx - Z * X) / (Grid1.Rows - 2) * xx - X 2) * 10000) / 10000Text3 = Int(Z - Text2 * X) / (Grid1.Rows - 2) * 10000) / 10000Text5.Text = Format(10 (Text3.Text), "0.0000")Text6.T

48、ext = Format(1 / (-Form1.Text2.Text) / 2.3, "0.0000")End SubPrivate Sub nihegongshi_Click()Text4.Visible = TrueLabel3.Visible = TrueText4.Text = "lg(q/t)=" & Text2 & "t2" & "+" & Text3End SubPrivate Sub yucejieguo_Click()save.Enabled = Truenian

49、chanliang.Enabled = Trueleijichanliang.Enabled = TrueGrid1.Visible = FalseGrid3.Visible = TrueLabel6.Visible = TrueLabel7.Visible = TrueLabel8.Visible = TrueLabel9.Visible = TrueText7.Visible = TrueText8.Visible = TrueText9.Visible = TrueText10.Visible = TrueLabel1.Visible = FalseLabel2.Visible = Fa

50、lseLabel3.Visible = FalseLabel5.Visible = FalseLabel4.Visible = FalseText2.Visible = FalseText3.Visible = FalseText4.Visible = FalseText5.Visible = FalseText6.Visible = FalseMSChart1.Visible = FalseMSChart2.Visible = FalseMSChart3.Visible = FalseDim g!, a!, c!Dim j(1 To 17) As StringDim h(1 To 17) A

51、s StringDim l(1 To 17) As StringFor i = 1 To Grid3.Rows - 1 Grid3.TextMatrix(i, 0) = Grid1.TextMatrix(i, 0) Grid3.TextMatrix(i, 1) = Grid1.TextMatrix(i, 1) Grid3.TextMatrix(i, 2) = Grid1.TextMatrix(i, 2) Grid3.TextMatrix(i, 3) = Grid1.TextMatrix(i, 3) Grid3.TextMatrix(i, 4) = Grid1.TextMatrix(i, 4)

52、Grid3.TextMatrix(i, 5) = Grid1.TextMatrix(i, 5)Next ig = 0For i = 1 To Grid3.Rows - 1l(i) = Val(Text2.Text) * t(i) 2 + Val(Text3.Text)h(i) = t(i) * 10 l(i)g = g + h(i)j(i) = gGrid3.TextMatrix(i, 0) = Grid1.TextMatrix(i, 0) Grid3.TextMatrix(i, 2) = Grid1.TextMatrix(i, 2) Grid3.TextMatrix(i, 1) = Grid

53、1.TextMatrix(i, 1) Grid3.TextMatrix(i, 3) = Int(h(i) * 10000) / 10000 Grid3.TextMatrix(i, 4) = Grid1.TextMatrix(i, 3)Grid3.TextMatrix(i, 5) = Int(j(i) * 10000) / 10000 Next i a = 10 (Text3.Text) c = 1 / Val(-Text2.Text) / 2.3 Text7.Text = Int(0.5 * a * c * 10000) / 10000 Text8.Text = Int(0.4289 * a

54、* c 0.5 * 10000) / 10000 Text9.Text = Int(0.707 * (c 0.5 * 10000) / 10000 Text10.Text = Int(0.5 * a * c * (1 - 10 (-0.5) / 0.52137) * 10000) / 10000End SubPrivate Sub rayleigh_Click() Label1.Visible = FalseLabel1.Visible = FalseLabel2.Visible = FalseLabel3.Visible = FalseLabel4.Visible = FalseLabel5

55、.Visible = FalseLabel6.Visible = FalseLabel7.Visible = FalseLabel8.Visible = FalseLabel9.Visible = FalseText1.Visible = FalseText2.Visible = FalseText3.Visible = FalseText4.Visible = FalseText5.Visible = FalseText6.Visible = FalseText7.Visible = FalseText8.Visible = FalseText9.Visible = FalseText10.

56、Visible = FalseGrid1.Visible = FalseGrid3.Visible = FalseMSChart2.Visible = FalseMSChart3.Visible = FalseMSChart1.Visible = TrueDim n!, i!, MyData#(), t!(1 To 16), j%For j = 1 To Grid1.Rows - 1 t(j) = Grid1.TextMatrix(j, 1)Next j n = Grid1.Rows - 1ReDim r(1 To n, 1 To n) As StringFor i = 1 To Grid1.

57、Rows - 1r(i, 2) = Log(q(i) / t(i) / Log(10)r(i, 1) = t(i) 2Next iReDim MyData(n - 1, 3) As Double '-x軸坐標值-Y軸坐標值- For i = 1 To Grid1.Rows - 1 '第一點數(shù)據(jù)的X軸坐標,Y軸坐標 MyData(i - 1, 0) = r(i, 1): MyData(i - 1, 1) = r(i, 2) MyData(i - 1, 3) = Val(Text2.Text) * r(i, 1) + Val(Text3.Text) MyData(i - 1, 2)

58、 = r(i, 1) Next i '波形圖外觀設置With MSChart1 .Plot.SeriesCollection(1).Pen.Width = 30 .Plot.SeriesCollection(1).Pen.Style = VtPenStyleSolid .Plot.Axis(VtChAxisIdX).ValueScale.Auto = False .Plot.Axis(VtChAxisIdY).ValueScale.Auto = False .Plot.Axis(VtChAxisIdX).ValueScale.Maximum = 300 .Plot.Axis(VtChA

59、xisIdY).ValueScale.Maximum = 3 .Plot.Axis(VtChAxisIdY).ValueScale.Minimum = 0 .Plot.Axis(VtChAxisIdX).ValueScale.Minimum = 0 .Plot.Axis(VtChAxisIdX).ValueScale.MajorDivision = 12 .Plot.Axis(VtChAxisIdY).ValueScale.MajorDivision = 10 .Plot.Axis(VtChAxisIdX).ValueScale.MinorDivision = 0 .Plot.Axis(VtChAxisIdY).ValueScale.MinorDivision = 0 .Plot.Axis(VtChAxisIdX).AxisGrid.MajorPen.Style = VtPenStyle