王昌勇 楊寶泉 高博禹 胡曉慶 朱國金 鄭榮才

(1.成都理工大學沉積地質(zhì)研究院 成都 610059;2.中國海洋石油總公司研究總院 北京 100027)

荔灣3-1井區(qū)珠江組深水扇高分辨率層序分析及應用①

王昌勇1楊寶泉2高博禹2胡曉慶2朱國金2鄭榮才1

(1.成都理工大學沉積地質(zhì)研究院 成都 610059;2.中國海洋石油總公司研究總院 北京 100027)

以高分辨率層序分析為指導思想,以巖芯和測井分析為基礎,結(jié)合地震資料對荔灣3-1井區(qū)珠江組下部深水扇沉積層序界面和海泛面成因類型進行綜合研究,探討了深水扇短期旋回結(jié)構(gòu)類型和分布規(guī)律,建立了主要產(chǎn)層時間-地層格架,并對地層格架中的小層砂體進行劈分和等時追蹤對比,討論了不同基準面旋回條件下儲層發(fā)育特征。通過層序-巖相古地理編圖,結(jié)合GR約束反演,刻畫了不同發(fā)育階段深水扇砂體展布特征,建立了深水扇水道砂體的時空間展布模式,對今后深水扇型油藏的砂體預測具有重要借鑒意義。

深水扇 高分辨率層序 層序結(jié)構(gòu) 珠江組 白云凹陷 荔灣井區(qū)

0 引言

圖1 荔灣3-1井區(qū)構(gòu)造位置圖(據(jù)龐雄等,2007)Fig.1 Regional geology of Liwan 3-1 area

近年來,在國內(nèi)外多個地區(qū)深水扇沉積體系中發(fā)現(xiàn)了商業(yè)性油氣藏,深水扇型油氣藏逐漸成為油氣勘探的一個新領域[1～4],我國南海北部珠江口盆地白云凹陷的荔灣3-1井區(qū)珠江組下部即發(fā)育此類深水扇型油氣藏(圖1),但該深水扇研究程度較低,儲層分布規(guī)律尚處于探索階段,砂體展布預測主要依賴于地震資料解釋。高昂的鉆采費用要求砂體預測需要達到很高精度,但由于地震資料分辨率的局限性,很難達到精確描述砂體三維形態(tài)的目的,迫切需要綜合鉆井、測井資料對小層砂體進行精細刻畫,基于此目的,本次研究綜合井、震資料對荔灣3-1井區(qū)珠江組深水扇沉積進行高分辨率層序地層分析。

1 層序界面的識別

高分辨率層序地層分析中的“界面”主要包括兩種類型:其一為層序的頂?shù)捉缑?即層序界面),其二為層序內(nèi)部海平面由升到降的轉(zhuǎn)折面(即海泛面)[5]。本次研究首先根據(jù)地震剖面特征、測井曲線形態(tài)及巖芯資料,在珠江組下部地層中識別出2個對應于構(gòu)造幕式活動變化的Ⅲ類層序界面[6]及1個最大海泛面,以此作為全區(qū)宏觀層序劃分和對比依據(jù),結(jié)合測井解釋和巖芯資料,將珠江組下部深水扇沉積劃分為1個長期旋回、3個中期和5個短期旋回層序(分別相當于Vail提出的Ⅲ級、Ⅳ級和Ⅴ級層序,表1),確定產(chǎn)層為MSC1中期旋回層序。

1.1 界面測井曲線特征及識別標志

LSC1長期旋回層序界面底界(亦即珠江組與珠海組分界)位于箱型或鐘型GR曲線底部,表現(xiàn)為大型沖刷面,該界面對應于距今23.8 Ma的白云運動[7],之下一般為珠海組前陸架邊緣三角洲泥巖[8],之上為珠江組深水扇水道砂巖;LSC1長期旋回層序界面頂界位于GR曲線指狀低值與齒狀較高值轉(zhuǎn)折處,表現(xiàn)為相關(guān)整合面,該界面為深水扇與盆地相變面,之下發(fā)育低密度濁流沉積的粉砂巖,之上為大套含抱球蟲灰質(zhì)泥巖;長期旋回層序最大海泛面LMFS表現(xiàn)為測井曲線單向脈動性移動達低幅極限位置后折向增高的轉(zhuǎn)換點位置,位于大段高GR泥巖中上部(圖2)。

表1 荔灣3-1井區(qū)珠江組高分辨率層序劃分方案Table1 High-resolution sequence stratigraphic division of Zhujiang Formation of Liwan 3-1 area

1.2 界面地震剖面特征及其識別標志

根據(jù)地震合成記錄標定結(jié)果,珠江組與珠海組分界的LSC1長期旋回層序底界在地震剖面上表現(xiàn)為強連續(xù)反射第二相位,為典型標準反射層和易于連續(xù)穩(wěn)定追蹤的反射界面。位于珠江組內(nèi)部的LSC1長期旋回層序頂界在地震剖面上具有中強反射的單一相位特征,橫向連續(xù)性雖然不太好,但特征明顯,可穩(wěn)定追蹤,是典型的地震層序界面(圖3)。

2 層序基本類型及其分布規(guī)律

由于受不同成因類型的界面控制,荔灣3-1井區(qū)珠江組下部地層基準面旋回層序發(fā)育2種基本結(jié)構(gòu)類型,分別反映不同的沉積環(huán)境和水動力條件。

2.1 向上“變深”非對稱型旋回結(jié)構(gòu)(A型)

此類層序只在短期旋回層序中發(fā)育,其頂、底界面均為規(guī)模較大的沖刷面,層序內(nèi)僅保留基準面上升半旋回沉積記錄,層序主體主要為水道化的細-粗粒砂巖構(gòu)成,與粉砂巖和泥巖組成向上變細“加深”的沉積序列(圖2)。此類層序一般形成于A/S<1的過補償條件下(其中A為可容空間增量,S為沉積物通量,下同),可細分為A/S?1低可容納空間(A1型)和A/S<1的高可容納空間(A2型)2種亞類型[9]。

2.2 對稱型旋回結(jié)構(gòu)(C型)

圖2 荔灣3-1井區(qū)珠江組下部短期旋回層序的基本類型(LW3-1-D井)Fig.2 Structure types of the short term cycle sequence of the Lower Zhujiang Formation in LW3-1 area

此類層序形成于沉積物供給速率等于或小于可容納空間增長率(A/S≥1)的弱補償至欠補償條件下[9],主要由細砂巖、粉砂巖及泥巖組成向上變細復變粗的沉積旋回,為一個相對完整的海侵海退序列,在短期旋回層序中最為常見。C型短期旋回層序界面多為小型沖刷面,少量為相關(guān)整合面,由海泛面將層序劈分為基準面上升和下降兩個半旋回,按上升與下降半旋回厚度變化,可細分為以上升半旋回為主的不完全對稱型(C1型)、上升半旋回與下降半旋回近于相等的對稱型(C2型)和下降半旋回為主的不完全對稱型(C3型)3個亞類型(圖2)。

圖3 LW3-1-B、LW3-1-D井層序界面地震剖面特征Fig.3 Seismic profile with sequence boundary ofWell LW3-1-B&LW3-1-D

圖4 荔灣3-1井區(qū)珠江組深水扇中期和短期旋回層序結(jié)構(gòu)類型、疊加樣式和分布模式Fig.4 Distribution,structure and stacking pattern correlation of themiddle-short term cycle sequence of Zhujiang Formation in Liwan3-1area

2.3 不同結(jié)構(gòu)類型層序分布模式

荔灣3-1井區(qū)珠江組中,不同結(jié)構(gòu)類型的短期旋回層序分布具有特定的規(guī)律性(圖4),區(qū)域上明顯受基準面上升期有效可容納空間向物源方向溯源遷移,而下降期間則向盆地方向順源遷移的雙向變化過程和有效可容納空間變化機制控制。中新世珠江組沉積期,荔灣3-1井區(qū)位于深水扇中扇附近,其不同結(jié)構(gòu)特征的短期旋回層序分布有如下幾個基本特點:

①A型短期旋回層序主要分布在水道活動區(qū),C型短期旋回層序主要分布在水道間漫溢沉積區(qū);

②在主水道下部以砂質(zhì)碎屑流沉積為主,下切侵蝕強烈和充填快速,以發(fā)育A1型短期旋回層序結(jié)構(gòu),在分支水道或水道側(cè)翼,流體能量減弱、沉積速率降低,主要發(fā)育A2型短期旋回層序結(jié)構(gòu);

③從主水道向分支水道及水道間漫溢方向,短期旋回層序的結(jié)構(gòu)類型呈現(xiàn)出A1型-A2型-C1型-C2型-C3型變化規(guī)律(圖4),層序結(jié)構(gòu)的垂向演化同樣存在這一特征(圖2)。

深水扇沉積體系中短期旋回層序結(jié)構(gòu)以A1型、A2型、C1型及C2型為主,缺少B型層序結(jié)構(gòu)類型,反映深水扇沉積物供給通量較為充足。位于大陸斜坡之上的峽谷是沉積物堆積的主要場所,基準面上升時期大量砂質(zhì)沉積物快速堆積充填峽谷,而基準面下降期大量沉積物則被搬運至更深海域,在研究區(qū)缺失或僅留下少量沉積記錄。同時,扇體內(nèi)部短期旋回層序的垂向演化序列具有向上A/S值逐漸增大的特征,說明這一沉積時期相對海平面逐漸升高,深水扇有向岸遷移的趨勢。

3 層序地層格架

層序地層格架,系指將同時代形成的巖層有序地納入相關(guān)年代的時間-地層對比格架中,在此基礎上進行等時地層對比和描述地層疊置樣式的地層學研究方法[10],在時間-地層格架中所標定的巖層具有同步沉積演化序列,能為有利層位及有利相帶預測提供可靠的地質(zhì)模型。

3.1 等時地層格架

圖5 荔灣3-1井區(qū)珠江組下部層序等時地層格架及砂體分布Fig.5 Sequence stratigraphic framework of the Lower Zhujiang Formation in LW3-1 area and isochronous correlation of small sand-body layers

按地層基準面原理將小層砂體納入各短期時間尺度發(fā)現(xiàn):LW3-1-A及LW3-1-C井明顯缺失珠海組末期沉積記錄,在珠江組內(nèi)部LW3-1-B及LW3-1-D井各短期下降半旋回沉積記錄有不同程度的缺失(圖5),反映珠江組深水扇發(fā)育的不同時期,水道活動區(qū)不斷發(fā)生遷移,深水扇發(fā)育初期,靠近盆地一側(cè)的LW 3-1-A及LW3-1-C井附近為主水道活動區(qū),下伏珠海組晚期沉積因遭受珠江組早期低位扇強烈侵蝕而缺失,深水扇發(fā)育中、晚期,主水道活動區(qū)遷移至更靠物源的LW3-1-B及LW3-1-D井附近,導致水道中主要發(fā)育A型短期旋回層序。

3.2 地層格架中的儲層發(fā)育特征

通過層序地層對比及GR約束地震反演成果發(fā)現(xiàn),研究區(qū)珠江組深水扇沉積砂體主要分布在中期旋回層序MSC1上升半旋回相域,其中LW3-1-A及LW3-1-C井砂體主要發(fā)育在SSC1上升半旋回相域, LW3-1-B及 LW3-1-D井砂體主要發(fā)育在 SSC2和SSC3上升半旋回相域,基準面下降半旋回相域幾乎不發(fā)育砂體(圖5)。上述特征表明:

①當基準面呈初始上升狀態(tài)時,深水扇水道沉積最為活躍,以發(fā)育向盆地低部位下超、向高部位上超的主動進積序列為特征,非常有利于儲集砂體的發(fā)育;

②伴隨海平面上升幅度逐漸加大靠近最大海泛面附近、有效可容納空間向物源方向遷移,深水扇向物源方向退縮,沉積物變細,形成向上變細的退積序列,此時水道容易發(fā)生改道,泥質(zhì)隔層增多,當海平面上升至中期海泛面后,為深水扇提供物源的陸架邊緣三角洲前緣將向岸退積越過坡折,導致深水扇水道突然消失,但來自前三角洲的物源仍能以低密度濁流的形式搬運到研究區(qū);

圖6 砂體平面展布特征圖(A。波阻抗屬性圖;B。振幅屬性圖;C.GR反演砂體圖;D。珠江組下部深水扇沉積相圖及砂體等厚圖,虛線框為圖E、F范圍;E.SSC1砂體等厚圖;F.SSC2+3砂體等厚圖)Fig.6 Sandstone distribution of deepwater fan(A.wave impedance attributemap;B.amplitude attributemap;C.sand body distribution by GR inversion;D.Sedimentary distribution and isolinemap of deepwater fan of Lower Zhujiang Formation,the dashed correspond to range of figure B&C;B.Isolinemap of SSC1;C.Isolinemap of SSC2&SSC 3)

③可容納空間變化對深水扇沉積的影響極為明顯,伴隨海平面的下降和有效可容納空間向盆地方向的遷移,研究區(qū)深水扇水道中的沉積物隨重力流搬運至更深海域,此時荔灣3-1井區(qū)以南可能為更為有利的儲層發(fā)育區(qū)。

由此可見,對應基準面升降過程的可容納空間變化機制,是控制深水扇沉積體系短期基準面旋回結(jié)構(gòu)和儲層有序分布的主要因素。

4 沉積層序編圖及砂體預測

4.1 層序-沉積相編圖

以層序地層對比為依據(jù),選取深水扇水道主要發(fā)育的MSC1上升半旋回相域為編圖單元,主要根據(jù)三維地震資料解釋和GR約束地震反演成果對珠江組下部深水扇砂體進行了預測,結(jié)果顯示:砂巖發(fā)育區(qū)具有強波谷振幅和低波阻抗屬性的特征(圖6A、圖6B),地震三維工區(qū)內(nèi)砂體主要呈北西-南東向或近南北向展布,研究區(qū)北部為水道砂體主要發(fā)育區(qū),南部主要為大陸斜坡相泥巖發(fā)育區(qū)(圖6D),大部分區(qū)域砂巖累積厚度大于20 m(圖6C、圖6D)。

通過層序地層對比發(fā)現(xiàn)荔灣3-1井區(qū)珠江組下部深水扇水道主要由兩個砂體構(gòu)成,其中砂體Ⅰ發(fā)育于短期旋回SSC1之中,砂體Ⅱ則橫跨SSC2和SSC3兩個短期旋回層序,以此,本次研究選取SSC1、SSC2 +3為編圖單元,以區(qū)內(nèi)4口鉆井鉆遇厚度為依據(jù),以地震反演預測的深水扇水道形態(tài)(圖6D)為指導,結(jié)合層序地層格架和小層砂體對比成果,分別編制研究區(qū)兩個主要小層砂體的厚度等值線圖。編圖發(fā)現(xiàn):研究區(qū)水道砂體分為兩支,其中砂體Ⅰ厚度較小,一般不超過10 m,近南北向砂體厚度相對較大(圖6E);砂體Ⅱ厚度可達30 m,北西-南東向砂體最為發(fā)育(圖6F)。

4.2 砂體預測

根據(jù)砂體厚度平面分布特征,結(jié)合層序地層對比成果,對小層砂體進行刻畫和描述,并預測其時空分布規(guī)律。層序地層對比發(fā)現(xiàn):LW3-1-C及LW3-1-A井明顯缺失部分珠海組末期沉積記錄,而LW3-1-B及LW3-1-D井無明顯的地層缺失(圖5),反映LW3-1-C及LW3-1-A井位于下切能力最強的主水道附近,結(jié)合砂體厚度分布特征(圖6E),確定SSC1沉積時期LW3-1-A井位于主水道中央,LW3-1-C井位于主水道側(cè)翼,LW3-1-B及LW3-1-D井位于分支水道側(cè)緣或水道間漫溢,主水道過LW3-1-A井呈近南北向分布,砂體具有向北減薄上超和向南加厚下超的特點(圖7)。同時,通過層序?qū)Ρ劝l(fā)現(xiàn):僅LW3-1-B井SSC1短期旋回層序發(fā)育A型層序,明顯缺失基準面下降半旋回沉積記錄(圖5),結(jié)合砂體厚度分布特征(圖6F),確定SSC2～SSC3沉積時期LW3-1-B井位于主水道中央,LW3-1-D井位于主水道側(cè)翼,LW3-1-A井位于分支水道側(cè)緣,LW3-1-C井位于水道間,主水道過LW3-1-B及LW3-1-D井呈北西-南東向展布,砂體Ⅱ具有向北西上超和向西、向北逐漸增厚的特征(圖7)。

圖7 深水扇水道砂體三維形態(tài)Fig.7 3D configuration of deepwater fan channel sand body

荔灣3-1井區(qū)深水扇水道砂體時空分布的三維形態(tài)(圖7),與Vail[11]提出的深水扇沉積模式極為相似,砂體Ⅰ相當于盆底扇水道砂體,推測在LW3-1-A井以南還有較遠距離延伸,砂體Ⅱ相當于斜坡扇水道砂體,預測在LW3-1-B井向西北的溯源方向,砂體的規(guī)模更大。

5 結(jié)語

實踐證明,高分辨率層序地層分析在陸相小層砂體對比和有利儲層相帶預測方面已經(jīng)取得了豐碩成果。而本次研究表明:結(jié)合地震資料、測井數(shù)據(jù)及巖芯資料對海相深水扇沉積進行高分辨率層序地層學分析同樣有很好的效果。根據(jù)短期旋回層序基本類型的分布規(guī)律,可以確定不同時期主水道位置,同時在時間-地層格架內(nèi),可以對小層砂體進行對比和追蹤,描述砂體在不同時間和空間的分布特征,并預測有利儲層發(fā)育相帶。

References)

1 彭大鈞,陳長民,龐雄,等。南海珠江口盆地深水扇系統(tǒng)的發(fā)現(xiàn)[J]。石油學報,2004,25(5):17-23[Peng Dajun,Chen Changmin, Pang Xiong,et al.Discovery of deep-water fan system in South China Sea[J]。Acta Petrolei Sinica,2004,25(5):17-23]

2 彭大鈞,陳長民,龐雄,等。從淺水陸架走向深水陸坡——南海深水扇系統(tǒng)的研究[J]。石油學報,2005,23(1):1-11[Peng Dajun, Pang Xiong,Chen Changmin,et al.From shallow-water shelf to deepwater slope:The research of deep-water fan systems in South China Sea [J]。Acta Sedimentologica Sinica,2005,23(1):1-11]

3 彭大均,龐雄,陳長民,等。南海珠江深水扇系統(tǒng)的形成特征與控制因素[J]。沉積學報,2006,24(1):10-19[Peng Dajun,Pang Xiong,Chen Changmin,et al.The characteristics and controlling factors for the formation of deep-water fan system in South China Sea[J]。Acta Sedimentologica Sinica,2006,24(1):10-19]

4 彭大均,龐雄,黃先律,等。南海珠江深水扇系統(tǒng)的形成模式[J]。石油學報,2007,28(5):7-1[Peng Dajun,Pang Xiong, Huang Xianlü,et al.Depositionalmodel of Pearl River deep-water fan system in South China Sea[J]。Acta Petrolei Sinica,2007,28(5): 7-11]

5 Cross T A,Lessenger M A.Sediment volume partitioning:rationale for stratigraphicmodel evaluation and high-resolution stratigraphic correlation[M]∥Gradstein F M,Sandvik K O,Milton N J,eds.Sequence Stratigraphy Concepts and Applications.NPF Special Publication, 1998,8:171-195

6 鄭榮才,彭軍,吳朝容。陸相盆地基準面旋回的級次劃分和研究意義[J]。沉積學報,2001,19(2):249-245[Zheng Rongcai,Peng Jun,Wu Chaorong.Grade division of base-level cycles of terrigenous basin and its implications[J]。Acta Sedimentologica Sinica,2001,19 (2):249-25]

7 龐雄,陳長民,邵磊,等。白云運動:南海北部漸新統(tǒng)-中新統(tǒng)重大地質(zhì)事件及其意義[J]。地質(zhì)論評,2007,53(2):145-15[Pang Xiong,Chen Changmin,Shao Lei,et al.Baiyun movement:a great tectonic eventon the Oligocene-Miocene boundary in the northern South China Sea and its implications[J]。Geological Review,2007,53(2): 145-152]

8 柳保軍,申俊,龐雄,等。珠江口盆地白云凹陷珠海組淺海三角洲沉積特征[J]。石油學報,2007,28(2):49-61[Liu Baojun,Shen Jun,Pang Xiong,et al.Characteristics of continental delta deposits in Zhuhai Formation of Baiyun Depressionin Pearl River Mouth Basin[J]。Acta Petrolei Sinica,2007,28(2):49-56]

9 鄭榮才,尹世民,彭軍?；鶞拭嫘亟Y(jié)構(gòu)與疊加樣式的沉積動力學分析[J]。沉積學報,2000,18(3):369-375[Zheng Rongcai,Yin Shiming,Peng Jun.Sedimentary dynamic analysis of sequence structure and stacking pattern of base-level cycle[J]。Acta Sedimentologica Sinica,2000,18(3):369-375]

10 鄭榮才,彭軍,彭光明,等。高分辨率層序分析在油藏開發(fā)工程中的應用——以百色盆地侖35塊那二段油藏為例[J]。沉積學報,2003,21(4):654-662[Zheng Rongcai,Peng Jun,Peng Guangming,et al.Analysis of high-rsolution sequence stratigraphy of the second member of Nadu Formation in Lun-35 block of Baise basin and its application in developmentof oil reservoir[J]。Acta Sedimentologica Sinica,2003,21(4):654-662]

11 Vail P R.Seismic-stratigraphy interpretation using sequence stratigraphy.Pt。1:Seismic stratigraphy interpretation procedures[C]∥Bally A W.Atlas of seismic stratigraphy.AAPG Stedies in Geology,27, 1987:1-10

High-Resolution Sequence Stratigraphy Analysis and Application in Deepwater Fan of Zhujiang Formation,Liwan3-1 Area

WANG Chang-yong1YANG Bao-quan2GAO Bo-yu2HU Xiao-qing2ZHU Guo-jin2ZHENG Rong-cai
(1.Institute of Sedimentary Geology,Chengdu University of Technology,Chengdu 610059; 2.Research Institute,CNOOC,Beijing 100027)

Through the comprehensive analysis of cores,well logs and seismic interpretation data,and under the guidance of the theory and technology of the high-resolution sequence stratigraphy,the characteristics of the different grade sequence interface and marine flooding surface in the Lower Zhujiang Formation Miocene in Liwan3-1 Block are studied.In this paper,the sequence structure and distribution of the short sequence of deep water fan are discussed, and the sequence and time stratigraphic framework ofmain reservoir is established.The small sand body layers in stratigraphic framework of the Zhujiang Formation in the study area are isochronously correlated and wedged,and the reservoir developmental features of different base-level cycles are discussed.It indicates that the canyon on the continental slope is an activity location of sedimentary process.The canyon was filled by sandy silt quickly in base-level ascending phase.A majority of sedimentwas carried tomuch deeper sea area than the study area,and only little sedi-ment preserved or even caused stratigraphic lacuna in the study area.The A/S value(the change of the accommodation space to sediment supply ratios)of short sequence increasing in a vertical direction,and it indicates that deepwater fan wasmoving towards to shore land.Themoving of themain channel of deepwater fan is an important reason for the stratigraphic lacuna under the channel.The deposit in the deepwater fan channel ismost active in base-level initial ascending phase,characteristics as active progradation sequence that have downlap and onlap,and itwill be beneficial to develop reservoir sand bodies.Effective accommodation space and deepwater fan moving towards to the source area accompany with sea level rising to themaximum flooding surface.It developed retrogradation sequence in this period and argillic layers grow in number because of the channelmoving frequently.Variation of the accommodation space corresponding to base-level ascending or descending is the major factor that control short base-level cycle types and reservoir distribution in deep water fan system.Accompanying with sea level descending and effective accommodation spacemoving towards to the deep-sea plain,sediment in the deepwater fan channel in the study areawas carried to deeper sea area by gravity current,and itwould bemore beneficial to develop reservoir sand bodies to the south of Liwan 3-1 area.The results of GR inversion indicate that sandstone characteristics as strong wave trough amplitude and low wave impedance attribute.The sandbodys stretch in N-Sor NW-SE directions.The channel sandbody developed in the northland of the study area while themudstone of continental slope developed in the southland.The sand body distribution of deep water fan in different stage is described by sequence-lithofacies paleogeographic mapping and GR inversion,and the temporal and spatial distribution model of channel sand body of deep water fan have been established.The results are useful in predicting sand body of deep water fan reservoir.

deep water fan;high-resolution sequence stratigraphy;sequence structure;Zhujiang Formation;Baiyun Sag;Liwan area

王昌勇 男 1981年出生 博士后 沉積學及石油地質(zhì) E-mail:qisiwoye_2001@163。com

P539.2

A

1000-0550(2011)06-1122-08

①國家重大專項(編號:2008ZX05056-02-02)資助。

2010-10-23;收修改稿日期:2011-01-10