瑞利(rayleigh)模型的应用

石油工程学院课程设计（报告）PAGE东北石油大学课程设计课程课程石油工程课程设计题目瑞利(Rayleigh)模型的应用院系石油工程学院专业班级油工班学生姓名学生学号指导教师王立军20年7月19日PAGEPAGE17东北石油大学课程设计任务书课程石油工程课程设计题目瑞利(Rayleigh)模型的应用专业石油工程姓名学号主要内容、基本要求、主要参考资料等主要内容：推导瑞利预测模型；根据油气田实际生产数据，进行线性回归，求得直线的截距和斜率；根据相关公式，确定瑞利模型常数a和c；计算油田年产油量；计算油田累积产油量；计算可采储量；计算最高年产量发生的时间；计算最高年产量；计算油田达到最高产量时的累积产量；绘制油田实际年产量与预测产量对比曲线；绘制油田实际累积产量与预测累积产量对比曲线。基本要求：基础数据:某油田或区块的开发数据。2)设计要求该专题设计最终要求是，学生通过自选上述基础数据，利用所学知识完成规定设计内容，编制相应软件，并提交规范设计报告。主要参考资料陈元千油藏工程实践[M].石油工业出版社.完成期限2021年7月1日～2021年7月19日指导教师王立军专业负责人王立军2021年7月19日目录TOC\o"1-2"\h\z\u1前言 31.1设计的目的意义 31.2设计的主要内容 32基础数据 43基础理论 53.1瑞利(Rayleigh)模型的应用模型的建立 53.2模型的求解方法 74设计结果 8认识与结论 12参考文献 13附录 141前言1.1设计的目的意义本课题的目的是让学生通过自选一组数据，利用所学专业知识在指导教师的指导下独立完成对某一油田或区块开发指标的预测。本课题要求学生对瑞利(Rayleigh)模型的关系式进行推导，其结果包括，对油气田的年产量、累积产量、可采储量、最高年产量、最高年产量发生的时间的预测方法。从而将理论知识和实际问题相结合，通过该专题设计的训练，加强学生理论知识运用能力，计算机技术应用能力及解决实际问题的工程应用能力。1.2设计的主要内容根据已有的基础数据，利用所学的专业知识，在指导教师指导下独立完成并提交一个油田或一个区块油田开发指标预测结果，设计主要内容如下：推导瑞利预测模型；根据油气田实际生产数据，进行线性回归，求得直线的截距和斜率；根据相关公式，确定瑞利模型常数a和c；计算油田年产油量；计算油田累积产油量；计算可采储量；计算最高年产量发生的时间；计算最高年产量；计算油田达到最高产量时的累积产量；绘制油田实际年产量与预测产量对比曲线；绘制油田实际累积产量与预测累积产量对比曲线。2基础数据中国华北地区的任丘雾迷山古潜山油田于1975年7月被发现，1976年投产。其油藏类型为块状底水潜山，孔、洞、缝发育，储层位于震旦系，以白云岩为主。油藏埋藏深度为2588~3510m，含油面积为56.9km2，油层有效厚度为272.0m.，有效孔隙度为6%，空气渗透率为106×10-3μm2。原始地层压力为32.6MPa，原始饱和压力为1.3MPa，原始气油比为4.4m3/t，地面原油密度为0.889g/cm3，地层油粘度为8.20mPa·s。开始上报的原油地质储量为53038×104t，可采储量为11760×104表2-1开发数据表年份开发时间(a)年产量(104t/a)累积产量(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基础理论瑞利(Rayleigh)模型能够很好地预测对油气田的年产量、累积产量、可采储量、最高年产量、最高年产量发生的时间，以及最高年产量对应的累计产量。3.1瑞利(Rayleigh)模型的应用模型的建立在数理统计中瑞利分布的分布密度表示为：(3-1)式中f(x)——瑞利分布的分布密度函数；x——分布变量，其值大于等于0；μ——控制分布峰值的参数。由式(3-l)表示的瑞利分布，在x从0到∞区间内的累积分布（或分布函数）等于1，因此可以得到关系式如下（3-2）对于油气田来说，开发时间从0到∞之间的累计产量可视为油气田的可采储量，可表示为(3-3)式中：NR为油气田的可采储量，104t（油），108m3（气）；NP（t）为累积产量，104t（油），108m3（气）；Q年产量，104t/a（油），108m3/a（气）；t为开发时间，a。将式（3-3）等号两端同除以NR得（3-4）由式（3-2）与式（3-4）对比可以看出，，因此，若将式（3-1）转换为预测油气田产量分布的模型时，需将该式的右端乘以可采储量。该可采储量可视为，将数理统计中的理论分布转为预测模型的转换常数。由此，式（3-1）可写为（3-5）设，则；又设和得（3-6）再作如下设定（3-7）则得到预测油气田产量的模型为（3-8）油气田的累积产量表示为（3-9）将式（3-8）代入式（3-9）进行积分得（3-10）再讲式（3-7）代入式（3-10），得到预测油气田累积产量的模型为（3-11）对式（3-8）两边取导数，整理得（3-12）当时，必然有（3-13）由式（3-13）可以得到最高年产量发生的时间为（3-14）式中：为最高年产量发生的时间，a将式（3-14）代入式（3-8）的得到预测油气田最高年产量的公式为（3-15）式中：为最高年产量，104t/a（油），108m3/a（气）。再将式（3-14）代入式（3-11）得到预测最高年产量对应的累积产量为（3-16）式中：为最高年产量对应的累积产量，104t/a（油），108m3/a（气）。将式（3-7）代入式（3-16）得（3-17）已知，因此式（3-17）又可写为（3-18）由式（3-18）可以看出，对于瑞利模型来说，当可采储量的40%时，油田进入递减期。3.2模型的求解方法由上述推导的结果可以看出，要想利用瑞利模型进行各项预测，就必须首先确定模型常数a和c的数值。为此，将式（3-18）等号两端先除以t再取常用对数，得（3-19）若设：（3-20）（3-21）则得：（3-22）由式（3-22）看出，Q/t与之间具有半对数直线关系。当根据油气田的实际生产数据，由式（3-22）进行线性回归之后，可以求得直线的截距和斜率的数值，并由式（3-22）和式（21）确定模型常数a和c的数值。4设计结果处理表1中的数据得到和相应的，按照（3-22）式的直线关系绘于图4-1，得到了一条很好的直线。经线性回归求得直线的截距α=2.5816；直线的斜率β=0.007528；直线的相关系数r=0.9945。将α和β的数值分别代入式（3-20）和式（3-21）求得a为387.59，c为57.76。图4-1该油田与的半对数关系将a和c的数值代入式（3-7）得到该油田的可采储量为11020×104t；将c的值代入式（3-14）得到最高年产量发生的时间为5.3692a；将a和c的数值代入式（3-15）得到该油田的最高年产量为1242.8×104t；将该油田可采储量的数值代入式（3-18）得到该当油田年产量达到最高年产量时的累积产量为4401.14×104t。将a和c的数值代入式（3-8）得，预测该油田年产量的公式为（3-23）将a和c的数值代入式（3-10）得到预测该油田累积产量的公式为（3-24）由式（3-23）和式（3-24）预测该油田的产量和累积产量如表4-1表4-1开发数据表年份开发时间(a)预测年产量(104t/a)预测累积产量(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根据表4-1中数据，在同一坐标系中绘制实际产量与预测产量及其对应时间t之间的关系曲线，绘于图4-2。图4-2实际与预测产量随时间变化曲线在同一坐标系中绘制实际累计产量与预测累计产量及其对应时间t之间的关系曲线，绘于图4-3。图4-3实际与预测累计产量随时间变化曲线认识与结论通过本次课程设计，得出如下认识与结论：（1）通过此次课程设计，基本掌握了瑞利(Rayleigh)模型的应用方法，并能初步对油田进行预测分析。（2）能熟练掌握、运用瑞利(Rayleigh)模型原理、推导及公式，可以对油气田的年产量、累计产量、可采储量、最高年产量、最高年产量发生的时间以及与最高年产量对应的累积产量进行预测。（3）通过计算机的编程学习与操作，使我对VB语言有了更进一步的认识和了解。同时，提高了使用计算机技应用的能力以及解决实际问题的工程应用能力，使理论与实践相结合。（4）这次课程设计，使我认识到学好计算机的重要性与实践操作的必要性，在老师的指导下和与同学的讨论中，不断地发现问题、不断地解决问题、不断地领悟、不断地获取，使自己的程序不断完善。（5）通过这次的课程设计，培养了我独立思考和动手操作的能力，也使我在日常的学习与讨论当中增长了许多课外知识，认识到数学模型对油气田的产量预测的重要性。虽然在课程设计的过程中遇到了一些问题，但经过反复的思考与检查找出原因所在并进行改正。参考文献[1]陈元千.瑞利模型的完善推导与应用[J].油气地质与采收率，2021，11（4）：39-41.[2]陈元千油藏工程实践[M].石油工业出版社,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轴坐标值Y轴坐标值Fori=0To15MyData(i,2)=b(i+2):MyData(i,3)=g(i+2)'本句代表了：第一点数据的X轴坐标，Y轴坐标NextiFori=0To15MyData(i,0)=b(i+2):MyData(i,1)=c(i+2)Nexti'波形图外观设置WithMSChart3.TitleText="实际与预测年产量随时间变化曲线"''设置图线的外观.Plot.SeriesCollection(1).Pen.Width=30.Plot.SeriesCollection(1).Pen.Style=VtPenStyleSolid''设置XY轴.Plot.Axis(VtChAxisIdX).ValueScale.Auto=False.Plot.Axis(VtChAxisIdY).ValueScale.Auto=False'//设置最大值.Plot.Axis(VtChAxisIdX).ValueScale.Maximum=18.Plot.Axis(VtChAxisIdY).ValueScale.Maximum=2021'//设置最小值.Plot.Axis(VtChAxisIdY).ValueScale.Minimum=0.Plot.Axis(VtChAxisIdX).ValueScale.Minimum=0'//.Plot.Axis(VtChAxisIdX).ValueScale.MajorDivision=12'X轴主要网格数量.Plot.Axis(VtChAxisIdY).ValueScale.MajorDivision=10'Y轴主要网格数量.Plot.Axis(VtChAxisIdX).ValueScale.MinorDivision=0'X轴次要网格数量.Plot.Axis(VtChAxisIdY).ValueScale.MinorDivision=0'Y轴次要网格数量.Plot.Axis(VtChAxisIdX).AxisGrid.MajorPen.Style=VtPenStyleDotted.Plot.Axis(VtChAxisIdY).AxisGrid.MajorPen.Style=VtPenStyleDottedMSChart3.Plot.AutoLayout=FalseMSChart3.Plot.UniformAxis=FalseMSChart3.chartType=VtChChartType2dXY'设置图形为二维散点图MSChart3.ChartData=MyData'数据DimIndexAsIntegerIndex=1'设置图标类型为非自动选择，这样可以设置为自己需要的类型MSChart3.Plot.SeriesCollection(Index).ShowLine=False'去除连线MSChart3.Plot.SeriesCollection.Item(Index).SeriesMarker.Auto=FalseWithMSChart3.Plot.SeriesCollection.Item(Index).DataPoints.Item(-1).Marker.Visible=True'散点的点是否可见.Size=80'散点的大小.Style=VtMarkerStyleFilledDiamond'选择菱形的.FillColor.Automatic=False'颜色不自动填充.FillColor.Set400,40,400'设置为黑色.Pen.VtColor.Set255,40,400'同样这个也必须设置，否则中间有默认的颜色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轴坐标值Y轴坐标值Fori=0To15MyData(i,2)=b(i+2):MyData(i,3)=h(i+2)'本句代表了：第一点数据的X轴坐标，Y轴坐标NextiFori=0To15MyData(i,0)=b(i+2):MyData(i,1)=d(i+2)Nexti'波形图外观设置WithMSChart2.TitleText="实际与预测累计产量随时间变化曲线"''设置图线的外观.Plot.SeriesCollection(1).Pen.Width=30.Plot.SeriesCollection(1).Pen.Style=VtPenStyleSolid''设置XY轴.Plot.Axis(VtChAxisIdX).ValueScale.Auto=False.Plot.Axis(VtChAxisIdY).ValueScale.Auto=False'//设置最大值.Plot.Axis(VtChAxisIdX).ValueScale.Maximum=18.Plot.Axis(VtChAxisIdY).ValueScale.Maximum=12021'//设置最小值.Plot.Axis(VtChAxisIdY).ValueScale.Minimum=0.Plot.Axis(VtChAxisIdX).ValueScale.Minimum=0'//.Plot.Axis(VtChAxisIdX).ValueScale.MajorDivision=12'X轴主要网格数量.Plot.Axis(VtChAxisIdY).ValueScale.MajorDivision=10'Y轴主要网格数量.Plot.Axis(VtChAxisIdX).ValueScale.MinorDivision=0'X轴次要网格数量.Plot.Axis(VtChAxisIdY).ValueScale.MinorDivision=0'Y轴次要网格数量.Plot.Axis(VtChAxisIdX).AxisGrid.MajorPen.Style=VtPenStyleDotted.Plot.Axis(VtChAxisIdY).AxisGrid.MajorPen.Style=VtPenStyleDottedMSChart2.Plot.AutoLayout=FalseMSChart2.Plot.UniformAxis=FalseMSChart2.chartType=VtChChartType2dXY'设置图形为二维散点图MSChart2.ChartData=MyData'数据DimIndexAsIntegerIndex=1'设置图标类型为非自动选择，这样可以设置为自己需要的类型MSChart2.Plot.SeriesCollection(Index).ShowLine=False'去除连线MSChart2.Plot.SeriesCollection.Item(Index).SeriesMarker.Auto=FalseWithMSChart2.Plot.SeriesCollection.Item(Index).DataPoints.Item(-1).Marker.Visible=True'散点的点是否可见.Size=80'散点的大小.Style=VtMarkerStyleFilledDiamond'选择菱形的.FillColor.Automatic=False'颜色不自动填充.FillColor.Set400,40,400'设置为黑色.Pen.VtColor.Set55,40,400'同样这个也必须设置，否则中间有默认的颜色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轴坐标值Y轴坐标值Fori=0ToUBound(a)-2MyData(i,0)=f(i+2):MyData(i,1)=e(i+2)'本句代表了：第一点数据的X轴坐标,Y轴坐标Nexti'波形图外观设置WithMSChart1.TitleText="Rayleigh曲线"''设置图线的外观.Plot.SeriesCollection(1).Pen.Width=30.Plot.SeriesCollection(1).Pen.Style=VtPenStyleSolid''设置XY轴.Plot.Axis(VtChAxisIdX).ValueScale.Auto=False.Plot.Axis(VtChAxisIdY).ValueScale.Auto=False'//设置最大值.Plot.Axis(VtChAxisIdX).ValueScale.Maximum=300.Plot.Axis(VtChAxisIdY).ValueScale.Maximum=3'//设置最小值.Plot.Axis(VtChAxisIdY).ValueScale.Minimum=0.Plot.Axis(VtChAxisIdX).ValueScale.Minimum=0'//.Plot.Axis(VtChAxisIdX).ValueScale.MajorDivision=20'X轴主要网格数量.Plot.Axis(VtChAxisIdY).ValueScale.MajorDivision=20'Y轴主要网格数量.Plot.Axis(VtChAxisIdX).ValueScale.MinorDivision=0'X轴次要网格数量.Plot.Axis(VtChAxisIdY).ValueScale.MinorDivision=0'Y轴次要网格数量.Plot.Axis(VtChAxisIdX).AxisGrid.MajorPen.Style=VtPenStyleDotted.Plot.Axis(VtChAxisIdY).AxisGrid.MajorPen.Style=VtPenStyleDottedMSChart1.Plot.AutoLayout=FalseMSChart1.Plot.UniformAxis=FalseMSChart1.chartType=VtChChartType2dXY'设置图形为二维散点图MSChart1.ChartData=MyData'数据DimIndexAsIntegerIndex=1'设置图标类型为非自动选择，这样可以设置为自己需要的类型MSChart1.Plot.SeriesCollection(Index).ShowLine=False'去除连线MSChart1.Plot.SeriesCollection.Item(Index).SeriesMarker.Auto=FalseWithMSChart1.Plot.SeriesCollection.Item(Index).DataPoints.Item(-1).Marker.Visible=True'散点的点是否可见.Size=80'散点的大小.Style=VtMarkerStyleFilledDiamond'选择菱形的.FillColor.Automatic=False'颜色不自动填充.FillColor.Set400,40,400'设置为黑色.Pen.VtColor.Set255,40,400'同样这个也必须设置，否则中间有默认的颜色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轴坐标值Y轴坐标值Fori=0ToUBound(a)-2MyData(i,0)=f(i+2):MyData(i,1)=e(i+2)'本句代表了：第一点数据的X轴坐标，Y轴坐标NextiMyData(0,2)=Text4.Text:MyData(0,3)=Text5.TextFori=1ToUBound(a)-2MyData(i,2)=Text6.Text:MyData(i,3)=Text7.TextNexti'波形图外观设置WithMSChart1.TitleText="Rayleigh曲线"''设置图线的外观.Plot.SeriesCollection(1).Pen.Width=30.Plot.SeriesCollection(1).Pen.Style=VtPenStyleSolid''设置XY轴.Plot.Axis(VtChAxisIdX).ValueScale.Auto=False.Plot.Axis(VtChAxisIdY).ValueScale.Auto=False'//设置最大值.Plot.Axis(VtChAxisIdX).ValueScale.Maximum=300.Plot.Axis(VtChAxisIdY).ValueScale.Maximum=3'//设置最小值.Plot.Axis(VtChAxisIdY).ValueScale.Minimum=0.Plot.Axis(VtChAxisIdX).ValueScale.Minimum=0'//.Plot.Axis(VtChAxisIdX).ValueScale.MajorDivision=20'X轴主要网格数量.Plot.Axis(VtChAxisIdY).ValueScale.MajorDivision=20'Y轴主要网格数量.Plot.Axis(VtChAxisIdX).ValueScale.MinorDivision=0'X轴次要网格数量.Plot.Axis(VtChAxisIdY).ValueScale.MinorDivision=0'Y轴次要网格数量.Plot.Axis(VtChAxisIdX).AxisGrid.MajorPen.Style=VtPenStyleDotted.Plot.Axis(VtChAxisIdY).AxisGrid.MajorPen.Style=VtPenStyleDottedMSChart1.Plot.AutoLayout=FalseMSChart1.Plot.UniformAxis=FalseMSChart1.chartType=VtChChartType2dXY'设置图形为二维散点图MSChart1.ChartData=MyData'数据.Plot.Axis(VtChAxisIdX,0).AxisTitle="t2".Plot.Axis(VtChAxisIdY,0).AxisTitle="LogQ/t"DimIndexAsIntegerIndex=1'设置图标类型为非自动选择，这样可以设置为自己需要的类型MSChart1.Plot.SeriesCollection(Index).ShowLine=False'去除连线MSChart1.Plot.SeriesCollection.Item(Index).SeriesMarker.Auto=FalseWithMSChart1.Plot.SeriesCollection.Item(Index).DataPoints.Item(-1).Marker.Visible=True'散点的点是否可见.Size=80'散点的大小.Style=VtMarkerStyleFilledDiamond'选择菱形的.FillColor.Automatic=False'颜色不自动填充.FillColor.Set400,40,400'设置为黑色.Pen.VtColor.Set400,40,400'同样这个也必须设置，否则中间有默认的颜色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="年份|开发时间|年产量|累积产量".FormatString=MsgdEndWithWithMSFlexGrid4.ClearMsgd="年份|开发时间|年产量|累积产量|预测年产量|预测累积产量".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'隐藏文本框，设置高度和宽度Text1.Visible=FalseText1.Height=.RowHeight(.Row)-(Screen.TwipsPerPixelY)*2Text1.Width=.ColWidth(.Col)'计算文本框左坐标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运行结果数据处理及线性回归年产量与累积产量的曲线PAGE50东北石油大学课程设计成绩评价表课程名称石油工程课程设计题目名称瑞利(Rayleigh)模型的应用学生姓名学号指导教师姓名王立军职称教授序号评价项目指标满分评分1工作量、工作态度和出勤率按期圆满的完成了规定的任务，难易程度和工作量符合教学要求，工作努力，遵守纪律，出勤率高，工作作风严谨，善于与他人合作。202课程设计质量课程设计选题合理，计算过程简练准确，分析问题思路清晰，结构严谨，文理通顺，撰写规范，图表完备正确。453创新工作中有创新意识，对前人工作有一些改进或有一定应用价值。54答辩能正确回答指导教师所提出的问题。30总分评语：指导教师：2021年7月19日

论大学生写作能力写作能力是对自己所积累的信息进行选择、提取、加工、改造并将之形成为书面文字的能力。积累是