基于GOCAD的成都地鐵聲波參數(shù)可視化應(yīng)用研究

楊期祥，王森，牟迪，林剛，巫錫勇，廖昕

(1.西南交通大學(xué) 地質(zhì)工程系，四川 成都 611756；2.中鐵二院工程集團(tuán)有限責(zé)任公司 地下鐵道設(shè)計(jì)研究院，四川 成都 610031)

?

基于GOCAD的成都地鐵聲波參數(shù)可視化應(yīng)用研究

楊期祥1，王森1，牟迪2，林剛2，巫錫勇1，廖昕1

(1.西南交通大學(xué) 地質(zhì)工程系，四川 成都 611756；2.中鐵二院工程集團(tuán)有限責(zé)任公司 地下鐵道設(shè)計(jì)研究院，四川 成都 610031)

摘要:利用GOCAD軟件的參數(shù)建模功能，將成都地鐵某段線路沿線聲波測(cè)試法測(cè)得的各層巖土體的縱、橫波傳播速度作為區(qū)域化隨機(jī)變量進(jìn)行屬性地質(zhì)模型的構(gòu)建，沿地鐵走向建立的三維地質(zhì)模型不僅可以體現(xiàn)沿線地層空間分布規(guī)律特征，而且可以直觀展現(xiàn)出地質(zhì)體的動(dòng)彈性力學(xué)特性，實(shí)現(xiàn)了地質(zhì)體空間結(jié)構(gòu)展布與力學(xué)特征的結(jié)合和可視化，并且通過(guò) GOCAD軟件提供的地質(zhì)統(tǒng)計(jì)分析方法來(lái)預(yù)測(cè)模型中任意點(diǎn)的參數(shù)值，更有助于觀察聲波參數(shù)在研究區(qū)域內(nèi)的變化趨勢(shì)，可簡(jiǎn)捷高效對(duì)地鐵沿線建筑土類型及場(chǎng)地級(jí)別進(jìn)行劃分，從而方便工程設(shè)計(jì)人員更加便捷詳細(xì)地獲得直觀、準(zhǔn)確的工程地質(zhì)信息。

關(guān)鍵詞：GOCAD；參數(shù)建模；三維地質(zhì)模型；空間分布；地質(zhì)統(tǒng)計(jì)

隨著城市地鐵建設(shè)快速發(fā)展，對(duì)于工程地質(zhì)勘察工作也提出了更高的要求。如何結(jié)合地質(zhì)勘察為地鐵的設(shè)計(jì)及施工提供詳實(shí)的地質(zhì)資料，如何準(zhǔn)確有效地評(píng)價(jià)建設(shè)場(chǎng)地的工程地質(zhì)條件，對(duì)于城市地鐵建設(shè)意義重大。而對(duì)于傳統(tǒng)的工程地質(zhì)分析方法，一方面，工程地質(zhì)資料的分析和解釋多以文字、圖表等二維、靜態(tài)的形式展示，描述空間地質(zhì)構(gòu)造起伏變化的直觀性差，不能充分揭示地層的空間分布規(guī)律；另一方面，傳統(tǒng)勘察方法在實(shí)際工作中操作難度較高。尤其在城市地鐵建設(shè)中，限于沿線高樓林立、地表建筑物及地下管線密集，實(shí)際工作中經(jīng)常遇到無(wú)勘察工作面，無(wú)法按設(shè)計(jì)進(jìn)行布點(diǎn)的情況出現(xiàn)。而隨著現(xiàn)代計(jì)算機(jī)技術(shù)的發(fā)展及三維地質(zhì)理論的逐漸成熟，三維地質(zhì)建模軟件GOCAD開始應(yīng)用于輔助傳統(tǒng)的地質(zhì)勘查工作[1-3]。利用GOCAD可以將野外勘察工作中得到的地質(zhì)資料經(jīng)過(guò)整理歸類后以數(shù)據(jù)庫(kù)的形式保存并以立體的、與真實(shí)地質(zhì)體具有同樣屬性的三維地質(zhì)模型呈現(xiàn)；三維地質(zhì)模型的建立還可以進(jìn)一步指導(dǎo)野外勘查工作，輔助勘察人員進(jìn)行合理的鉆孔布置，使勘察工作在實(shí)際應(yīng)用中操作性更強(qiáng)，甚至可以利用鉆孔資料插值擬合出周圍地層數(shù)據(jù)，直接在地質(zhì)模型中獲得部分控制孔的勘察信息，而無(wú)需實(shí)際鉆孔；更重要的是利用GOCAD參數(shù)建模可以實(shí)現(xiàn)地質(zhì)體物理力學(xué)參數(shù)的可視化、直觀化[4-6]。建立的屬性參數(shù)模型不僅能夠用于地質(zhì)信息的查詢，而且可以通過(guò)多種地質(zhì)統(tǒng)計(jì)分析方法對(duì)地質(zhì)模型中任一點(diǎn)的物理力學(xué)參數(shù)值進(jìn)行預(yù)測(cè)以及觀察該參數(shù)在研究區(qū)域的變化趨勢(shì)，從而為設(shè)計(jì)施工提供全面詳實(shí)的地質(zhì)數(shù)據(jù)，有利于設(shè)計(jì)人員對(duì)該區(qū)域工程地質(zhì)條件進(jìn)行綜合分析評(píng)價(jià)。目前，國(guó)內(nèi)將GOCAD應(yīng)用于工程中的相關(guān)研究相對(duì)較少，大都用于建立簡(jiǎn)單的小范圍三維立體模型[7-9]，基于GOCAD的大區(qū)域參數(shù)建模還鮮有研究，張聚興等[10-11]將GOCAD 技術(shù)應(yīng)用于城市環(huán)境地質(zhì)評(píng)價(jià)取得了較好的效果，但對(duì)于參數(shù)模型的方法及原理并未作出明確解釋。董慧超等[12-13]將GOCAD三維建模應(yīng)用于城市地鐵建設(shè)工程，但所建立的模型較粗糙，未能結(jié)合地鐵走向?qū)δＰ瓦M(jìn)行優(yōu)化與改進(jìn)。基于以上研究，結(jié)合地鐵工程走向，選取有效參數(shù)將其應(yīng)用于GOCAD進(jìn)行空間參數(shù)立體化展示，是本文擬解決的關(guān)鍵問(wèn)題。 在城市地鐵施工中，尤其是對(duì)于成都這種典型的砂卵石地層，與巖土體的彈性波速有關(guān)的動(dòng)彈性力學(xué)參數(shù)是設(shè)計(jì)、施工中一項(xiàng)十分重要的影響指標(biāo)[14]。該指標(biāo)不僅可以反映地層地質(zhì)特征以及物理力學(xué)性質(zhì)，實(shí)現(xiàn)對(duì)圍巖體的工程地質(zhì)分類及穩(wěn)定性綜合評(píng)價(jià)，測(cè)得動(dòng)彈性指標(biāo)后結(jié)合室內(nèi)物理力學(xué)實(shí)驗(yàn)進(jìn)行分析比較，還能保證測(cè)量方法的多樣性以及測(cè)量數(shù)據(jù)的可靠性、正確性。目前常用的測(cè)量方法是原位聲波測(cè)試法，即彈性波法，原位聲波測(cè)試是利用脈沖震源激發(fā)彈性波,彈性波在巖體中傳播時(shí)，其聲學(xué)參數(shù)如速度、振幅、頻譜等將受到介質(zhì)體的巖性、結(jié)構(gòu)、風(fēng)化程度等影響而發(fā)生變化[15]。這些變化將直接反映出巖體的地質(zhì)構(gòu)造和物理力學(xué)性質(zhì)。相較于靜力法，原位動(dòng)態(tài)法的優(yōu)點(diǎn)是原位動(dòng)態(tài)、測(cè)試面廣、設(shè)備輕便、周期短、經(jīng)濟(jì)且易于掌握。目前，工程中對(duì)于將原位聲波測(cè)試得出的動(dòng)彈性力學(xué)參數(shù)大都只是作為靜態(tài)的試驗(yàn)數(shù)據(jù)進(jìn)行整理和分析，并未考慮將該參數(shù)與三維地質(zhì)模型結(jié)合分析。通過(guò)GOCAD參數(shù)建模不僅可以實(shí)現(xiàn)動(dòng)彈性力學(xué)參數(shù)在空間上的三維分布，及時(shí)提供信息指導(dǎo)現(xiàn)場(chǎng)施工，還可以結(jié)合地表實(shí)際情況進(jìn)行建筑場(chǎng)地的穩(wěn)定性評(píng)價(jià)，甚至能夠結(jié)合地區(qū)抗震等級(jí)，對(duì)建筑場(chǎng)地進(jìn)行抗震設(shè)計(jì)。本文收集整理成都地鐵某段線路的勘察鉆孔資料和室內(nèi)土工試驗(yàn)數(shù)據(jù)，利用GOCAD軟件建立三維地質(zhì)模型，在此基礎(chǔ)上，將地鐵沿線聲波測(cè)試法測(cè)得的各層巖土體的縱、橫波傳播速度作為區(qū)域化變量賦予地質(zhì)體模型中并建立參數(shù)模型，實(shí)現(xiàn)了地質(zhì)體的空間結(jié)構(gòu)展布與物理力學(xué)性質(zhì)有效結(jié)合，可為類似地鐵工程地質(zhì)建模提供借鑒。

1三維地質(zhì)模型的構(gòu)建

為建立參數(shù)模型，首先應(yīng)當(dāng)實(shí)現(xiàn)三維地質(zhì)體模型的構(gòu)建，三維地質(zhì)模型可以分為構(gòu)造模型(structural modeling)和三維儲(chǔ)層?xùn)鸥窠Y(jié)構(gòu)模型(3D Reservoir Grid Construction)，構(gòu)造模型主要用來(lái)模擬地層面、斷層面的形態(tài)和相互位置關(guān)系，可以用于地質(zhì)構(gòu)造分析以及地震勘探中地震反演的重要手段；三維儲(chǔ)層?xùn)鸥窠Y(jié)構(gòu)建模則可以將地質(zhì)體的巖層特性和資源分布等特性參數(shù)導(dǎo)入，生成地質(zhì)體參數(shù)模型。因此，本次建模采用的是三維儲(chǔ)層?xùn)鸥窠Y(jié)構(gòu)模型。通常三維地質(zhì)建模主要有以下幾個(gè)步驟：點(diǎn)模型的構(gòu)建，面模型的構(gòu)建和體模型的構(gòu)建。而點(diǎn)模型及面模型的構(gòu)建通常包括數(shù)據(jù)分析整理及離散數(shù)據(jù)插值擬合2步。

根據(jù)勘察鉆探揭露分析，依據(jù)該場(chǎng)地巖土層特性，可將地鐵沿線土層主要分為4層，如表1所示。

表1　地鐵沿線地層巖性

通過(guò)收集和整理沿地鐵走向總共8個(gè)站的112個(gè)鉆孔資料，并按照如下步驟依次進(jìn)行GOCAD三維地質(zhì)建模。

1)建立點(diǎn)模型。首先根據(jù)鉆孔分層資料提取各個(gè)鉆孔分層點(diǎn)的(X，Y，Z)數(shù)據(jù)導(dǎo)入sufer軟件進(jìn)行 Kriging 插值，以點(diǎn)的形式(pointset)導(dǎo)入 GOCAD 中，這時(shí)GOCAD生成的是一些三維的坐標(biāo)點(diǎn)。為直觀清晰地反映導(dǎo)入的坐標(biāo)點(diǎn)與實(shí)際線路的關(guān)系，根據(jù)點(diǎn)的分布繪制出一條擬合曲線。如圖1～2所示分別為成都地鐵某段線路走向示意圖和導(dǎo)入的原鉆孔布置點(diǎn)擬合曲線 ，對(duì)比可知，地鐵走向與原鉆孔控制點(diǎn)擬合曲線十分吻合，證明收集和整理的鉆孔點(diǎn)數(shù)據(jù)是有效和準(zhǔn)確的。

圖1　地鐵線路走向示意圖(上為平面圖，下為縱斷面圖)Fig.1 Stretch figure of the subway line(the above is plane figure, the below is longitudinal profile)

圖2　鉆孔控制點(diǎn)擬合曲線Fig.2 Fitting curveof borehole control points

2)建立面模型。以各個(gè)鉆孔分層點(diǎn)作為控制點(diǎn)，形成 Delaunay 三角網(wǎng)格化的曲面，并利用 GOCAD 軟件里的 DSI 插值算法優(yōu)化各層初始曲面；在地層尖滅或者缺失區(qū)域獲得其交線，然后利用交線約束修改原尖滅、缺失地層曲面從而最終建立各地層的層面，由于第5層卵石土厚度較大，鉆孔最底部揭露地層即為該層，因此本次建模將該層添加為底面，且所有的Z坐標(biāo)都統(tǒng)一定為鉆孔揭露的該層最大深度。如圖3～4所示。

為了直觀地呈現(xiàn)沿地鐵走向的三維地質(zhì)模型，需要對(duì)已建立的各層初始曲面沿地鐵走向線進(jìn)行切割，即先將地鐵走向線的Z方向進(jìn)行延伸，生成側(cè)面作為切割面，然后沿交線對(duì)各層初始曲面進(jìn)行切割，得到切割后的各層曲面，如圖5～6所示。

圖3　面層插值數(shù)據(jù)Fig.3 Surface interpolation data

圖4　各層初始曲面Fig.4 Initial surface of each layer

圖5　各層初始曲面切割圖Fig.5 Cutting figure of initial surface

圖6　切割后的各層曲面Fig.6 Each surface after cutting

3)建立體模型。GOCAD中提供了2種體模型的建立方法，一種是實(shí)體模型對(duì)象(Solid)，另一種是網(wǎng)格模型對(duì)象(SGrid)。Solid模型相對(duì)簡(jiǎn)單，主要是對(duì)地質(zhì)體幾何特征和空間結(jié)構(gòu)的模擬，不能賦予模型任何物理力學(xué)參數(shù)信息，因此在實(shí)際工程中應(yīng)用較少。而SGrid模型卻能彌補(bǔ)上述缺點(diǎn)，不僅可以建立三維模型，還能導(dǎo)入地質(zhì)體物理力學(xué)參數(shù)信息，并且SGrid模型的柵格節(jié)點(diǎn)信息和拓?fù)潢P(guān)系還可以導(dǎo)入到FLAC以及ANSYS等數(shù)值軟件中進(jìn)行后處理分析，因?yàn)檫@些優(yōu)勢(shì)SGrid模型應(yīng)用相對(duì)廣泛。因此，本次參數(shù)建模采用了SGrid模型建立體模型。SGrid模型的建立步驟可以簡(jiǎn)單分為以下幾步：1)選擇頂?shù)酌妫謩e設(shè)置頂、底面；2)設(shè)置頂?shù)酌嬷g的連接；3)設(shè)置平面網(wǎng)格方向；4)設(shè)置縱向上網(wǎng)格方向；5)創(chuàng)建中間層單元；6)定義縱向上網(wǎng)格數(shù)；7)定義平面上網(wǎng)格數(shù)；8)創(chuàng)建三維地質(zhì)網(wǎng)格模型，參照上述步驟建立的三維地質(zhì)網(wǎng)格體模型如圖7所示。

圖7　由SGrid法做出的三維地質(zhì)模型Fig.7 Three-dimensional geological modeling generated by the method of SGrid

2三維參數(shù)模型的構(gòu)建

2.1區(qū)域化隨機(jī)變量的確定

區(qū)域化隨機(jī)變量也稱為區(qū)域化變量(Regionalized Variable)，當(dāng)一個(gè)變量呈現(xiàn)空間分布時(shí)，就稱之為區(qū)域化變量。它是一個(gè)以空間點(diǎn)的三維坐標(biāo)(x，y，z)為自變量的隨機(jī)場(chǎng)，它的顯著特征是變量的隨機(jī)性以及結(jié)構(gòu)性。假設(shè)空間中任一區(qū)域化隨機(jī)變量F(u)為空間點(diǎn)u(x，y，z)的函數(shù)，其中(x，y，z)為點(diǎn)u的三維坐標(biāo)，則任一點(diǎn)m處的區(qū)域化隨機(jī)變量即可表示為F(m)，這就是區(qū)域化變量的隨機(jī)性體現(xiàn)；而空間中任2點(diǎn)a和b的區(qū)域化隨機(jī)變量F(a)和F(b)通常具有某種相關(guān)性，即空間中同種屬性的區(qū)域化隨機(jī)變量具有一定的聯(lián)系，這就是區(qū)域化隨機(jī)變量的結(jié)構(gòu)性體現(xiàn)，例如同種土層的同種物理力學(xué)特征在空間上呈相似性分布。基于GOCAD的參數(shù)建模所采用的參數(shù)必須是區(qū)域化隨機(jī)變量，本次參數(shù)建模采用的參數(shù)是通過(guò)原位聲波測(cè)試法得出的動(dòng)彈性力學(xué)參數(shù)如橫波傳播速度，縱波傳播速度，以及動(dòng)剪切波速等。這些參數(shù)既具有隨機(jī)性又具有結(jié)構(gòu)性，所以可以被看作是區(qū)域化隨機(jī)變量。根據(jù)前人的研究分析可知[14]，橫波傳播速度只與介質(zhì)密度和剪切模量有關(guān)，反映了巖體抗剪切的能力；縱波傳播速度不僅與壓縮模量和介質(zhì)密度有關(guān)，而且也與剪切模量有關(guān)，反映了巖體壓縮與抗拉伸的能力。

2.2建立參數(shù)模型

通過(guò)對(duì)地鐵沿線每個(gè)車站區(qū)間選取2～3個(gè)代表性鉆孔中不同深度的土層進(jìn)行原位聲波測(cè)試，確定每個(gè)土層的動(dòng)彈性力學(xué)參數(shù)，并把直接量測(cè)得出的橫波傳播速度、縱波傳播速度值經(jīng)過(guò)整理后作為取樣點(diǎn)的一個(gè)屬性值進(jìn)行導(dǎo)入，即利用軟件的Reservoir Properties 功能選擇已建立好的SGrid模型，建立1個(gè)或多個(gè)屬性參數(shù)值將其添加到該網(wǎng)格模型中，根據(jù)該參數(shù)點(diǎn)類型選擇離散性或連續(xù)性，利用軟件中的DSI插值算法進(jìn)行屬性多維網(wǎng)格插值，根據(jù)實(shí)際需要選擇一種賦值到網(wǎng)格體的方法(最近原則，算術(shù)平均，幾何平均，諧和平均，反距離平均，加權(quán)平均)，然后利用 GOCAD 中的地質(zhì)統(tǒng)計(jì)Kriging插值法(簡(jiǎn)單克里金，普通克里金，外部漂移克里金，貝葉斯克里金，協(xié)克里金，趨勢(shì)克里金) 進(jìn)行地質(zhì)統(tǒng)計(jì)模擬估算參數(shù)的分布，使這些參數(shù)分布以三維體表現(xiàn)出來(lái)形成參數(shù)模型，可以形象直觀地得到某參數(shù)在整個(gè)地質(zhì)體中的變化情況。下面以橫波傳播速度，縱波傳播速度建立的參數(shù)模型為例，見圖8～9。

圖8　橫波速度參數(shù)模型Fig.8 Parameter model of shear velocity

圖9　縱波速度參數(shù)模型Fig.9 Parameter model of p-wave velocity

由圖8～9對(duì)比可看出，在空間分布上，橫波速度與縱波速度大小成正相關(guān)，即橫波速度大的地段縱波速度相應(yīng)也較大。由色標(biāo)示值還可以看出二者均隨地層深度增大，顏色趨于速度大的色區(qū)，即波速逐漸增大，地層愈密實(shí)愈堅(jiān)硬。而同一地層的波速基本上保持同一色區(qū)，即保持在一個(gè)相同的波速范圍之內(nèi)。為驗(yàn)證所建模型與實(shí)際情況是否相符，以橫波速度參數(shù)模型為例，統(tǒng)計(jì)了現(xiàn)場(chǎng)所測(cè)得的各試驗(yàn)點(diǎn)的橫波傳播速度值，與建立的橫波速度模型進(jìn)行對(duì)比分析，并依據(jù)《建筑抗震設(shè)計(jì)規(guī)范》[16](GB50011—2010)中第4.1.3條規(guī)定，將場(chǎng)地土類型劃分見表2。

表2　橫波傳播速度統(tǒng)計(jì)對(duì)比及場(chǎng)地土類型劃分

由表2可知，參數(shù)模型與實(shí)測(cè)數(shù)據(jù)吻合度較高，說(shuō)明所建模型能較好地預(yù)測(cè)橫波傳播速度在地鐵沿線上的空間分布情況，并可以根據(jù)模型色標(biāo)識(shí)別場(chǎng)地土類型，大大提高了勘察效率。同理，還可以根據(jù)各層土的橫波速度得出等效剪切波速，參照規(guī)范判別標(biāo)準(zhǔn)對(duì)沿線地鐵區(qū)間進(jìn)行建筑場(chǎng)地類別劃分。

另外，參數(shù)模型建立之后，GOCAD還具有強(qiáng)大的屬性建模后期處理功能，包括地質(zhì)統(tǒng)計(jì)(計(jì)算平均值、中值，最小值、最大值、比例等)、概率計(jì)算(計(jì)算屬性不同范圍內(nèi)的實(shí)現(xiàn)概率)、統(tǒng)計(jì)計(jì)算三維屬性，在二維圖上顯示等功能。通過(guò)對(duì)建立的聲波波速參數(shù)模型進(jìn)行分析，不僅可直接用于參數(shù)信息統(tǒng)計(jì)與查詢，還能結(jié)合室內(nèi)試驗(yàn)進(jìn)行建筑場(chǎng)地的穩(wěn)定性分級(jí)，也可間接將建立的SGrid模型導(dǎo)入FLAC(3D)和Midas/GTS等軟件中進(jìn)行數(shù)值計(jì)算，還可以方便地導(dǎo)入到MPGIS等地理信息系統(tǒng)軟件中，建立區(qū)域化的數(shù)據(jù)庫(kù)，對(duì)于工程地質(zhì)工作而言，具有十分重要的應(yīng)用價(jià)值。

3結(jié)論

1)利用有限數(shù)量的鉆孔控制點(diǎn)通過(guò)Kriging插值生成擴(kuò)大范圍的曲面后，針對(duì)地層尖滅或者缺失區(qū)域，可先在交錯(cuò)面上生成交線，然后利用交線約束修改原尖滅、缺失地層曲面，從而最終建立各地層的層面。

2)對(duì)于已建立的地層層面，若需要按照實(shí)際情況進(jìn)行切割，可先生成一條沿切割方向的曲線，并通過(guò)Z向延伸得到側(cè)面，將該側(cè)面作為切割面使之與初始層面相交，然后沿交線對(duì)各層初始曲面進(jìn)行切割，得到切割后的各層曲面。

3)對(duì)于符合區(qū)域化隨機(jī)變量的參數(shù)可以導(dǎo)入GOCAD進(jìn)行參數(shù)建模，本文選取了縱、橫波傳播速度導(dǎo)入GOCAD建立了三維地質(zhì)參數(shù)模型。通過(guò)模型參數(shù)與實(shí)測(cè)數(shù)據(jù)的對(duì)比分析，證明該參數(shù)模型能較好地反映縱、橫波速度在空間上的分布情況，表達(dá)效果較理想，能有效預(yù)測(cè)巖土體性質(zhì)分布趨勢(shì)，進(jìn)行場(chǎng)地土類別劃分及場(chǎng)地等級(jí)評(píng)價(jià)等工作。

4)在構(gòu)建好波速參數(shù)模型的基礎(chǔ)上，若能與室內(nèi)實(shí)測(cè)的各類型巖性波速試驗(yàn)進(jìn)行比較，則可以使模型參數(shù)實(shí)現(xiàn)量化與可視化，對(duì)于今后成都地鐵數(shù)據(jù)庫(kù)的建立與完善將發(fā)揮重要的輔助作用。

5)本次所建三維地質(zhì)模型，限于鉆孔數(shù)量有限和地質(zhì)分層較粗略，未能呈現(xiàn)出砂卵石層中透鏡體性質(zhì)以及空間分布情況，有待進(jìn)一步細(xì)化與加強(qiáng)；而在結(jié)合地質(zhì)構(gòu)造進(jìn)行地質(zhì)解釋等方面也有待進(jìn)一步改進(jìn)。

參考文獻(xiàn)：

[1] Marcus A. Development of a 3D GIS basedon the 3D modeler GOCAD[C]//IAMG Meeting,2001:256-268.

[2] Bell F G. Geological-evolutionary model of a gravity-induced slope deformation in the carbonate central apennines (Italy) [J]. Quarterly Journal of Engineering Geology & Hydrogeology, 2004,37(1): 31-47.

[3] 張洋洋, 周萬(wàn)蓬, 吳志春, 等. 三維地質(zhì)建模技術(shù)發(fā)展現(xiàn)狀及建模實(shí)例[J]. 東華理工大學(xué)學(xué)報(bào)(社會(huì)科學(xué)版), 2013, 32(3): 403-409.

ZHANG Yangyang, ZHOU Wanpeng, WU Zhichun,et al. The development status of 3D geological modeing technology and modeling instances[J].Journal of East China Institute of Technology(Social Science) , 2011, 28(4): 69-73.

[4] 黃靜莉, 王清. 基于 GOCAD 的三維工程地質(zhì)地層建模研究[J]. 吉林地質(zhì), 2013, 32(3): 111-114.

HUANG Jingli, WANG Qing. Research of three-dimensional engineering geology trata modeling based on GOCAD[J]. Jilin Geology, 2013, 32(3): 111-114.

[5] 邵國(guó)波, 郭艷, 蔡冰. GOCAD 軟件在三維地質(zhì)建模中的應(yīng)用[J]. 山東交通科技, 2012 (5): 25-27.

SHAO Guobo, GUO Yan, CAI Bing.GOCAD software application in three-dimensional geologicalmodeling method[J]. Shandong Transportation Technology, 2012 (5): 25-27.

[6] 向中林, 王妍, 王潤(rùn)懷, 等. 基于鉆孔數(shù)據(jù)的礦山三維地質(zhì)建模及可視化過(guò)程研究[J]. 地質(zhì)與勘探, 2009 (1): 75-81.

XIANG Zonglin, WANG Yan, WANG Runhuai, et al. 3D geological modeling and visualization process of mines based on borehole data[J]. Geology and Exploration, 2009 (1): 75-81.

[7] 張燕飛, 朱杰勇, 張威. 基于 GOCAD 的三維地質(zhì)模型構(gòu)建[J]. 河北工程大學(xué)學(xué)報(bào)(自然科學(xué)版), 2011, 28(4): 69-73.

ZHANG Yanfei, ZHU Jieyong, ZHANG Wei.3D geological modeling method based on GOCAD[J].Journal of Hebei University of Engineering(Nature Science Edition), 2011, 28(4): 69-73.

[8] 童亮, 胡卸文, 汪雪瑞. GOCAD 在某橋基邊坡三維地質(zhì)建模中的應(yīng)用研究[J]. 地質(zhì)災(zāi)害與環(huán)境保護(hù), 2009, 20(1): 60-65.

TONG Liang, HU Xiewen, WANG Xuerui. Study on application of three-dimensional modeling to a bridge slope by using the software named gocad[J]. Journal of Geological Hazards and Environment Preservation, 2009, 20(1): 60-65.

[9] 賈新會(huì), 王小兵, 曹文廣. 基于 GOCAD 平臺(tái)的地質(zhì)三維建模技術(shù)在水電工程中的應(yīng)用[J]. 西北水電, 2012 (4): 18-22.

JIA Xinhui, WANG Xiaobing, CAO Wenguang. Application of technology on 3D geological modeling based on GOCAD in hydropower engineering[J]. Northwest Water Power, 2012 (4): 18-22.

[10] 張聚興.三維地質(zhì)建模在工程地質(zhì)環(huán)境質(zhì)量評(píng)價(jià)中的應(yīng)用研究[D]. 北京: 中國(guó)地質(zhì)大學(xué), 2006.

ZANG Juxing. The study on application of three-dimensional modeling in evaluation of engineering geological environmental quality[D]. Beijing: China University of Geosciences, 2006.

[11] 宋大各, 胡卸文. GOCAD 技術(shù)在城市環(huán)境地質(zhì)評(píng)價(jià)中的應(yīng)用[J]. 地質(zhì)災(zāi)害與環(huán)境保護(hù), 2009, 20(1): 55-59.

SONG Dage,HU Xiewen. The application of gocad technology in urban environmental geology evaluation[J]. Journal of Geological Hazards and Environment Preservation, 2009, 20(1): 55-59.

[12] 董慧超. 北京地鐵九號(hào)線三維地質(zhì)建模研究[D]. 北京: 中國(guó)地質(zhì)大學(xué), 2011.

DONG Huichao. Geological modeling and analyzing of Beijing subway line 9[D].Beijing: China University of Geosciences, 2011.

[13] 王明. GOCAD 在鄭州地鐵 1 號(hào)線三維地質(zhì)建模中的應(yīng)用研究[D]. 成都：西南交通大學(xué), 2013.

WANG Ming. The study on application of three-dimensionalmodeling in Zhengzhou Subway Line 1 by usingthe software named GOCAD[D]. Chengdu: Southwest Jiaotong University, 2013.

[14] 劉晶海. 聲波測(cè)試技術(shù)在工程地質(zhì)中的應(yīng)用[J]. 吉林地質(zhì), 1997, 16(2): 71-77.

LIU Jinghai. An application of sound wave measuring technique in engineering geology[J]. Jilin Geology, 1997, 16(2): 71-77.

[15] 熊健, 胡樹林, 左勖, 等. 基于鉆孔電視和聲波測(cè)試的某地下水封洞庫(kù)巖體質(zhì)量初步評(píng)價(jià)[J]. 工程勘察, 2014, 42(8): 24-29.

XIONG Jian, HU Shuling, ZUO Xu,et al. Preliminary evaluation for rock mass quality of certain underground water enclosed cavern based on borehole television and sonic testing[J].Geological Investigation & Surveying, 2014, 42(8): 24-29.

[16] 馬中高, 解吉高. 巖石的縱, 橫波速度與密度的規(guī)律研究[J]. 地球物理學(xué)進(jìn)展, 2006, 20(4): 905-910.

MA Zhonggao, XIE Jiegao. Relationship among compressional wave,shear wave velocities and density of rocks[J]. Progress in Geophysics, 2006, 20(4): 905-910.

[17] 2010 G B, 建筑抗震設(shè)計(jì)規(guī)范[S].

2010 G B, Code for seismic design of buildings[S].

[18] 錢洪,唐榮昌. 成都平原的形成與演化[J].四川地震,1997(3):1-7.

QIAN Hong, TANG Rongchang. On the formation and evolution of the Chengdu plain[J]. Sichuan Earthquake,1997(3):1-7.

(編輯陽(yáng)麗霞)

Study on GOCAD-based visualization of acoustic parameters applied to the subway project in Chengdu

YANG Qixiang1, WANG Sen1, MOU Di2, LIN Gang2, WU Xiyong1, LIAO Xin1

(1.Department of Geological Engineering, Southwest Jiaotong University, Chengdu 611756,China;2.Metro Design and Research Institute, China Railway Eryuan Engineering Group Co. Ltd, Chengdu 610031, China)

Abstract:Based on the parametric modeling capabilities of GOCAD, this paper gets the data of the shear wave velocity and the longitudinal wave velocity which were measured by acoustic test method on each stratum along the subway in Chengdu .They were used as the regionalized random variables to build attributive geologic model. The established 3D geological model can not only reflect the characteristics of stratum spatial distribution but also intuitively show the dynamic elastic mechanics characteristics of geologic body.It achieved combination and visualization of the distribution of spatial structure and mechanical properties. In addition, the use of kinds of geological statistical analysis methods provided by GOCAD software canpredict any parameter values of arbitrary point in the model and well observe the change trend of acoustic parameters in the study area, which can be simple and efficient to evaluate the soil type and site level along the subway, assisting designers in achieving intuitive and accurate engineering geological information conveniently.

Key words：GOCAD; parametric modeling; 3D geological model; spatial distribution; geological statistics

中圖分類號(hào)：TU443

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1672-7029(2016)03-0523-06

通訊作者：廖昕(1986-)，男，湖南衡陽(yáng)人，博士，從事特殊巖土體工程特性研究；E-mail:xinliao@swjtu.edu.cn

基金項(xiàng)目：國(guó)家自然基金科學(xué)基金資助項(xiàng)目(41172261)

收稿日期：2015-06-19