專利名稱：一種高分辨率巖性框架繪制方法
技術(shù)領(lǐng)域：
本發(fā)明屬于油氣地質(zhì)勘探技術(shù)領(lǐng)域，特別涉及一種高分辨率巖性框架繪制方法，基于沉積粒度變化和反映沉積動(dòng)力學(xué)特征，實(shí)現(xiàn)能在巖性剖面上反映更多沉積學(xué)信息的高分辨率巖性框架繪制方法和技術(shù)。
背景技術(shù)：
在油氣地質(zhì)勘探研究中，柱狀圖是最基本的圖件之一。隨著計(jì)算機(jī)技術(shù)的飛速發(fā)展，手工繪圖正逐步被計(jì)算機(jī)繪圖所取代。目前，繪圖軟件所繪制的巖性剖面中將所有巖性均用矩形框架來(lái)表示。對(duì)于陸源碎屑巖來(lái)說(shuō)，這樣既不能表示出其粒度變化情況，又不能反映出它形成時(shí)的沉積動(dòng)力學(xué)特征。在層序地層學(xué)，尤其是高分辨率層序地層學(xué)研究中，更低級(jí)次基準(zhǔn)面旋回層序的劃分和精細(xì)的沉積微相研究，對(duì)柱狀圖中巖性剖面的繪制提出了更高的要求，要求在巖性剖面上反映更多的沉積學(xué)信息。實(shí)際上，陸源碎屑巖的粒度一般都不是均一的，而是漸變的，大部分均具有向上變細(xì)的正粒序，如曲流河邊灘砂巖和辮狀河心灘砂礫巖等；少部分具有向上變粗的逆粒序，如決口扇和河口壩粉砂巖-細(xì)砂巖等。

發(fā)明內(nèi)容
為了克服上述現(xiàn)有技術(shù)的不足，本發(fā)明的目的采用一種高分辨率巖性框架繪制方法，采用曲邊矩形代替矩形來(lái)繪制陸源碎屑巖巖性框架，既能表示出其粒度變化情況，又能反映出它形成時(shí)的沉積動(dòng)力學(xué)特征。為了實(shí)現(xiàn)上述目的，本發(fā)明采用的技術(shù)和方案是一種高分辨率巖性框架繪制方法，包括如下步驟
第一步，陸源碎屑巖的粒度越粗，其堆積速率越快，粒度變化也越明顯，所以，可按陸源碎屑巖粒度粗細(xì)分配其曲邊矩形框架中曲邊的寬度定義粉砂巖、細(xì)(中、粗)砂巖、礫巖框架總寬度，其中，粉砂巖寬度為3. 0mm，細(xì)砂巖、中砂巖和粗砂巖為4. 0mm，礫巖為5. Omm ；
第二步，曲邊矩形框架中的曲邊可看作由均分曲邊高度和寬度的曲邊右上角曲線①和曲邊右下角曲線②兩部分組成，兩曲線的函數(shù)表達(dá)式可由基本初等三角函數(shù)經(jīng)適當(dāng)?shù)淖儞Q而生成，建立二維XOY坐標(biāo)系，繪制曲邊矩形框架中右上角曲線①；首先由函數(shù)Y=C0s(x)，X G (0，/2)開始變換；其次，該函數(shù)應(yīng)滿足曲邊矩形框架中右上角曲線①的高度和寬度所限定的條件振幅等于曲線①高度的2倍，周期等于曲線①寬度的4倍，即A=b、T=2 / a ;則變換為y=bcos ( a x), x G (0, /2);再沿y軸向下平移b,則變換為曲線①y=bcos (ax)_b,xG(0, n /2)；
第三步，建立二維X1O1Y1坐標(biāo)系，繪制曲邊矩形框架中右下角曲線②；由函數(shù)yi=sin(X1)jX1 g(0, /2)開始變換，同樣，函數(shù)也應(yīng)滿足曲邊矩形中右下角曲線②高度和寬度所限定的條件振幅等于曲線①高度的2倍，周期等于曲線①寬度的4倍，A=Ji / a，T=4b;則變換為71=^288111((31/213) X1), X1 G (0,b));至此，曲邊矩形框架中右下角曲線②的繪制函數(shù)表達(dá)式已生成；
第四步，曲邊矩形中右下角曲線②和右上角曲線①的函數(shù)表達(dá)式是在不同坐標(biāo)系下生成的，為了將兩曲線統(tǒng)一到XOY坐標(biāo)系中，應(yīng)對(duì)坐標(biāo)系X1O1Y1進(jìn)行變換；首先，X1軸沿Y1軸向左平移/2a,則曲邊右下角曲線②的函數(shù)表達(dá)式變?yōu)閅1= Ji /2asin(Ji /2bxl) + Ji /2a;X1 G (0，b);其次，Y1軸沿X1軸向上平移b，則曲邊右下角曲線②的函數(shù)表達(dá)式變?yōu)閅1=H/2asinO/2b(xl-b)) + Ji/2a; X1 G (b, 2b);最后，將坐標(biāo)系X1O1Y1逆時(shí)針旋轉(zhuǎn)/2與XOY坐標(biāo)系重合，則曲邊右下角曲線②的函數(shù)表達(dá)式變?yōu)閥=-2b/ n arcsin (2a / n x~l) -b, x G ( n /2a, n/a),至此，曲邊右下角曲線②和曲邊右上角曲線①的函數(shù)表達(dá)式在XOY坐標(biāo)系中銜接起來(lái)；
第五步，在X軸上將曲邊30等分，共31個(gè)點(diǎn)，曲邊上各點(diǎn)的橫坐標(biāo)，曲邊寬度(mm)與其巖性代碼相等，并將單位換算為0.1 mm ;
第六步，粒度變化為正粒序的情況，采用方法如下
曲邊矩形右上角曲線繪制 y=bcos(ax)~b, x G (0, 3i 2a)，
曲邊矩形右下角曲線繪制 y=-2b Ti arcsin (2a n x_l) ~h, x G (3i 2a, a)；
第七步粒度變化為逆粒序的情況，采用技術(shù)方法如下
曲邊矩形右上角曲線繪制 y=2b Ti arcsin (2a n (x_l) _b x G (Ji 2a, a)
曲邊矩形右下角曲線繪制 y=-bcos(ax)_b; x G (0, 3i 2a)；
第八步，粒度從下往上沒有變化情況，直接用矩形框；
第九步，VC++程序中，F(xiàn)RAME是視圖類中定義的巖性曲邊矩形函數(shù)，其主要將第一步到第八步驟的方法用于生成曲邊框架上各點(diǎn)的坐標(biāo)；OnDraw是視圖類的繪圖函數(shù)，通過調(diào)用繪圖函數(shù)OnDraw內(nèi)部的Polygon函數(shù)進(jìn)行巖性框架的繪制高分辨率巖性框架。本發(fā)明的有益效果是
本發(fā)明一種高分辨率巖性框架繪制方法既能表示出其粒度變化情況，又能反映出它形成時(shí)的沉積動(dòng)力學(xué)特征。它將傳統(tǒng)的地學(xué)和現(xiàn)代的計(jì)算機(jī)結(jié)合起來(lái)，不僅可以把科研人員從繁瑣、繁重的室內(nèi)基礎(chǔ)工作中解放出來(lái)，而且還能夠創(chuàng)造可觀的經(jīng)濟(jì)效益，從而促進(jìn)地學(xué)的進(jìn)一步發(fā)展。


附圖1為本發(fā)明采用一種高分辨率巖性框架繪制技術(shù)方法處理模型圖，其中附圖1中圖1(a)為正粒序中，曲邊矩形框架中均分曲邊高度和寬度的曲邊右上角曲線①和曲邊右下角曲線②兩部分組成；(b) XOY坐標(biāo)系中，正粒序情況中，繪制曲邊矩形框架中右上角曲線①，(c) X1O1Y1坐標(biāo)系中，正粒序情況中，繪制曲邊矩形框架中右下角曲線②，(d)正粒序情況中，X1軸沿Y1軸向左平移/2a，曲邊矩形框架中右下角曲線②繪制圖，(e)正粒序情況中，Y1軸沿X1軸向上平移b，曲邊矩形框架中右下角曲線②繪制圖，(f)逆粒序中，曲邊右上角曲線①和曲邊右下角曲線②繪制。圖2a為本發(fā)明實(shí)施例一套8層不同粒度巖性段的高分辨率巖性柱狀效果圖；圖2b為地層古沉積環(huán)境水動(dòng)力特征。具體實(shí)施方法
下面結(jié)合附圖對(duì)本發(fā)明進(jìn)一步詳細(xì)說(shuō)明。
參見附圖1，一種高分辨率巖性框架繪制方法，包括如下步驟
第一步，陸源碎屑巖的粒度越粗，其堆積速率越快，粒度變化也越明顯，所以，可按陸源碎屑巖粒度粗細(xì)分配其曲邊矩形框架中曲邊的寬度定義粉砂巖、細(xì)(中、粗)砂巖、礫巖框架總寬度，其中，粉砂巖寬度為3. 0mm，細(xì)砂巖、中砂巖和粗砂巖為4. 0mm，礫巖為5. Omm ；
第二步，曲邊矩形框架中的曲邊可看作由均分曲邊高度和寬度的曲邊右上角曲線①和曲邊右下角曲線②兩部分組成，兩曲線的函數(shù)表達(dá)式可由基本初等三角函數(shù)經(jīng)適當(dāng)?shù)淖儞Q而生成。建立二維XOY坐標(biāo)系，繪制曲邊矩形框架中右上角曲線①(附圖1中a、b)，由函數(shù)Y=COS (x),x G (0，/2)開始變換，該函數(shù)應(yīng)滿足曲邊矩形框架中右上角曲線①的高度和寬度所限定的條件振幅等于曲線①高度的2倍，周期等于曲線①寬度的4倍，即A=b、T=2 / a ;則變換為y=bcos ( a x), x G (0, /2);再沿y軸向下平移b,則變換為曲線①y=bcos (ax)_b,xG(0, n /2)；
第三步，建立二維X1O1Y1坐標(biāo)系，繪制曲邊矩形框架中右下角曲線②(附圖1中c);由函數(shù)yi=sin Cx1), X1 G(0, /2)開始變換；同樣，函數(shù)也應(yīng)滿足曲邊矩形中右下角曲線②高度和寬度所限定的條件振幅等于曲線①高度的2倍，周期等于曲線①寬度的4倍，A= Ji/a，T=4b;則變換為 Y1=Ji/2asin(0/2b) X1), X1 G (0，b)(圖1 (c));至此，曲邊矩形框架中右下角曲線②的繪制函數(shù)表達(dá)式已生成；
第四步，曲邊矩形右下角曲線②和右上角曲線①的函數(shù)表達(dá)式是在不同坐標(biāo)系下生成的，為了將兩曲線統(tǒng)一到XOY坐標(biāo)系中，應(yīng)對(duì)坐標(biāo)系X1O1Y1進(jìn)行變換；首先，X1軸沿Y1軸向左平移/2a (附圖1中d)，則曲邊右下角曲線②的函數(shù)表達(dá)式變?yōu)閅1=JIZ^asin(JiAbxI)+ Ji/2a; X1 G (0，b);其次，Y1 軸沿 X1 軸向上平移 b (附圖1 中 e)，則曲邊右下角曲線②的函數(shù)表達(dá)式變?yōu)閅1=H/2asin O/2b (xl-b) ) + 3i/2a; X1 G (b, 2b)；最后，將坐標(biāo)系X1O1Y1逆時(shí)針旋轉(zhuǎn)/2與XOY坐標(biāo)系重合，則曲邊右下角曲線②的函數(shù)表達(dá)式變?yōu)閥= -2b/arcsin(2a / x-l)_b, x G O/2a, n/a),至此，曲邊右下角曲線②和曲邊右上角曲線①的函數(shù)表達(dá)式在XOY坐標(biāo)系中銜接起來(lái)；
第五步，在X軸上將曲邊30等分，共31個(gè)點(diǎn)，曲邊上各點(diǎn)的橫坐標(biāo)，曲邊寬度(mm)與其巖性代碼相等，并將單位換算為0.1 mm ;
第六步，粒度變化為正粒序的情況，采用技術(shù)方法如下
曲邊矩形右上角曲線繪制(附圖1中b中曲線①) y=bcos(ax) _bX G (0, 31 2a)
曲邊矩形右下角曲線繪制(附圖1中b中曲線②) y=-2b Ti arcsin (2a n x_l) _b; x G (Ji 2a, a)
第七步，粒度變化為逆粒序的情況，采用技術(shù)方法如下
曲邊矩形右上角曲線繪制(附圖1中f中曲線①) y=2b arcsin (2a (x_l)_b, x G (Ji 2a, a)
曲邊矩形右下角曲線繪制(附圖1中f中曲線②) y=-bcos(ax)~b, x G (0, 3i 2a)
第八步，粒度從下往上沒有變化情況，直接用矩形框；
第九步，VC++程序中，F(xiàn)RAME是視圖類中定義的巖性曲邊矩形函數(shù)，其主要將第一步到第八步驟的方法用于生成曲邊框架上各點(diǎn)的坐標(biāo)；OnDraw是視圖類的繪圖函數(shù)，通過調(diào)用繪圖函數(shù)OnDraw內(nèi)部的Polygon函數(shù)進(jìn)行巖性框架的繪制高分辨率巖性框架。實(shí)施例
在實(shí)際應(yīng)用中，可將表示陸源碎屑巖巖性框架的曲邊矩形做成一個(gè)函數(shù)，在曲邊上提取31點(diǎn)，與曲邊矩形左邊框兩端點(diǎn)構(gòu)成一個(gè)封閉的曲邊矩形，然后用繪圖函數(shù)進(jìn)行繪制。在VC++源程序中，F(xiàn)RAME是視圖類中定義的曲邊矩形函數(shù)，用于生成曲邊矩框架上各點(diǎn)的坐標(biāo)；OnDraw是視圖類中的繪圖函數(shù)，通過調(diào)用繪圖函數(shù)Polygon進(jìn)行巖性框架的繪制。FRAME函數(shù)各形參的意義為thichness為所繪制巖性的厚度；rock為該巖性代碼；grade為該巖性的粒度變化情況；origin為表示該巖性框架的曲邊矩形左上角頂點(diǎn)的坐標(biāo)。需要說(shuō)明的是，程序設(shè)定的映像模式為MM_L0METRIC，即窗口上一個(gè)邏輯單位映像到視口上為0.1mm,且視口坐標(biāo)系的Y軸向上，X軸向右；存放曲邊矩形框架上各點(diǎn)坐標(biāo)的POINT型一維數(shù)組變量sandframe (有33個(gè)元素)定義為public型，以便于函數(shù)間數(shù)據(jù)共享。以一套共計(jì)8個(gè)層位不同粒度巖性段為例，其中，I套逆粒序地層，I套粒序無(wú)變化地層和6套正粒序地層，巖性有礫巖、細(xì)(中、粗)砂巖和泥巖。第一步，陸源碎屑巖的粒度越粗，其堆積速率越快，粒度變化也越明顯，所以，可按陸源碎屑巖粒度粗細(xì)分配其曲邊矩形框架中曲邊的寬度定義粉砂巖、細(xì)(中、粗)砂巖、礫巖框架總寬度，其中，粉砂巖寬度為3. Omm,細(xì)砂巖、中砂巖和粗砂巖為4. Omm,礫巖為
5.Omm ；
第二步，曲邊矩形框架中的曲邊可看作由均分曲邊高度和寬度的曲邊右上角曲線①和曲邊右下角曲線②兩部分組成，兩曲線的函數(shù)表達(dá)式可由基本初等三角函數(shù)經(jīng)適當(dāng)?shù)淖儞Q而生成。建立二維XOY坐標(biāo)系，繪制曲邊矩形框架中右上角曲線①(附圖1中a、b);由函數(shù)Y=Cos (x), X G (0, /2)開始變換；該函數(shù)應(yīng)滿足曲邊矩形框架中右上角曲線①的高度和寬度所限定的條件振幅等于曲線①高度的2倍，周期等于曲線①寬度的4倍，即A=b、T=2 / a ;則變換為y=bcos ( a x), x G (0, /2);再沿y軸向下平移b,則變換為曲線①y=bcos (ax)_b,xG(0, n /2)；
第三步建立二維X1O1Y1坐標(biāo)系，繪制曲邊矩形框架中右下角曲線②(附圖1中c);由函數(shù)yi=sin Cx1), X1 G(0, /2)開始變換；同樣，函數(shù)也應(yīng)滿足曲邊矩形中右下角曲線②高度和寬度所限定的條件振幅等于曲線①高度的2倍，周期等于曲線①寬度的4倍，A= Ji/a，T=4b ;貝丨』變換為 Y1= /2asin(0/2b) X1), X1 G (O，b)(圖1 (c));至此，曲邊矩形框架中右下角曲線②的繪制函數(shù)表達(dá)式已生成；
第四步，曲邊矩形右下角曲線②和右上角曲線①的函數(shù)表達(dá)式是在不同坐標(biāo)系下生成的；為了將兩曲線統(tǒng)一到XOY坐標(biāo)系中，應(yīng)對(duì)坐標(biāo)系X1O1Y1進(jìn)行變換；首先，X1軸沿Y1軸向左平移/2a (圖1 (d))，則曲邊右下角曲線②的函數(shù)表達(dá)式變?yōu)閅1=JiAasin(JiAbxi)+ Ji/2a;X1 G (0，b);其次，Y1 軸沿 X1 軸向上平移 b (圖1 (e)中)，則曲邊右下角曲線②的函數(shù)表達(dá)式變?yōu)閅1=H/2asin O/2b (xl-b)) +/2a; X1 G (b, 2b)；最后，將坐標(biāo)系X1O1Y1逆時(shí)針旋轉(zhuǎn)/2與XOY坐標(biāo)系重合，則曲邊右下角曲線②的函數(shù)表達(dá)式變?yōu)閥= -2b/arcsin(2a / x-l)_b, x G O/2a, n/a),至此，曲邊右下角曲線②和曲邊右上角曲線①的函數(shù)表達(dá)式在XOY坐標(biāo)系中銜接起來(lái)；
第五步，在X軸上將曲邊30等分，共31個(gè)點(diǎn)，曲邊上各點(diǎn)的橫坐標(biāo)，曲邊寬度(mm)與其巖性代碼相等，并將單位換算為0.1 mm ;
第六步，粒度變化為正粒序的情況，采用方法如下
曲邊矩形右上角曲線繪制(附圖1中b中曲線①) y=bcos(ax) _b x G (0, 3i 2a)
曲邊矩形右下角曲線繪制(附圖1中b中曲線②) y=-2b Ti arcsin (2a n x_l) _b; x G (Ji 2a, a)
第七步，粒度變化為逆粒序的情況，采用方法如下
曲邊矩形右上角曲線繪制(附圖1中f中曲線①) y=2b Ti arcsin (2a n (x_l) _b; x G (Ji 2a, a)
曲邊矩形右下角曲線繪制(附圖1中f中曲線②) y=-bcos(ax)_b; x G (0, 3i 2a)；
第八步，粒度從下往上沒有變化情況，直接用矩形框；
第九步，VC++程序中，F(xiàn)RAME是視圖類中定義的巖性曲邊矩形函數(shù)，其主要將第一步到第八步驟的方法用于生成曲邊框架上各點(diǎn)的坐標(biāo)；OnDraw是視圖類的繪圖函數(shù)，通過調(diào)用繪圖函數(shù)OnDraw內(nèi)部的Polygon函數(shù)進(jìn)行巖性框架的繪制高分辨率巖性框架。圖2是一套8個(gè)層位不同粒度巖性段的高分辨率巖性柱狀效果圖。較好地表示出了砂礫巖的粒度變化情況，同時(shí)，也反映出了其形成時(shí)的沉積動(dòng)力學(xué)特征，從下往上，古沉積環(huán)境從辮狀河到泛濫平原水動(dòng)力條件由強(qiáng)變?nèi)酰簽E平原中的決口扇是在水動(dòng)力突然增大的條件下形成的。
權(quán)利要求
1.一種高分辨率巖性框架繪制方法，其特征在于，包括如下步驟 第一步，陸源碎屑巖的粒度越粗，其堆積速率越快，粒度變化也越明顯，所以，可按陸源碎屑巖粒度粗細(xì)分配其曲邊矩形框架中曲邊的寬度定義粉砂巖、細(xì)(中、粗)砂巖、礫巖框架總寬度，其中，粉砂巖寬度為3. Omm,細(xì)砂巖、中砂巖和粗砂巖為4. Omm,礫巖為.5.Omm ； 第二步，曲邊矩形框架中的曲邊由均分曲邊高度和寬度的曲邊右上角曲線①和曲邊右下角曲線②兩部分組成，兩曲線的函數(shù)表達(dá)式可由基本初等三角函數(shù)經(jīng)適當(dāng)?shù)淖儞Q而生成，建立二維XOY坐標(biāo)系，繪制曲邊矩形框架中右上角曲線①；首先由函數(shù)Y=COS (X)，X e (O, /2)開始變換；其次，該函數(shù)應(yīng)滿足曲邊矩形框架中右上角曲線①的高度和寬度所限定的條件振幅等于曲線①高度的2倍，周期等于曲線①寬度的4倍，即A=b、Τ=2 π /α ;則變換為y=bcos ( α χ), X e (O, π/2);再沿y軸向下平移b,則變換為曲線①y=bcos(ax)-b,xe(0, π /2)； 第三步，建立二維X1O1Y1坐標(biāo)系，繪制曲邊矩形框架中右下角曲線②；由函數(shù)yi=sin(X1)jX1 e(0, π/2)開始變換，同樣，函數(shù)也應(yīng)滿足曲邊矩形中右下角曲線②高度和寬度所限定的條件振幅等于曲線①高度的2倍，周期等于曲線①寬度的4倍，A= π / a，T =4b;則變換為Y1= π/2asin(0/2b) X1), X1 e (0, b);至此，曲邊矩形框架中右下角曲線②的繪制函數(shù)表達(dá)式已生成； 第四步，曲邊矩形中右下角曲線②和右上角曲線①的函數(shù)表達(dá)式是在不同坐標(biāo)系下生成的，為了將兩曲線統(tǒng)一到XOY坐標(biāo)系中，應(yīng)對(duì)坐標(biāo)系X1O1Y1進(jìn)行變換；首先，X1軸沿Y1軸向左平移π /2a,則曲邊右下角曲線②的函數(shù)表達(dá)式變?yōu)閅1= Ji /2asin(Ji /2bxl) + Ji /2a;X1 e (O, b);其次，Y1軸沿X1軸向上平移b，則曲邊右下角曲線②的函數(shù)表達(dá)式變?yōu)閅1=H/2asinO/2b(xl-b)) + Ji/2a; X1 e (b, 2b);最后，將坐標(biāo)系X1O1Y1逆時(shí)針旋轉(zhuǎn)π/2與XOY坐標(biāo)系重合，則曲邊右下角曲線②的函數(shù)表達(dá)式變?yōu)閥=-2b/ narcsin (2a / η χ-l) -b, χ e ( ji /2a, π /a)； 至此，曲邊右下角曲線②和曲邊右上角曲線①的函數(shù)表達(dá)式在XOY坐標(biāo)系中銜接起來(lái)； 第五步，在X軸上將曲邊30等分，共31個(gè)點(diǎn)，曲邊上各點(diǎn)的橫坐標(biāo)，曲邊寬度(mm)與其巖性代碼相等，并將單位換算為O.1 mm ; 第六步，粒度變化為正粒序的情況，采用方法如下 曲邊矩形右上角曲線繪制 y=bcos(ax)~b, X e (O, Ji 2a)， 曲邊矩形右下角曲線繪制 y=-2b τι arcsin (2a η χ-l) -b, X e ( ji 2a, π a)； 第七步粒度變化為逆粒序的情況，采用技術(shù)方法如下 曲邊矩形右上角曲線繪制 y=2b τι arcsin (2a η (χ-l) -b X e ( ji 2a, n a)曲邊矩形右下角曲線繪制 y=-bcos(ax)_b;X e (O, Ji 2a)； 第八步，粒度從下往上沒有變化情況，直接用矩形框； 第九步，VC++程序中，F(xiàn)RAME是視圖類中定義的巖性曲邊矩形函數(shù)，其主要將第一步到第八步驟的方法用于生成曲邊框架上各點(diǎn)的坐標(biāo)；OnDraw是視圖類的繪圖函數(shù)，通過調(diào)用繪圖函數(shù)OnDraw內(nèi)部的Polygon函數(shù)進(jìn)行巖性框架的繪制高分辨率巖性框架。
2.根據(jù)權(quán)利要求1所述的，其特征在于，包括以下步驟第一步，陸源碎屑巖的粒度越粗，其堆積速率越快，粒度變化也越明顯，所以，可按陸源碎屑巖粒度粗細(xì)分配其曲邊矩形框架中曲邊的寬度定義粉砂巖、細(xì)(中、粗)砂巖、礫巖框架總寬度，其中，粉砂巖寬度為3. Omm,細(xì)砂巖、中砂巖和粗砂巖為4. Omm, 5樂巖為5. Omm ； 第二步，曲邊矩形框架中的曲邊可看作由均分曲邊高度和寬度的曲邊右上角曲線①和曲邊右下角曲線②兩部分組成，兩曲線的函數(shù)表達(dá)式可由基本初等三角函數(shù)經(jīng)適當(dāng)?shù)淖儞Q而生成,建立二維XOY坐標(biāo)系，繪制曲邊矩形框架中右上角曲線①；由函數(shù)Y=cos(x)，X e(O, π/2)開始變換；該函數(shù)應(yīng)滿足曲邊矩形框架中右上角曲線①的高度和寬度所限定的條件振幅等于曲線①高度的2倍，周期等于曲線①寬度的4倍，即A=b、Τ=2π/α ;則變換為y=bcos ( α X)，X e (O, Ji/2);再沿y軸向下平移b,則變換為曲線①y=bcos ( α χ)-b,X e (O, Ji /2)； 第三步，建立二維X1O1Y1坐標(biāo)系，繪制曲邊矩形框架中右下角曲線②；由函數(shù)yi=sin(X1)jX1 e(0, π/2)開始變換；同樣，函數(shù)也應(yīng)滿足曲邊矩形中右下角曲線②高度和寬度所限定的條件振幅等于曲線①高度的2倍，周期等于曲線①寬度的4倍，A= π / a，T=4b ;則變換為Y1=JiAasin ((Ji/2b) X1), X1 e (O, b);至此，曲邊矩形框架中右下角曲線②的繪制函數(shù)表達(dá)式已生成； 第四步，曲邊矩形右下角曲線②和右上角曲線①的函數(shù)表達(dá)式是在不同坐標(biāo)系下生成的；為了將兩曲線統(tǒng)一到XOY坐標(biāo)系中，應(yīng)對(duì)坐標(biāo)系X1O1Y1進(jìn)行變換；首先，X1軸沿Y1軸向左平移π /2a,則曲邊右下角曲線②的函數(shù)表達(dá)式變?yōu)閥^Ti^asinCTi^bx^ + Ti/2a-,X1 e (O, b);其次，Y1軸沿X1軸向上平移b，則曲邊右下角曲線②的函數(shù)表達(dá)式變?yōu)閅1=H/2asinO/2b(xl-b)) + Ji/2a; X1 e (b, 2b);最后，將坐標(biāo)系X1O1Y1逆時(shí)針旋轉(zhuǎn)π/2與XOY坐標(biāo)系重合，則曲邊右下角曲線②的函數(shù)表達(dá)式變?yōu)閥=-2b/ η arcsin (2a / η χ-l) -b, x e ( ji /2a, π/a),至此，曲邊右下角曲線②和曲邊右上角曲線①的函數(shù)表達(dá)式在XOY坐標(biāo)系中銜接起來(lái)； 第五步，在X軸上將曲邊30等分，共31個(gè)點(diǎn)，曲邊上各點(diǎn)的橫坐標(biāo)，曲邊寬度mm與其巖性代碼相等，并將單位換算為O.1mm ； 第六步，粒度變化為正粒序的情況，采用方法如下 曲邊矩形右上角曲線繪制 y=bcos(ax) _b χ e (O, Ji 2a) 曲邊矩形右下角曲線繪制 y=-2b τι arcsin (2a η χ-l) -b;X e (JI 2a, n a) 第七步，粒度變化為逆粒序的情況，采用方法如下 曲邊矩形右上角曲線繪制 y=2b τι arcsin (2a η (χ-l) -b; χ e ( ji 2a, π a) 曲邊矩形右下角曲線繪制 y=-bcos(ax)_b; χ e (O, Ji 2a)； 第八步，粒度從下往上沒有變化情況，直接用矩形框； 第九步，VC++程序中，F(xiàn)RAME是視圖類中定義的巖性曲邊矩形函數(shù)，其將第一步到第八步驟的方法用于生成曲邊框架上各點(diǎn)的坐標(biāo)；OnDraw是視圖類的繪圖函數(shù)，通過調(diào)用繪圖 函數(shù)OnDraw內(nèi)部的Polygon函數(shù)進(jìn)行巖性框架的繪制高分辨率巖性框架。
全文摘要
一種高分辨率巖性框架繪制方法，包括如下步驟第一，分配其曲邊矩形框架中曲邊的寬度，定義粉砂巖、細(xì)砂巖、礫巖框架總寬度；第二，建立二維X0Y坐標(biāo)系，繪制曲邊矩形框架中右上角曲線①；第三，建立二維X101Y1坐標(biāo)系，繪制曲邊矩形框架中右下角曲線②；第四，將兩曲線統(tǒng)一到XOY坐標(biāo)系中,應(yīng)對(duì)坐標(biāo)系X1O1Y1進(jìn)行變換；第五，在X軸上將曲邊30等分,共31個(gè)點(diǎn)；第六，粒度變化為正粒序的情況時(shí)的，繪制曲邊矩形右上角和右下角的曲線；第七，粒度變化為逆粒序的情況，繪制曲邊矩形右上角、右下角的曲線；第八，粒度從下往上沒有變化，直接用矩形框；第九，將第一步到第八步驟生成曲邊框架上各點(diǎn)的坐標(biāo)，調(diào)用繪圖函數(shù)OnDraw內(nèi)部的Polygon函數(shù)繪制高分辨率巖性框架。
文檔編號(hào)G06F17/50GK103020372SQ20121056429
公開日2013年4月3日 申請(qǐng)日期2012年12月24日 優(yōu)先權(quán)日2012年12月24日
發(fā)明者魏欽廉, 肖玲, 高紅燦 申請(qǐng)人:西安石油大學(xué)