測(cè)井曲線自動(dòng)分層問(wèn)題

129431、長(zhǎng)、長(zhǎng)、長(zhǎng)、長(zhǎng)、長(zhǎng)、長(zhǎng)、長(zhǎng)、長(zhǎng)、長(zhǎng)、長(zhǎng)738182、規(guī)律，依次定名為長(zhǎng)、長(zhǎng)、長(zhǎng)、長(zhǎng)、長(zhǎng)、長(zhǎng)、長(zhǎng)、長(zhǎng)、長(zhǎng)、長(zhǎng)31、長(zhǎng)32層。通常這些工作都是通過(guò)人工來(lái)進(jìn)行的，這就是所謂人工分層方法。該方(DEN)(AC)(CNL)(GR)、自然(SP)(RT)171813

1DEVi(井斜)AZImCAL(井徑)PORWPORT(總孔隙度)POR(孔隙度)PORRPORFPERM(絕對(duì)滲透率)SWSH(泥質(zhì)含量)SXOPOW(含水孔隙度)CARBFW(產(chǎn)水率)BULKCALC

CL(粘土體積)PORXPORH(油氣重量)DEPTHDENRILDRILMR4.0(4MSPGRAC(聲波測(cè)井)RMLCNLRT(電阻率測(cè)井)WARMFA(視泥漿電阻率)（一）C=a1b1+a2b2+a3b3+a4b4+a5b5+......+a67CBACBA

C（分層；下面介紹 方法11號(hào)井為標(biāo)準(zhǔn)井，根據(jù)此井的各種測(cè)井曲線數(shù)據(jù)，一號(hào)井分為167CURVENAME=DEN,RILD,RILM,RLL8,SP1,R4.0,SP,DEVi,AZIm,GR,AC,RML,RMN,CAL,CNL,R4.0%,cnl%,GR%,AC%,RILD%,RILM%,RLL8%,SP1%,SP%,den%,CAL%,RML%,RMN%,PORW,PORT,PORR,PORF,PERM,SW,SH,SXO,POR,SPC,RT,PW,POW,PO,YO,YOW,YW,CARB,CENG,FO,FOW,FW,CID,CIM,CIL8,THK,PORD,PORC,PP,BULK,HF,PF,CALC,CL,PORX,PORH,RWA,RMFA層--深度關(guān)系（17當(dāng)然也可以取組1311

1（1）67形式為15行67 的矩陣，形式簡(jiǎn)寫(xiě)如 25B

45 1CT0=（0,1,2,3,4,5,6,7,8,9,10,11,12,13,14）0--143192，當(dāng)然完全可以把C0

01 1 C233 BA即c1 2bba1b2,67 2c322A

k1 k2 k3

A0

a1a2a3 c1c2bbb2,67 2c3

CBA0BC2。origin

C0—150.5-0.50.5-1.5但實(shí)際上受一些因素的影響（A0到0.51CC12、C1C1+kd（k=1,2,...14)所對(duì)應(yīng)的深度值作為以后各層的C1，C1=0；0-12，1-23C0的取值決定的。針對(duì)4號(hào)井進(jìn) yy層函數(shù)---x深度4yy層函數(shù)值x自動(dòng)分層數(shù)據(jù)與人工分層數(shù)據(jù)對(duì)比自動(dòng)分層數(shù) 人工分層數(shù)長(zhǎng)31長(zhǎng)31長(zhǎng)32長(zhǎng)32長(zhǎng)33長(zhǎng)33長(zhǎng)41長(zhǎng)41長(zhǎng)42長(zhǎng)42長(zhǎng)61長(zhǎng)61長(zhǎng)62長(zhǎng)62長(zhǎng)63長(zhǎng)63長(zhǎng)71長(zhǎng)71長(zhǎng)72長(zhǎng)72長(zhǎng)73長(zhǎng)73長(zhǎng)81長(zhǎng)81長(zhǎng)82長(zhǎng)82長(zhǎng)91長(zhǎng)91長(zhǎng)92長(zhǎng)92excelExcelA0散點(diǎn)插值求解程序：294330368410 575697 814257290326369414458 612682717765213251 422455 578621663700740 276311 398438483518 680720760803 492 618 741780 900.4934972 334368410454495528576615652690770810].]forforj=1:1:6

XYXY長(zhǎng)31底 長(zhǎng)63底深的分層結(jié)果（其中8-13通過(guò)上面插值獲得長(zhǎng)長(zhǎng)長(zhǎng)長(zhǎng)長(zhǎng)長(zhǎng)長(zhǎng)長(zhǎng)12 456789長(zhǎng)71底 長(zhǎng)92底深分層結(jié)長(zhǎng)長(zhǎng)長(zhǎng)長(zhǎng)長(zhǎng)長(zhǎng)長(zhǎng)123456789層--深度關(guān)系（8131針 號(hào)井進(jìn)行自動(dòng)分層處自動(dòng)分層數(shù) 插值分層數(shù)長(zhǎng)31長(zhǎng)31長(zhǎng)32長(zhǎng)32長(zhǎng)33長(zhǎng)33長(zhǎng)41長(zhǎng)41長(zhǎng)42長(zhǎng)42長(zhǎng)61長(zhǎng)61長(zhǎng)62長(zhǎng)62長(zhǎng)63長(zhǎng)63長(zhǎng)71長(zhǎng)71長(zhǎng)72長(zhǎng)72長(zhǎng)73長(zhǎng)73長(zhǎng)81長(zhǎng)81長(zhǎng)82長(zhǎng)82長(zhǎng)91長(zhǎng)91長(zhǎng)92長(zhǎng)9215解時(shí)，810

44YYAxis0--XAxis55YYAxis50-- XAxis該模型分層時(shí)的錯(cuò)誤 很大一部分是有一下因素造成的 不同井擁有的測(cè)井量不同，如果特 選取不當(dāng)，對(duì)應(yīng)的數(shù)值在其他井測(cè)井

特解時(shí)帶 的人工分層數(shù)據(jù)值，

A07.0使 機(jī)械工數(shù)學(xué)模型 ：高等教 數(shù)學(xué)建模與數(shù)學(xué)實(shí)驗(yàn) ：高等教育，楊雄文測(cè)井曲線自動(dòng)分層模型設(shè)計(jì)與實(shí)現(xiàn)A04040#include<iostream>#include<cmath>#definerealfloat;#defineNZ67usingnamespacestd;doublejun(floata[][7000],intj,intstart,int{inti;doubleadd=0;return(add/(end-}intgetsp(char*s,char{inti;}int{charstr[1000];FILE*p=fopen("1.txt","r");floatfscanf(p,"STDEP=%f",&h0);getc(p);floathe;fscanf(p,"ENDEP=%f",&he);getc(p);floatrlev;fscanf(p,"RLEV=%f",&rlev);getc(p);inti,j;floatarryy[NZ][7000];intsum=(he-}

{fscanf(p,"%f",&(arryy[j][i]));if(arryy[j][i]==-doublearr[100]={doubleresult[100][100];inth=0,hend=0;{result[j][i]=jun(arryy,i,h,hend);//}{fprintf(p,"%f}jun.txt15*67下面用讀入均值數(shù)據(jù)…………..]A0=[-------------------00-0000-0--]CBA0#include<iostream>#include<cmath>#definerealfloat#include<cstring>#include<cstdio>#include<stdlib.h>intNZ;charvar[67][100]={"#DEPTH","DEN","RILD","RILM","RLL8","SP1","R4.0","SP","DEVi","AZIm","GR","AC","RML","CAL","CNL","R4.0%","cnl%","AC%","RILD%","RILM%","RLL8%","SP1%","den%","CAL%","RML%","RMN%","PORW","PORR","PORF","PERM","SW","SH","SXO","POR","SPC","RT","PW","POW","PO","YO","YOW","YW","CARB","CENG","FO","FOW","FW","CID","CIM","PORD","PORC","PP","BULK","HF","PF","CALC","CL","PORX","PORH","RWA","RMFA"chars[100][100];intcmp(char*s){intn;returnn;}intsel(int{inti=cmp(s[n]);returni;}usingnamespacedoublejun(floata[][1000],intj,intstart,int{inti;doubleadd=0;return(add/(end-}intcharstr[1000];intfilen;1-13FILE*p=fopen(str,"r");floatfscanf(p,"STDEP= floathe;fscanf(p,"ENDEP= floatrlev; inti;intn=1;floatarryy[67][7500];intj;floatarr[100];intsum=0;p=fopen("shu.txt","r");//打開(kāi)A0特解文件shu.txt并特doubleresult[8000];{}{}}}內(nèi)插法畫(huà)15層圖程序 294330368410 495.3 575 697 771.8814 257290326369414458497.5 582612682717765 213251 377.9422455502578621663700740 276311 398438483518 680720760803 492 618 698741 934972 334368410454495528576615652690733770 ].]for surf(X,