王亞奇，湯 軍，祝志超

(長江大學(xué) 地球科學(xué)學(xué)院，湖北 武漢 430100)

0 前 言

在油氣地質(zhì)研究過程中，有大量測繪數(shù)據(jù)、地層巖性數(shù)據(jù)、地球物探數(shù)據(jù)、油氣田生產(chǎn)數(shù)據(jù)等，然而數(shù)據(jù)量大、數(shù)據(jù)分散、缺乏有效分析方法、數(shù)據(jù)格式眾多、數(shù)據(jù)傳輸效率慢、數(shù)據(jù)重復(fù)這幾個問題位列數(shù)據(jù)使用者最頭疼問題前幾名[1]。面對復(fù)雜的數(shù)據(jù)，我們需要建立數(shù)據(jù)模型對其進(jìn)行儲存、管理、分析和決策，生成并輸出各種油氣地質(zhì)信息，從而為地質(zhì)研究提供新的評價因素，為油田企業(yè)生產(chǎn)提供重要的決策依據(jù)。

時空數(shù)據(jù)模型的研究從上個世紀(jì)60年代開始，Gail Langran[2]在1992年總結(jié)了時空立方體、快照序列、基態(tài)修正和時空復(fù)合4種時空數(shù)據(jù)模型[2]。按照時間在數(shù)據(jù)模型中所起的作用，大體可分為3個時期[3]：第1個時期側(cè)重記錄實(shí)體時態(tài)變化，稱為時態(tài)快照時期，此時期提出的時空數(shù)據(jù)模型包括時空立方體模型[4]、快照序列模型[5]，以及基態(tài)修正模型[2]、離散格網(wǎng)單元列表模型[6]、時空復(fù)合模型[7]、非第一范式時空數(shù)據(jù)模型[8]等；第2個時期側(cè)重表達(dá)實(shí)體變化的前后關(guān)系，稱為對象變化時期，此時期提出的時空數(shù)據(jù)模型包括面向?qū)ο蟮臅r空數(shù)據(jù)模型[9]、基于圖論的時空數(shù)據(jù)模型[10]、面向過程的時空數(shù)據(jù)模型[11]等；第3個時期側(cè)重描述實(shí)體變化語義，稱為事件與活動時期，此時期提出的時空數(shù)據(jù)模型包括基于事件的模型[12]、時空三域模型[13]、基于本體的時空數(shù)據(jù)模型[14]等?，F(xiàn)有的模型已經(jīng)在各行業(yè)中進(jìn)行了很多應(yīng)用，但在通用性方面還有不足，在地質(zhì)構(gòu)造方面沒有得到有效的應(yīng)用。

在油氣地質(zhì)研究與決策過程中，地質(zhì)構(gòu)造的研究具有極其重要的意義。其研究的主要內(nèi)容包括構(gòu)造的形態(tài)、產(chǎn)狀、分布和組合形式等，更重要的是我們需要研究構(gòu)造的形成條件、形成機(jī)制、形成時間、先后順序和演變規(guī)律。其中，構(gòu)造類型、空間位置、構(gòu)造變化、變化方式及時間等等都可以進(jìn)行數(shù)字化管理，而傳統(tǒng)的GIS數(shù)據(jù)模型一般只能描述現(xiàn)實(shí)世界的靜態(tài)現(xiàn)象[15]，難以滿足以上要求。因此，我們需要通過地震、測井等數(shù)據(jù)分析，設(shè)計一個地質(zhì)構(gòu)造時空數(shù)據(jù)模型。

該模型的建立，可以有助于闡明構(gòu)造在空間上的相互關(guān)系和時間上的發(fā)育順序，明晰構(gòu)造的運(yùn)動規(guī)律及其動力來源，探明其客觀規(guī)律。進(jìn)而指導(dǎo)生產(chǎn)實(shí)踐，對解決礦產(chǎn)分布、水文地質(zhì)、工程地質(zhì)、地震及環(huán)境地質(zhì)等方面問題提供必要的幫助。

1 地質(zhì)構(gòu)造時空數(shù)據(jù)模型的設(shè)計

1.1 時空變化表達(dá)

復(fù)雜的地理現(xiàn)象是由不同數(shù)量的地理對象構(gòu)成的，其表現(xiàn)為多個地理對象隨時間變化并相互作用，地理對象之間的相互作用可以用事件來傳遞。事件是指地理對象變化達(dá)到某種臨界狀態(tài)時生成的，并且傳遞給其他相關(guān)的地理對象，在某種條件下又可以驅(qū)動相關(guān)地理對象發(fā)生相應(yīng)的變化，地理對象的變化通過該對象的狀態(tài)序列來記錄。例如褶皺演化過程的地理對象是各個巖層，事件是巖層擠壓和巖層拉張，以上兩個事件是由巖層受力發(fā)生的彎曲變形屬性變化達(dá)到某種臨界狀態(tài)時引起的，使用狀態(tài)記錄巖層某一時刻的屬性變化。

根據(jù)以上分析，本文提出一個構(gòu)造時空數(shù)據(jù)模型，用于記錄與管理在復(fù)雜地理現(xiàn)象變化過程中涉及的時空數(shù)據(jù)。首先，對該模型中相關(guān)要素的概念進(jìn)行說明。

地理對象：各種構(gòu)造類型的抽象表達(dá)，由一般屬性和可變屬性組成。

事件：地理對象顯著變化的一次過程。

狀態(tài)：各種構(gòu)造類型可變屬性在某一時刻所表現(xiàn)出來的形態(tài)，用來表現(xiàn)構(gòu)造的時空變化。

圖層：具有共同結(jié)構(gòu)和功能的地質(zhì)構(gòu)造集合。

構(gòu)造時空數(shù)據(jù)模型中各個要素及其相互關(guān)系，如圖1所示。

圖1 構(gòu)造時空數(shù)據(jù)概念模型

時空過程是構(gòu)造演化的總稱，包含著有限多個地理對象和事件。地理對象表現(xiàn)為其所包含的一般屬性和可變屬性，其中我們把可變屬性記錄在狀態(tài)序列中，圖2所示。

圖2 地理對象的狀態(tài)

每個狀態(tài)記錄著該地理對象可變化部分某個時刻的快照。在表達(dá)地理對象某一時刻的整體狀態(tài)時，可以通過時態(tài)屬性查詢相應(yīng)的可變屬性和不可變屬性，將其合并到一起。

1.2 建立構(gòu)造時空數(shù)據(jù)模型

構(gòu)造時空數(shù)據(jù)模型應(yīng)能夠表達(dá)構(gòu)造類型、空間屬性、時態(tài)屬性語義。

根據(jù)上節(jié)中對地理現(xiàn)象時空變化相關(guān)要素及其互關(guān)系的分析，本節(jié)建立面向?qū)ο蟮臉?gòu)造時空數(shù)據(jù)模型，圖3給出了該模型的簡圖。

圖3 構(gòu)造時空數(shù)據(jù)模型簡圖

2 地質(zhì)構(gòu)造時空數(shù)據(jù)模型應(yīng)用實(shí)例

在構(gòu)造地質(zhì)學(xué)里，構(gòu)造可以分為褶皺、斷層、節(jié)理與劈理等。構(gòu)造時空數(shù)據(jù)模型和相關(guān)軟件中用到的所有數(shù)據(jù)都需要指定相關(guān)的數(shù)據(jù)類型。油田上，要建立一個地區(qū)的構(gòu)造時空數(shù)據(jù)模型，就必須有褶皺、斷層的相關(guān)字段項(xiàng)，如“斷層類型”、“斷層走向”、“斷層傾向”等。

同時，無論是勘探還是開發(fā)，都涉及到了大量原始數(shù)據(jù)，為方便做出智能預(yù)測和高效管理，以實(shí)現(xiàn)對最終油藏的有利預(yù)測，這就需要明確構(gòu)造狀態(tài)中各種屬性的定義。

2.1 斷層的時空數(shù)據(jù)模型

構(gòu)造帶涉及到斷層及其相關(guān)的一系列參數(shù)，首先，要確定每一屬性字段是否為可變屬性，表1給出了圖3中斷層狀態(tài)屬性的定義，其中，屬性類型值n取0表示該屬性為一般屬性，不隨時間的推移而變化；n取1表示可變屬性，該屬性可能隨時間推移而發(fā)生變化。

表1 斷層屬性數(shù)據(jù)結(jié)構(gòu)

以獅子溝構(gòu)造帶為例，平面圖如圖4所示。

圖4 柴達(dá)木盆地獅子溝構(gòu)造帶平面圖(據(jù)吳嬋等，2013改)

構(gòu)造帶的演化模式可以分為4個階段(圖5)：①路樂河組(E1+2)—下干柴溝組(E32)沉積末期，構(gòu)造帶以近南北向拉張為主，形成少量的基底正斷層；②上干柴溝組(N1)—下油砂山組(N21)，構(gòu)造帶處于坳陷期，地層沉積平穩(wěn)，沒有發(fā)生明顯的變形，斷裂發(fā)育較少，深部地層發(fā)育楔形構(gòu)造；③上油砂山組(N22)沉積初期，構(gòu)造帶受近南北向的壓扭作用；④獅子溝組(N23)沉積初期，造帶再次強(qiáng)烈擠壓，裂活動。

圖5 獅子溝構(gòu)造帶構(gòu)造演化剖面圖(據(jù)隋立偉等，2014改)

由剖面圖整理得不同時期構(gòu)造狀態(tài)原始數(shù)據(jù)如表2。

表2 獅子溝各沉積時期斷層屬性

由上文圖表可知在路樂河—下干柴溝組沉積時期，只有兩個正斷層，則輸入構(gòu)造狀態(tài)表格，例如①斷層構(gòu)造狀態(tài)表格中，輸入“編號11001；發(fā)生時間為路樂河—下干柴溝組沉積時期；一般屬性：斷層名稱；可變屬性：斷層類型、斷層走向、斷層傾向、穿越地層、斷距變化”。

其中，編號第1位為表格代碼，“1”表示構(gòu)造狀態(tài)表格，“2”表示構(gòu)造類型表格等；第2位為構(gòu)造類型編碼，“1”表示斷層，“2”表示褶皺等；后3位表示該構(gòu)造狀態(tài)編碼。

同時，注意到在不同時期，斷層發(fā)生變化，以①斷層類型和穿越地層為例，如表3。

表3 斷層歷史區(qū)

同理，可整理出各時期③斷層～⑦斷層的構(gòu)造狀態(tài)，由所有構(gòu)造狀態(tài)的集合共同生成構(gòu)造類型表格，如圖6所示。

圖6 構(gòu)造狀態(tài)集合生成構(gòu)造類型

其中，在不同時期發(fā)生構(gòu)造變化導(dǎo)致生成不同斷層，可輸入事件表格，例如輸入“編號31001；發(fā)生時間為下油砂山組沉積時期；生成構(gòu)造：③斷層”。

如上文所述，整理出斷層事件集合，同時把具有共同結(jié)構(gòu)和行為特征的構(gòu)造類型輸入同一圖層，在圖層中能夠隨時添加和移除構(gòu)造類型。

斷層、褶皺都為構(gòu)造類型，其中還可總結(jié)節(jié)理、劈理等不同的構(gòu)造類型分別添加在不同的圖層中，通過圖層對相同類別的構(gòu)造進(jìn)行統(tǒng)一樣式設(shè)置和屬性定義。同時，各種構(gòu)造類型發(fā)生緩慢或顯著的變化，在構(gòu)造地質(zhì)學(xué)中，一般構(gòu)造變化均以百萬年計算，也偶有劇烈構(gòu)造變化(如地震、火山爆發(fā)等)，形成新的構(gòu)造類型，記錄在事件中。事件是模型的重要組成部分，記錄了生成此事件的構(gòu)造類型及生成時間。

圖層中包含各種構(gòu)造類型的集合，事件中包含構(gòu)造發(fā)生的變化及新生成的構(gòu)造集合，由其組成圖層集合和事件集合，共同包含在構(gòu)造演化過程中(見圖3)。

2.2 模型的時空數(shù)據(jù)查詢過程

用戶在某一時刻對空間數(shù)據(jù)進(jìn)行查詢操作時，例如查詢①斷層在上油砂山組沉積時期的穿越地層，就要在數(shù)據(jù)庫中選取該構(gòu)造狀態(tài)，繼而選取構(gòu)造狀態(tài)發(fā)生的時間：上油砂山組沉積時期，再從該構(gòu)造狀態(tài)中選取屬性穿越地層，就可知該對象的屬性信息和事件信息。這樣，我們可以通過建立的構(gòu)造時空數(shù)據(jù)模型對構(gòu)造演化進(jìn)行數(shù)據(jù)存儲和規(guī)律研究。例如調(diào)取過去某一時間段內(nèi)各構(gòu)造類型、狀態(tài)屬性、事件和當(dāng)時構(gòu)造地區(qū)的平面圖(如圖4)，就可以嘗試推演該地區(qū)更早時間段的構(gòu)造平面圖和該構(gòu)造的生成機(jī)理，進(jìn)而可以對該構(gòu)造帶進(jìn)行三維構(gòu)造地質(zhì)建模的探索。

3 結(jié) 語

1)分析了復(fù)雜地理現(xiàn)象時空變化的特點(diǎn)，提出了一種構(gòu)造時空數(shù)據(jù)模型，將構(gòu)造演化過程、不同構(gòu)造類型、構(gòu)造變化、構(gòu)造狀態(tài)等要素有機(jī)地結(jié)合在一起。并且通過對柴達(dá)木盆地獅子溝構(gòu)造帶數(shù)據(jù)的輸入，驗(yàn)證了模型的可行性和易用性。該模型符合相對通用的地質(zhì)構(gòu)造時空數(shù)據(jù)模型的特點(diǎn)。而采用面向?qū)ο蟮脑O(shè)計思想，也使得模型具有良好的可實(shí)現(xiàn)性和擴(kuò)展性，為油田在各個領(lǐng)域的擴(kuò)展應(yīng)用研究提供基本框架。

2)建立構(gòu)造時空數(shù)據(jù)模型，通過對構(gòu)造資料的處理及分析，不僅可以解決歷史數(shù)據(jù)和現(xiàn)時數(shù)據(jù)的整合問題，且實(shí)現(xiàn)了基于時空變化語義的查詢和分析?？梢杂糜谟筒孛枋鲋?，完善油藏的精細(xì)描述；將構(gòu)造的時空數(shù)據(jù)合力組織、存儲、管理起來，有利于完善數(shù)字油田的建設(shè)。

總體來說，構(gòu)造演化與其他地理對象不同，構(gòu)造的變化通常以百萬年記，因此構(gòu)造的時空變化是一個非常漫長的過程。如何依據(jù)構(gòu)造時空數(shù)據(jù)模型，在事件中加入變量對構(gòu)造演化進(jìn)行準(zhǔn)確的動態(tài)模擬是接下來需要研究和探討的內(nèi)容。

[1] 高倩.多源數(shù)字油田數(shù)據(jù)與油田數(shù)據(jù)融合及方法研究[D].西安：長安大學(xué),2015.

[2] Langran G. Time in geographic information systems[M]. London:Taylor & Francis Ltd, 1992.

[3] Michael Worboys. Event-oriented approaches to geographic phenomena[J]. International Journal of Geographical Information Science, 2005, 19(1):1-28.

[4] H?gerstraand T. What about people in reglonal science?[J]. Papers of the Regional Science Association, 2010, 24(1):6-21.

[5] Armstrong M P. Temporality in spatial database[C]// Proceedings of GIS/LIS. San Antonio:[s.n.], 1988:886-899.

[6] Gail Langran. A review of temporal database research and its use in GIS applications[J]. International Journal of Geographical Information Science, 1989, 3(3):215-232.

[7] Langran G, Chrisman N. A framework for temporal geographic information[J]. Cartographica the International Journal for Geographic Information & Geovisualization, 1988, 25(3):1-14.

[8] 陳軍, 陳尚超, 唐治鋒. 用非第一范式關(guān)系表達(dá)GIS時態(tài)屬性數(shù)據(jù)[J]. 武漢大學(xué)學(xué)報(信息科學(xué)版), 1995, 20(1):12-17.

[9] Worboys M F. A model for spatio-temporal information[C]// Proceedings of the 5th Interantional Symposium on Spatial Data Handling: 2. Columbia:University of South Carolina, 1992:602-611.

[10] 尹章才, 李霖, 艾自興. 基于圖論的時空數(shù)據(jù)模型研究[J]. 測繪學(xué)報, 2003, 32(2):168-172.

[11] Pang Y C. Development of process-based model for dynamic interaction process in spatio-temporal GIS[J]. Geoinformatica, 1999, 6(4):323-344.

[12] Peuquet D J, Duan N. An event-based spatiotemporal data model (ESTDM) for temporal analysis of geographical data[J]. International Journal of Geographical Information, 1995, 9(1):7-24.

[13] Yuan M. Representation of wildfire in geographic information systems[D]. Buffalo: State University of New York, 1994.

[14] Galton A. Desiderata for a spatio-temporal geo-ontology[C]// Spatial Information Theory: Foundations of Geographic Information Science. Berlin/Heidelberg: Springer, 2001:1-12.

[15] Peuquet D J. A conceptual framework and comparison of spatial data models[J]. Cartographica the International Journal for Geographic Information & Geovisualization, 1984, 21(4):66-113.