石油工程學(xué)院課程設(shè)計(jì)（報(bào)告）PAGE東北石油大學(xué)課程設(shè)計(jì)課程課程石油工程課程設(shè)計(jì)題目瑞利(Rayleigh)模型的應(yīng)用院系石油工程學(xué)院專業(yè)班級(jí)油工班學(xué)生姓名學(xué)生學(xué)號(hào)指導(dǎo)教師王立軍20年7月19日PAGEPAGE17東北石油大學(xué)課程設(shè)計(jì)任務(wù)書課程石油工程課程設(shè)計(jì)題目瑞利(Rayleigh)模型的應(yīng)用專業(yè)石油工程姓名學(xué)號(hào)主要內(nèi)容、基本要求、主要參考資料等主要內(nèi)容：推導(dǎo)瑞利預(yù)測(cè)模型；根據(jù)油氣田實(shí)際生產(chǎn)數(shù)據(jù)，進(jìn)行線性回歸，求得直線的截距和斜率；根據(jù)相關(guān)公式，確定瑞利模型常數(shù)a和c；計(jì)算油田年產(chǎn)油量；計(jì)算油田累積產(chǎn)油量；計(jì)算可采儲(chǔ)量；計(jì)算最高年產(chǎn)量發(fā)生的時(shí)間；計(jì)算最高年產(chǎn)量；計(jì)算油田達(dá)到最高產(chǎn)量時(shí)的累積產(chǎn)量；繪制油田實(shí)際年產(chǎn)量與預(yù)測(cè)產(chǎn)量對(duì)比曲線；繪制油田實(shí)際累積產(chǎn)量與預(yù)測(cè)累積產(chǎn)量對(duì)比曲線?；疽螅夯A(chǔ)數(shù)據(jù):某油田或區(qū)塊的開發(fā)數(shù)據(jù)。2)設(shè)計(jì)要求該專題設(shè)計(jì)最終要求是，學(xué)生通過自選上述基礎(chǔ)數(shù)據(jù)，利用所學(xué)知識(shí)完成規(guī)定設(shè)計(jì)內(nèi)容，編制相應(yīng)軟件，并提交規(guī)范設(shè)計(jì)報(bào)告。主要參考資料陳元千油藏工程實(shí)踐[M].石油工業(yè)出版社.完成期限2021年7月1日～2021年7月19日指導(dǎo)教師王立軍專業(yè)負(fù)責(zé)人王立軍2021年7月19日目錄TOC\o"1-2"\h\z\u1前言 31.1設(shè)計(jì)的目的意義 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é)生通過自選一組數(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í)際問題相結(jié)合，通過該專題設(shè)計(jì)的訓(xùn)練，加強(qiáng)學(xué)生理論知識(shí)運(yùn)用能力，計(jì)算機(jī)技術(shù)應(yīng)用能力及解決實(shí)際問題的工程應(yīng)用能力。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)容如下：推導(dǎo)瑞利預(yù)測(cè)模型；根據(jù)油氣田實(shí)際生產(chǎn)數(shù)據(jù)，進(jìn)行線性回歸，求得直線的截距和斜率；根據(jù)相關(guān)公式，確定瑞利模型常數(shù)a和c；計(jì)算油田年產(chǎn)油量；計(jì)算油田累積產(chǎn)油量；計(jì)算可采儲(chǔ)量；計(jì)算最高年產(chǎn)量發(fā)生的時(shí)間；計(jì)算最高年產(chǎn)量；計(jì)算油田達(dá)到最高產(chǎn)量時(shí)的累積產(chǎn)量；繪制油田實(shí)際年產(chǎn)量與預(yù)測(cè)產(chǎn)量對(duì)比曲線；繪制油田實(shí)際累積產(chǎn)量與預(yù)測(cè)累積產(chǎn)量對(duì)比曲線。2基礎(chǔ)數(shù)據(jù)中國華北地區(qū)的任丘霧迷山古潛山油田于1975年7月被發(fā)現(xiàn)，1976年投產(chǎn)。其油藏類型為塊狀底水潛山，孔、洞、縫發(fā)育，儲(chǔ)層位于震旦系，以白云巖為主。油藏埋藏深度為2588~3510m，含油面積為56.9km2，油層有效厚度為272.0m.，有效孔隙度為6%，空氣滲透率為106×10-3μm2。原始地層壓力為32.6MPa，原始飽和壓力為1.3MPa，原始?xì)庥捅葹?.4m3/t，地面原油密度為0.889g/cm3，地層油粘度為8.20mPa·s。開始上報(bào)的原油地質(zhì)儲(chǔ)量為53038×104t，可采儲(chǔ)量為11760×104表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.05198814209.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ì)于油氣田來說，開發(fā)時(shí)間從0到∞之間的累計(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)量的模型為（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)的累積產(chǎn)量，104t/a（油），108m3/a（氣）。將式（3-7）代入式（3-16）得（3-17）已知，因此式（3-17）又可寫為（3-18）由式（3-18）可以看出，對(duì)于瑞利模型來說，當(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)行線性回歸之后，可以求得直線的截距和斜率的數(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)量為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.67197841118.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í)間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é)論通過本次課程設(shè)計(jì)，得出如下認(rèn)識(shí)與結(jié)論：（1）通過此次課程設(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）通過計(jì)算機(jī)的編程學(xué)習(xí)與操作，使我對(duì)VB語言有了更進(jìn)一步的認(rèn)識(shí)和了解。同時(shí)，提高了使用計(jì)算機(jī)技應(yīng)用的能力以及解決實(shí)際問題的工程應(yīng)用能力，使理論與實(shí)踐相結(jié)合。（4）這次課程設(shè)計(jì)，使我認(rèn)識(shí)到學(xué)好計(jì)算機(jī)的重要性與實(shí)踐操作的必要性，在老師的指導(dǎo)下和與同學(xué)的討論中，不斷地發(fā)現(xiàn)問題、不斷地解決問題、不斷地領(lǐng)悟、不斷地獲取，使自己的程序不斷完善。（5）通過這次的課程設(shè)計(jì)，培養(yǎng)了我獨(dú)立思考和動(dòng)手操作的能力，也使我在日常的學(xué)習(xí)與討論當(dāng)中增長了許多課外知識(shí)，認(rèn)識(shí)到數(shù)學(xué)模型對(duì)油氣田的產(chǎn)量預(yù)測(cè)的重要性。雖然在課程設(shè)計(jì)的過程中遇到了一些問題，但經(jīng)過反復(fù)的思考與檢查找出原因所在并進(jìn)行改正。參考文獻(xiàn)[1]陳元千.瑞利模型的完善推導(dǎo)與應(yīng)用[J].油氣地質(zhì)與采收率，2021，11（4）：39-41.[2]陳元千油藏工程實(shí)踐[M].石油工業(yè)出版社,2021附錄程序PrivateSubCommand2_Click()Dima()AsStringDimb()AsStringDimc()AsStringDimd()AsStringDimiAsIntegeri=1Open"d:\aa.txt"ForInputAs#1OnErrorGoTohDoWhileNotEOF(1)ReDimPreservea(i)ReDimPreserveb(i)ReDimPreservec(i)ReDimPreserved(i)Input#1,a(i),b(i),c(i),d(i)i=i+1Looph:Close#1WithMSFlexGrid1.Rows=UBound(a)Fori=1ToUBound(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)NextiEndWithEndSubPrivateSubCommand1_Click()Open"d:\aa.txt"ForOutputAs#1WithMSFlexGrid1k=.RowsFori=0Tok-1strTemp=""Forj=0To.Cols-1DoEventsstrTemp=strTemp&.TextMatrix(i,j)&","NextjPrint#1,Left(strTemp,Len(strTemp)-1)NextiEndWithClose#1EndSubPrivateSubCommand3_Click()Dima()AsStringDimb()AsStringDimc()AsStringDimd()AsStringDime()AsStringDimf()AsStringDimiAsIntegeri=1Open"d:\aa.txt"ForInputAs#1OnErrorGoTohDoWhileNotEOF(1)ReDimPreservea(i)ReDimPreserveb(i)ReDimPreservec(i)ReDimPreserved(i)Input#1,a(i),b(i),c(i),d(i)i=i+1Looph:Close#1n=UBound(a)ReDime(1Ton)AsStringReDimf(1Ton)AsStringFori=2ToUBound(a)e(i)=Str(Log(Val(c(i))/Val(b(i)))/Log(Val(10)))f(i)=Str(Val(b(i))^2)NextiWithMSFlexGrid3MSFlexGrid3.TextMatrix(0,0)="年份"MSFlexGrid3.TextMatrix(0,1)="LogQ/t"MSFlexGrid3.TextMatrix(0,2)="t2".Rows=UBound(a)Fori=2ToUBound(a)MSFlexGrid3.TextMatrix(i-1,0)=a(i)MSFlexGrid3.TextMatrix(i-1,1)=e(i)MSFlexGrid3.TextMatrix(i-1,2)=f(i)NextiEndWithEndSubPrivateSubCommand4_Click()k=10^Val(Text2.Text)l=1/(2.303*Val(Text3.Text))Dima()AsStringDimb()AsStringDimc()AsStringDimd()AsStringDimg()AsSingleDimiAsIntegeri=1Open"d:\aa.txt"ForInputAs#1OnErrorGoTohDoWhileNotEOF(1)ReDimPreservea(i)ReDimPreserveb(i)ReDimPreservec(i)ReDimPreserved(i)Input#1,a(i),b(i),c(i),d(i)i=i+1Looph:Close#1n=UBound(a)ReDimg(2Ton+2)AsSingleFori=2Tong(i)=k*Val(b(i))*Exp(-(Val(b(i))^2/l))NextiWithMSFlexGrid4.Rows=UBound(a)Fori=0Ton-2MSFlexGrid4.TextMatrix(i+1,4)=g(i+2)NextiEndWithReDimMyData(n-2,3)AsDouble'x軸坐標(biāo)值Y軸坐標(biāo)值Fori=0To15MyData(i,2)=b(i+2):MyData(i,3)=g(i+2)'本句代表了：第一點(diǎn)數(shù)據(jù)的X軸坐標(biāo)，Y軸坐標(biāo)NextiFori=0To15MyData(i,0)=b(i+2):MyData(i,1)=c(i+2)Nexti'波形圖外觀設(shè)置WithMSChart3.TitleText="實(shí)際與預(yù)測(cè)年產(chǎn)量隨時(shí)間變化曲線"''設(shè)置圖線的外觀.Plot.SeriesCollection(1).Pen.Width=30.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=2021'//設(shè)置最小值.Plot.Axis(VtChAxisIdY).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(VtChAxisIdY).ValueScale.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.ChartData=MyData'數(shù)據(jù)DimIndexAsIntegerIndex=1'設(shè)置圖標(biāo)類型為非自動(dòng)選擇，這樣可以設(shè)置為自己需要的類型MSChart3.Plot.SeriesCollection(Index).ShowLine=False'去除連線MSChart3.Plot.SeriesCollection.Item(Index).SeriesMarker.Auto=FalseWithMSChart3.Plot.SeriesCollection.Item(Index).DataPoints.Item(-1).Marker.Visible=True'散點(diǎn)的點(diǎn)是否可見.Size=80'散點(diǎn)的大小.Style=VtMarkerStyleFilledDiamond'選擇菱形的.FillColor.Automatic=False'顏色不自動(dòng)填充.FillColor.Set400,40,400'設(shè)置為黑色.Pen.VtColor.Set255,40,400'同樣這個(gè)也必須設(shè)置，否則中間有默認(rèn)的顏色EndWithEndWithEndSubPrivateSubCommand5_Click()k=10^Val(Text2.Text)l=1/(2.303*Val(Text3.Text))Dima()AsStringDimb()AsStringDimc()AsStringDimd()AsStringDimh()AsSingleDimiAsIntegeri=1Open"d:\aa.txt"ForInputAs#1OnErrorGoTohDoWhileNotEOF(1)ReDimPreservea(i)ReDimPreserveb(i)ReDimPreservec(i)ReDimPreserved(i)Input#1,a(i),b(i),c(i),d(i)i=i+1Looph:Close#1n=UBound(a)ReDimh(2Ton+2)AsSingleFori=2Tonh(i)=k*l/2*(1-Exp(-(Val(b(i))^2/l)))NextiWithMSFlexGrid4.Rows=UBound(a)Fori=0Ton-2MSFlexGrid4.TextMatrix(i+1,5)=h(i+2)NextiEndWithReDimMyData(n-2,3)AsDouble'x軸坐標(biāo)值Y軸坐標(biāo)值Fori=0To15MyData(i,2)=b(i+2):MyData(i,3)=h(i+2)'本句代表了：第一點(diǎn)數(shù)據(jù)的X軸坐標(biāo)，Y軸坐標(biāo)NextiFori=0To15MyData(i,0)=b(i+2):MyData(i,1)=d(i+2)Nexti'波形圖外觀設(shè)置WithMSChart2.TitleText="實(shí)際與預(yù)測(cè)累計(jì)產(chǎn)量隨時(shí)間變化曲線"''設(shè)置圖線的外觀.Plot.SeriesCollection(1).Pen.Width=30.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=12021'//設(shè)置最小值.Plot.Axis(VtChAxisIdY).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(VtChAxisIdY).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)圖MSChart2.ChartData=MyData'數(shù)據(jù)DimIndexAsIntegerIndex=1'設(shè)置圖標(biāo)類型為非自動(dòng)選擇，這樣可以設(shè)置為自己需要的類型MSChart2.Plot.SeriesCollection(Index).ShowLine=False'去除連線MSChart2.Plot.SeriesCollection.Item(Index).SeriesMarker.Auto=FalseWithMSChart2.Plot.SeriesCollection.Item(Index).DataPoints.Item(-1).Marker.Visible=True'散點(diǎn)的點(diǎn)是否可見.Size=80'散點(diǎn)的大小.Style=VtMarkerStyleFilledDiamond'選擇菱形的.FillColor.Automatic=False'顏色不自動(dòng)填充.FillColor.Set400,40,400'設(shè)置為黑色.Pen.VtColor.Set55,40,400'同樣這個(gè)也必須設(shè)置，否則中間有默認(rèn)的顏色EndWithEndWithEndSubPrivateSubCommand6_Click()Dima()AsStringDimb()AsStringDimc()AsStringDimd()AsStringDime()AsStringDimf()AsStringDimiAsIntegeri=1Open"d:\aa.txt"ForInputAs#1OnErrorGoTohDoWhileNotEOF(1)ReDimPreservea(i)ReDimPreserveb(i)ReDimPreservec(i)ReDimPreserved(i)Input#1,a(i),b(i),c(i),d(i)i=i+1Looph:Close#1n=UBound(a)ReDime(2Ton)AsStringReDimf(2Ton)AsStringFori=2ToUBound(a)e(i)=Str(Log(Val(c(i))/Val(b(i)))/Log(Val(10)))f(i)=Str(Val(b(i))^2)NextiReDimMyData(n-2,1)AsDouble'x軸坐標(biāo)值Y軸坐標(biāo)值Fori=0ToUBound(a)-2MyData(i,0)=f(i+2):MyData(i,1)=e(i+2)'本句代表了：第一點(diǎn)數(shù)據(jù)的X軸坐標(biāo),Y軸坐標(biāo)Nexti'波形圖外觀設(shè)置WithMSChart1.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(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).ValueScale.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=VtPenStyleDotted.Plot.Axis(VtChAxisIdY).AxisGrid.MajorPen.Style=VtPenStyleDottedMSChart1.Plot.AutoLayout=FalseMSChart1.Plot.UniformAxis=FalseMSChart1.chartType=VtChChartType2dXY'設(shè)置圖形為二維散點(diǎn)圖MSChart1.ChartData=MyData'數(shù)據(jù)DimIndexAsIntegerIndex=1'設(shè)置圖標(biāo)類型為非自動(dòng)選擇，這樣可以設(shè)置為自己需要的類型MSChart1.Plot.SeriesCollection(Index).ShowLine=False'去除連線MSChart1.Plot.SeriesCollection.Item(Index).SeriesMarker.Auto=FalseWithMSChart1.Plot.SeriesCollection.Item(Index).DataPoints.Item(-1).Marker.Visible=True'散點(diǎn)的點(diǎn)是否可見.Size=80'散點(diǎn)的大小.Style=VtMarkerStyleFilledDiamond'選擇菱形的.FillColor.Automatic=False'顏色不自動(dòng)填充.FillColor.Set400,40,400'設(shè)置為黑色.Pen.VtColor.Set255,40,400'同樣這個(gè)也必須設(shè)置，否則中間有默認(rèn)的顏色EndWithEndWithEndSubPrivateSubCommand7_Click()Dima()AsStringDimb()AsStringDimc()AsStringDimd()AsStringDime()AsStringDimf()AsStringDimiAsIntegeri=1Open"d:\aa.txt"ForInputAs#1OnErrorGoTohDoWhileNotEOF(1)ReDimPreservea(i)ReDimPreserveb(i)ReDimPreservec(i)ReDimPreserved(i)Input#1,a(i),b(i),c(i),d(i)i=i+1Looph:Close#1n=UBound(a)ReDime(2Ton)AsStringReDimf(2Ton)AsStringFori=2ToUBound(a)e(i)=Str(Log(Val(c(i))/Val(b(i)))/Log(Val(10)))f(i)=Str(Val(b(i))^2)NextiReDimMyData(n-2,3)AsDouble'x軸坐標(biāo)值Y軸坐標(biāo)值Fori=0ToUBound(a)-2MyData(i,0)=f(i+2):MyData(i,1)=e(i+2)'本句代表了：第一點(diǎn)數(shù)據(jù)的X軸坐標(biāo)，Y軸坐標(biāo)NextiMyData(0,2)=Text4.Text:MyData(0,3)=Text5.TextFori=1ToUBound(a)-2MyData(i,2)=Text6.Text:MyData(i,3)=Text7.TextNexti'波形圖外觀設(shè)置WithMSChart1.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(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).ValueScale.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=VtPenStyleDotted.Plot.Axis(VtChAxisIdY).AxisGrid.MajorPen.Style=VtPenStyleDottedMSChart1.Plot.AutoLayout=FalseMSChart1.Plot.UniformAxis=FalseMSChart1.chartType=VtChChartType2dXY'設(shè)置圖形為二維散點(diǎn)圖MSChart1.ChartData=MyData'數(shù)據(jù).Plot.Axis(VtChAxisIdX,0).AxisTitle="t2".Plot.Axis(VtChAxisIdY,0).AxisTitle="LogQ/t"DimIndexAsIntegerIndex=1'設(shè)置圖標(biāo)類型為非自動(dòng)選擇，這樣可以設(shè)置為自己需要的類型MSChart1.Plot.SeriesCollection(Index).ShowLine=False'去除連線MSChart1.Plot.SeriesCollection.Item(Index).SeriesMarker.Auto=FalseWithMSChart1.Plot.SeriesCollection.Item(Index).DataPoints.Item(-1).Marker.Visible=True'散點(diǎn)的點(diǎn)是否可見.Size=80'散點(diǎn)的大小.Style=VtMarkerStyleFilledDiamond'選擇菱形的.FillColor.Automatic=False'顏色不自動(dòng)填充.FillColor.Set400,40,400'設(shè)置為黑色.Pen.VtColor.Set400,40,400'同樣這個(gè)也必須設(shè)置，否則中間有默認(rèn)的顏色EndWithEndWithText3.Text=(Text5.Text-Text7.Text)/(Text6.Text-Text4.Text)Text2.Text=Text4.Text*Text3.Text+Text5.TextEndSubPrivateSubCommand8_Click()Text4.Text=7.14Text5.Text=2.53Text6.Text=241.43Text7.Text=0.76EndSubPrivateSubCommand9_Click()k=10^Val(Text2.Text)l=1/(2.303*Val(Text3.Text))Text8.Text=k*l/2Text9.Text=0.4289*k*l^(1/2)Text10.Text=0.707*l^(1/2)Text11.Text=0.4*k*l/2Dima()AsStringDimb()AsStringDimc()AsStringDimd()AsStringDimiAsIntegeri=1Open"d:\aa.txt"ForInputAs#1OnErrorGoTohDoWhileNotEOF(1)ReDimPreservea(i)ReDimPreserveb(i)ReDimPreservec(i)ReDimPreserved(i)Input#1,a(i),b(i),c(i),d(i)i=i+1Looph:Close#1WithMSFlexGrid4.Rows=UBound(a)Fori=1ToUBound(a)MSFlexGrid4.TextMatrix(i-1,0)=a(i)MSFlexGrid4.TextMatrix(i-1,1)=b(i)MSFlexGrid4.TextMatrix(i-1,2)=c(i)MSFlexGrid4.TextMatrix(i-1,3)=d(i)NextiEndWithEndSubPrivateSubForm_Load()DimMsgdAsStringWithMSFlexGrid1.ClearMsgd="年份|開發(fā)時(shí)間|年產(chǎn)量|累積產(chǎn)量".FormatString=MsgdEndWithWithMSFlexGrid4.ClearMsgd="年份|開發(fā)時(shí)間|年產(chǎn)量|累積產(chǎn)量|預(yù)測(cè)年產(chǎn)量|預(yù)測(cè)累積產(chǎn)量".FormatString=MsgdEndWithText1.Visible=FalseWithMSFlexGrid1.ColWidth(0)=1000'表格行寬.ColWidth(1)=1200.ColWidth(2)=1200.ColWidth(3)=1200EndWithWithMSFlexGrid4.ColWidth(0)=1000'表格行寬.ColWidth(1)=1200.ColWidth(2)=1200.ColWidth(3)=1200.ColWidth(4)=1200.ColWidth(5)=1200EndWithSkinH_AttachEndSubSubShowTextBox()WithMSFlexGrid1'隱藏文本框，設(shè)置高度和寬度Text1.Visible=FalseText1.Height=.RowHeight(.Row)-(Screen.TwipsPerPixelY)*2Text1.Width=.ColWidth(.Col)'計(jì)算文本框左坐標(biāo)Text1.Left=.CellLeft+.LeftText1.Top=.CellTop+.TopText1.Visible=TrueText1.SetFocusEndWithEndSubPrivateSubMSChart1_MouseDown(ButtonAsInteger,ShiftAsInteger,XAsSingle,YAsSingle)IfButton<>2ThenExitSubText4.Text=(X-1185)/21-16.42857Text5.Text=(3600-Y)/1080EndSubPrivateSubMSChart1_Mousemove(ButtonAsInteger,ShiftAsInteger,XAsSingle,YAsSingle)IfButton<>2ThenExitSubText6.Text=(X-1185)/21-16.42857Text7.Text=(3600-Y)/1080EndSubPrivateSubMSFlexGrid1_KeyPress(KeyAsciiAsInteger)DimcharAsStringIfKeyAscii=13ThenText1.Text=MSFlexGrid1.TextText1.SelStart=Len(Text1.Text)Elsechar=Chr$(KeyAscii)Text1.Text=charText1.SelStart=1EndIfShowTextBoxKeyAscii=0EndSubPrivateSubMSFlexGrid1_LeaveCell()MSFlexGrid1.Text=Text1.TextIfMSFlexGrid1.Col=0AndMSFlexGrid1.Row<>0And_MSFlexGrid1.Row=MSFlexGrid1.Rows-1AndMSFlexGrid1.Text<>""ThenMSFlexGrid1.Rows=MSFlexGrid1.Rows+1EndIfEndSubPrivateSubMSFlexGrid1_MouseDown(ButtonAsInteger,ShiftAsInteger,_XAsSingle,YAsSingle)Text1.Visible=FalseEndSubPrivateSubMSFlexGrid1_RowColChange()Text1.Text=MSFlexGrid1.TextEndSubPrivateSubMSFlexGrid1_DblClick()IfMSFlexGrid1.Row>0AndMSFlexGrid1.Col=0ThenMSFlexGrid1_KeyPress13EndSubPrivateSubText1_KeyPress(KeyAsciiAsInteger)IfKeyAscii=13ThenMSFlexGrid1.Text=Text1.TextText1.Visible=FalseMSFlexGrid1.SetFocusIfMSFlexGrid1.Col<(MSFlexGrid1.Cols-1)ThenMSFlexGrid1.Col=MSFlexGrid1.Col+1ElseIfMSFlexGrid1.Row<MSFlexGrid1.Rows-1ThenMSFlexGrid1.Row=MSFlexGrid1.Row+1MSFlexGrid1.Col=0EndIfKeyAscii=0EndIfEndSub模塊程序：PublicDeclareFunctionSkinH_AttachLib"SkinH_VB6.dll"()AsLong運(yùn)行結(jié)果數(shù)據(jù)處理及線性回歸年產(chǎn)量與累積產(chǎn)量的曲線PAGE50東北石油大學(xué)課程設(shè)計(jì)成績?cè)u(píng)價(jià)表課程名稱石油工程課程設(shè)計(jì)題目名稱瑞利(Rayleigh)模型的應(yīng)用學(xué)生姓名學(xué)號(hào)指導(dǎo)教師姓名王立軍職稱教授序號(hào)評(píng)價(jià)項(xiàng)目指標(biāo)滿分評(píng)分1工作量、工作態(tài)度和出勤率按期圓滿的完成了規(guī)定的任務(wù)，難易程度和工作量符合教學(xué)要求，工作努力，遵守紀(jì)律，出勤率高，工作作風(fēng)嚴(yán)謹(jǐn)，善于與他人合作。202課程設(shè)計(jì)質(zhì)量課程設(shè)計(jì)選題合理，計(jì)算過程簡(jiǎn)練準(zhǔn)確，分析問題思路清晰，結(jié)構(gòu)嚴(yán)謹(jǐn)，文理通順，撰寫規(guī)范，圖表完備正確。453創(chuàng)新工作中有創(chuàng)新意識(shí)，對(duì)前人工作有一些改進(jìn)或有一定應(yīng)用價(jià)值。54答辯能正確回答指導(dǎo)教師所提出的問題。30總分評(píng)語：指導(dǎo)教師：2021年7月19日

論大學(xué)生寫作能力寫作能力是對(duì)自己所積累的信息進(jìn)行選擇、提取、加工、改造并將之形成為書面文字的能力。積累是