劉紅衛(wèi), 謝興兵, 徐連三, 周 磊, 嚴(yán)良俊,徐聯(lián)澤,5, 鄒立帆, 徐 楊, 王富強, 周 倩
(1.湖北省地質(zhì)局 武漢水文地質(zhì)工程地質(zhì)大隊,湖北 武漢 430051; 2.湖北省地?zé)崮苎芯客茝V中心,湖北 武漢 430051;3.長江大學(xué) 地球物理與石油資源學(xué)院,湖北 武漢 430100; 4.油氣資源與勘探技術(shù)教育部重點實驗室,湖北 武漢 430100;5.武漢地質(zhì)工程勘察院有限公司,湖北 武漢 430051)
地球物理勘探是探尋深地特征、找油氣找礦找熱的重要手段,電磁法是其中最主要的方法之一[1]。近年來,電磁法勘探理論、技術(shù)創(chuàng)新發(fā)展迅速,中國在地面電磁法、地空電磁法、海洋電磁法、井中電磁法和航空電磁法等領(lǐng)域的研究取得了突破性進(jìn)展,研發(fā)了一批應(yīng)用于海陸空不同環(huán)境下的大功率、新型成套探測裝備和軟件系統(tǒng),并在應(yīng)用中取得顯著成效[2]。在中深層地?zé)峥碧椒矫妫S著中深層地?zé)峥碧綄ι疃取⒕纫蟮牟粩嘣鲩L,頻率域電磁法(MT/AMT/CSAMT)、廣域電磁法(WFEM)、大功率電性源時間域電磁法(LOTEM)、電性源短偏移距瞬變電磁法(SOTEM)、時頻電磁法(TFEM)等電磁勘探方法得到廣泛研發(fā)應(yīng)用,并向大功率、多手段、精細(xì)探測與解譯方向拓展[3-5]。盛勇等[6]在大別山某谷地采用微動結(jié)合放射性氡氣(Rn)測量、剖面、高密度電阻率測深、激電測深、瞬變電磁測深(TEM)、音頻大地電磁測深(AMT)等7種方法,最終圈定了熱儲位置并定孔驗證,成功打出了熱水。曾何勝等[7]利用大功率廣域電磁法,在復(fù)雜電磁干擾區(qū)大深度地?zé)豳Y源勘查中取得良好的效果。譚章坤等[8]結(jié)合實例說明了廣域電磁法在深部地?zé)峥碧街械挠行院蛢?yōu)越性。袁剛等[9]利用高精度大地電磁方法在昭通市低溫?zé)崴M(jìn)行了連續(xù)電導(dǎo)率成像法勘探,經(jīng)鉆孔鉆探驗證,證明了該方法的有效性。
除盆地層狀地?zé)嵬猓芏鄮r漿巖、變質(zhì)巖、深覆蓋碳酸鹽巖區(qū)地?zé)岢蕩罨蚓植繉訝顝?fù)合帶狀方式賦存,地?zé)崴畠Υ婵臻g、導(dǎo)流通道為不規(guī)則網(wǎng)絡(luò)型式,同一儲水層段、同一條導(dǎo)水導(dǎo)熱斷裂在橫向、縱向上富水性和導(dǎo)流能力差別大,亟需開展精細(xì)勘探:既要找到深部富水層段、導(dǎo)水導(dǎo)熱通道,又要找到其中的優(yōu)勢滲流通道,以指導(dǎo)地?zé)徙@井到達(dá)深部有利“甜點”目標(biāo)靶區(qū)。這個精度有時需要達(dá)到深度3 000 m以淺、異常塊段尺寸為50~100 m,甚至深度更深、異常塊段尺寸更小的精度。受經(jīng)濟(jì)、技術(shù)條件的制約,目前大多數(shù)深層地?zé)峥碧诫y以達(dá)到該精度,這是地?zé)峥碧街兄档藐P(guān)注的問題,故利用現(xiàn)有技術(shù)和裝備,試驗探索大深度、厚覆蓋、強電磁干擾區(qū)地?zé)峋?xì)勘探技術(shù),具有較強的實用性和指導(dǎo)意義。
本文結(jié)合武漢某地?zé)峥碧窖芯宽椖浚诘刭|(zhì)調(diào)查、廣域電磁探測基礎(chǔ)上,采用大功率電性源時間域電磁法(LOTEM)開展中深層地?zé)峋?xì)探測技術(shù)試驗研究,主要對該方法在中深層地?zé)峋?xì)探測應(yīng)用方面的可行性,數(shù)據(jù)采集、處理和解譯等進(jìn)行試驗、分析和討論,利用地層(地?zé)醿Α?dǎo)層段)電阻率和極化率差異,獲取研究區(qū)試驗勘探測線的深度—電阻率剖面和極化率剖面,并進(jìn)行綜合解譯,大致勾畫出深部地?zé)岙惓^(qū)和優(yōu)勢靶點位置,指導(dǎo)深層地?zé)峥碧竭x區(qū)選點和鉆探設(shè)計。
大功率電性源時間域電磁法以長接地導(dǎo)線(2~5 km)向地下發(fā)送大功率長時窗雙極性占空比為1∶1方波,在一次場斷電后,采用類地震多次疊加技術(shù)觀測感應(yīng)二次場(?Bz/?t、Ex、Ey)隨時間的衰減,地下感應(yīng)二次場的強弱、隨時間衰減的快慢與地下所要探測的地質(zhì)異常體的大小、位置、產(chǎn)狀和導(dǎo)電性能等因素密切相關(guān)。該方法具有穿透高阻層能力強、對低阻靈敏、人工源多次疊加技術(shù)能壓制隨機干擾、大功率長時窗全波形觀測勘探深度大(可達(dá)5 km)及觀測精度高等優(yōu)點,在遼河油田強干擾區(qū)地?zé)峥辈椤⒔瓭h油田潭口油田剩余油檢測等項目中應(yīng)用效果良好[10-12]。大功率電性源時間域電磁法野外勘探布設(shè)示意見圖1,其發(fā)射波形及二次場響應(yīng)特征見圖2。


圖1 野外勘探布設(shè)示意圖Fig.1 Schematic diagram of field deployment

a.發(fā)射波形;b.高低電阻率的二次場響應(yīng)衰減曲線圖2 發(fā)射波形及不同電阻率的二次場響應(yīng)特征Fig.2 Emission waveform and secondary field response characteristics of different resistivity
區(qū)內(nèi)褶皺主要呈近EW向展布,中部為背斜,向東傾伏;北側(cè)、南側(cè)為向斜。褶皺多被斷層錯斷。
區(qū)內(nèi)斷裂主要發(fā)育有三個方向:NE、NW和近EW向,斷層性質(zhì)以逆斷層、壓扭斷層為主,斷裂規(guī)模、性質(zhì)及組合關(guān)系較復(fù)雜。
根據(jù)研究區(qū)周邊已知地?zé)峋饕貙拥碾娮杪始皫r性(表1),建立如圖3所示的層狀地電模型。

表1 地?zé)峋饕貙与娮杪始皫r性Table 1 Main formation resistivity and lithology of geothermal well

圖3 層狀地電模型圖Fig.3 Layered geoelectric model diagram
研究區(qū)及其周邊地層三疊系、石炭系、奧陶系上統(tǒng)—寒武系上統(tǒng)以灰?guī)r、白云巖為主,呈高阻,電阻率在500~4 000 Ω·m;泥盆系以石英砂巖、礫巖為主,呈中低阻,電阻率在200~500 Ω·m;志留系、二疊系以砂巖、粉砂巖、泥巖、炭質(zhì)泥巖、炭質(zhì)灰?guī)r為主,呈低阻,電阻率在20~200 Ω·m;破碎帶電阻率在50~260 Ω·m。推測炭質(zhì)層和深部含水段極化率較高。
據(jù)區(qū)域及鄰近資料,研究區(qū)具備中深層地?zé)豳x存條件[14],地?zé)崴难a給源頭主要為周邊中低山區(qū)大氣降水,降水入滲后向深部徑流,沿斷裂、深部巖溶系統(tǒng)運移,并賦存于地下深部巖石、斷裂和巖溶空隙(管道)中,最終沿區(qū)域深大斷裂、盆地排泄。由于研究區(qū)及周邊大部分地區(qū)上部被近千米厚的泥巖等地層覆蓋,深部熱水不能直接出露地表,需要鉆鑿深井揭露。
根據(jù)圖3建立的層狀地電模型和研究區(qū)中深層地?zé)豳x存條件,通過Blender和Tetgen軟件建立起伏地形地質(zhì)—地電模型(圖4)。其中發(fā)射沿著y方向,發(fā)射電流100 A,發(fā)射偶極矩長2 km,偏移距2~4 km,低阻異常體位于地下-1 400 m且沿著一定的傾向延伸到地下-1 650 m,長約300 m、寬300 m、厚250 m,電阻率為50 Ω·m,并基于非結(jié)構(gòu)矢量有限元法三維正演模擬技術(shù)[15],對該地電模型進(jìn)行了三維網(wǎng)格剖分和正演模擬(圖5),得到不同偏移距電場響應(yīng)圖(圖6)和相對異常響應(yīng)圖(圖7)。

圖4 三維地電模型圖Fig.4 Three-dimensional geoelectric model

圖5 三維地電模型網(wǎng)格剖分圖Fig.5 Grid division of three-dimensional geoelectric model
根據(jù)圖6可知,在現(xiàn)有發(fā)射條件和地電條件下,不同收發(fā)距時電場的二次衰減場值對于層狀、連續(xù)含水層段的響應(yīng)明顯,但對于含水層段中低阻異常體難以識別,而用電場相對異常響應(yīng)能清晰地反映低阻異常體的位置(深度)(圖7),說明研究區(qū)采集電場分量用于解譯和識別地?zé)醿蛹捌渲袃?yōu)勢滲流通道的空間位置、分布是可行的。
試驗測線(SY01)垂直于主要地質(zhì)構(gòu)造走向布設(shè),與先期探測的廣域電磁測線(GY01)基本平行,測線長度為2 475 m,觀測點點距為50 m(局部加密至25 m),發(fā)射源(與測線垂直)長度為2 400 m,偏移距為1.6~4.0 km(圖8)。

(左為無低阻異常體;右為有低阻異常體)圖6 不同偏移距電場Ey響應(yīng)圖Fig.6 Electric field Ey response with different offset

圖7 不同偏移距低阻異常體電場相對異常響應(yīng)圖Fig.7 Electric field relative anomaly response of low resistance geophysical anomaly with different offset

圖8 測線位置示意圖Fig.8 Location of measuring line
大功率電性源時間域電磁法數(shù)據(jù)采集方式如圖1所示,采取徑向方式,即測線垂直于發(fā)射偶極,接收采集正交的Ex、Ey和Bz三個分量,接收點距50 m,核心區(qū)點距25 m,接收偶極矩長50 m;發(fā)射采用大功率SuperTx-150發(fā)射系統(tǒng),最大發(fā)射功率200 kW,最大發(fā)射電流100 A,發(fā)射波形為占空比1∶1雙極性方波(TD50)(周期為4 s)。接收系統(tǒng)采用加拿大鳳凰地球物理公司多功能電法采集V8系統(tǒng),測點定位采用天寶R8高精度GPS測量設(shè)備。
大功率電性源時間域電磁法數(shù)據(jù)處理采用UTMPro軟件,UTMPro軟件為自主研發(fā)的瞬變電磁數(shù)據(jù)處理軟件,主要包括數(shù)據(jù)管理、處理、繪圖、輸出等功能模塊。數(shù)據(jù)管理模塊主要包括原始時間序列數(shù)據(jù)的參數(shù)和波形瀏覽、發(fā)射和采集參數(shù)校正等功能;數(shù)據(jù)處理模塊是UTMPro軟件的核心,功能主要包括數(shù)據(jù)疊前對時和方向校正、切割、濾波、統(tǒng)計功能,數(shù)據(jù)的多周期疊加、雙極性波形的正負(fù)疊加功能,數(shù)據(jù)疊后的濾波、加窗功能,及疊前和疊后數(shù)據(jù)的傅里葉變換等;數(shù)據(jù)繪圖模塊主要包括原始波形的成圖、疊加數(shù)據(jù)的成圖、傅里葉頻譜成圖等;數(shù)據(jù)輸出模塊可以自定義,按照不同的格式輸出疊前和疊后不同處理階段的數(shù)據(jù)。
研究區(qū)大功率電性源時間域電磁法點位周邊有公路、高壓電網(wǎng)和變壓器、變電站、湖泊、村鎮(zhèn)、大型娛樂設(shè)施等,由于電磁干擾嚴(yán)重,為了獲取合格的野外數(shù)據(jù),針對部分測點的環(huán)境,采取了延長采集時間、按規(guī)范偏移測點位置等措施。測量期間,布設(shè)質(zhì)量檢查點3個,各檢查點均方相對誤差平均值均<5%,滿足行業(yè)標(biāo)準(zhǔn)要求。
為對比和定性分析每個測點原始衰減曲線的變化情況,繪制了所有測點的蒙太奇圖(圖9)和擬斷面圖(圖10)。

圖9 所有測點歸一化電場衰減曲線蒙太奇圖Fig.9 Montage diagram of normalized electric field attenuationcurve of all measuring points
由于大功率發(fā)射系統(tǒng)SuperTx-150發(fā)射一次場信號強(發(fā)射電流>100 A),既保證了接收信號的高信噪比,也保證了后期資料處理的可靠性。通過圖9-圖10可知,所有測點數(shù)據(jù)的二次衰減場在500 ms內(nèi)均達(dá)到優(yōu)秀級別,且具有較好的一致性。
大功率電性源時間域電磁法的反演方法很多[16],本文采用降維單參數(shù)法[17](類同于MT反演方法)。降維單參數(shù)法是根據(jù)一維正演計算的視電阻率理論(擬合)值與實測值差別依次調(diào)整每一個電性層的厚度、埋深和電阻率,循環(huán)進(jìn)行直至“虛擬全區(qū)視電阻率”的實測值與正演計算的理論(擬合)值在所有時間點上的擬合總體相對誤差小于預(yù)先指定的一個很小的數(shù)。為了消除淺部地形和不均勻異常體靜態(tài)偏移的影響,通過研究區(qū)時域垂直磁場進(jìn)行了二次衰減場響應(yīng)曲線的校正,保證了反演初始數(shù)據(jù)的可靠性。
依據(jù)上述反演步驟和方法,得到研究區(qū)一維反演視電阻率剖面(圖11)、二維反演視電阻率剖面(圖12)、一維反演視極化率剖面(圖13)。

圖10 所有測點歸一化電場衰減擬斷面圖Fig.10 Quasi section of normalized electric field attenuation of all measuring points

圖11 一維反演視電阻率剖面圖Fig.11 One-dimensional inversion resistivity profile

圖12 二維反演視電阻率剖面圖Fig.12 Two-dimensional inversion resistivity profile

圖13 一維反演視極化率剖面圖Fig.13 One-dimensional inversion polarization profile
綜合分析一維、二維反演視電阻率剖面和一維反演視極化率剖面可知:
(1) 圖11-圖12反演視電阻率剖面均反映淺部(地表以下-600~-1 200 m)電阻率為低阻,電阻率<200 Ω·m;下伏地層在測線南北兩端區(qū)域表現(xiàn)為中—高阻;測線中部有明顯的似“V”形、寬度較大、較陡立的低阻異常帶。
(2) 一維反演視電阻率剖面解譯結(jié)果反映在測線1 700~2 400 m段、深度-500~-800 m范圍存在水平狀、視電阻率為10 Ω·m的低阻帶,推測為志留系底部和龍馬溪組炭質(zhì)頁巖(異常①區(qū));測線1 400~1 900 m段、深度-600~-1 600 m的范圍存在一南傾、視電阻率<10 Ω·m的低阻帶(異常②區(qū)),推測為破碎帶(含水帶),該異常段特征與圖10結(jié)果基本一致;在該異常段的南側(cè)、測線900~1 400 m段中部及1 800~2 000 m段深部(異常③區(qū)),亦有相對異常反映。
(3) 二維反演視電阻率剖面解譯結(jié)果反映在測線中部存在一個較大的低阻異常段,總體南傾,異常段由三個異常區(qū)組成:中部偏北的異常區(qū)規(guī)模較小,分布在測線1 500~1 700 m段、深度-1 000~-1 400 m的位置(異常②區(qū)),有向深部延伸的趨勢;主異常區(qū)位于測線900~1 400 m段、深度-1 400~-1 900 m的位置(異常③區(qū));中部異常區(qū)分布在測線900~1 400 m段、深度-800~-1 400 m的位置(異常④區(qū)),亦呈現(xiàn)低阻、相對破碎特征,異常體總體傾向與北邊基本一致,但未直接連通。
(4) 一維和二維解譯結(jié)果在測線北段深部高電阻率層頂面深度相差較大(相差400~500 m),且北邊異常段的位置、深度亦略有差別。
(5) 視極化率異常主要反映在測線北段1 600~2 450 m(深度-800 m及以下),推測其中2 100~2 450 m 段(異常①區(qū))高極化為炭質(zhì)層引起;1 600~2 100 m 段異常區(qū)(異常②區(qū))可能為含較高礦化度水所致,該異常總體向南側(cè)深部延伸,與該段一維解譯結(jié)果北邊異常基本吻合,但延伸的深度較難確定。
結(jié)合地質(zhì)構(gòu)造、地層、地?zé)豳x存“源、蓋、儲、導(dǎo)”等要素綜合分析,研究區(qū)中段具備中深層地?zé)豳x存條件,有較好找熱前景。
(1) 地質(zhì)構(gòu)造。研究區(qū)地質(zhì)構(gòu)造較為復(fù)雜,以褶皺為主,斷層較為發(fā)育。受多期次地質(zhì)活動影響,斷層以由北向南的逆沖為主,兼具對沖(反沖)特征,同時存在由南向北(早期、深部)和由北向南(后期、淺部)逆沖運動,并受新華夏系剪切走滑運動改造,整體在近EW向主體構(gòu)造格架下被NE向和NNW向斷層切割,淺部發(fā)育多條以逆斷層性質(zhì)為主的順層淺表斷層。研究區(qū)主體為背斜,南北兩側(cè)為向斜,褶皺軸線近EW向,總體向東傾伏,軸面略南傾。測線所處位置基本為背斜核部,據(jù)物探結(jié)果并結(jié)合地質(zhì)調(diào)查分析,推測研究區(qū)中部有一條規(guī)模較大斷層通過,斷層走向近EW向,南傾,并與鄰近多條規(guī)模不等的NNW、NE向斷層交會。
(2) 地?zé)豳x存條件。從圖11-圖12反演視電阻率剖面明顯看出,剖面中部存在較為復(fù)雜的斷裂體系,高阻層可解譯為灰?guī)r、白云巖層,低阻層應(yīng)為泥巖、砂巖層和第四系沉積蓋層。在測線中段高阻層段中的低阻異常應(yīng)是斷層破碎帶的反映,該段電阻率在10~50 Ω·m,異常范圍較大,應(yīng)該是含水(層、帶)的體現(xiàn)。分析視極化率剖面,向深部延伸的高極化帶含水可能性較大。
從地質(zhì)構(gòu)造、地層、地?zé)豳x存條件綜合分析,淺部-900~-1 200 m深度內(nèi)主要為三疊系—志留系地層,為全覆蓋蓋層;-1 000~-2 000 m深度內(nèi)主要為寒武系—奧陶系灰?guī)r、白云巖,因地處背斜核部,節(jié)理裂隙較發(fā)育,有利于巖溶、風(fēng)化作用發(fā)育,加之多期次、多方向、不同性質(zhì)斷裂的影響,有利于與外界連接并形成儲水、導(dǎo)水空間和網(wǎng)絡(luò);研究區(qū)恒溫帶溫度為18℃,地溫梯度在2℃/100 m左右,推測地下-1 400~-2 000 m深度地溫約為46~58℃,總體具備找熱條件。
研究區(qū)地?zé)豳x存條件及地?zé)徙@井綜合分析見圖14,其底圖為廣域電磁法測線GY01(長度2 700 m)二維解譯成果圖。從圖14上看,廣域電磁法測線探查到在測線中部有一寬約600 m、深度在-1 100~-2 000 m、傾向SW的低阻異常區(qū),推測為斷裂及受斷裂影響富水有利區(qū),但是該異常區(qū)范圍大,導(dǎo)水、富水主要空間和通道不清晰。綜合廣域電磁法和大功率電性源時間域電磁法解譯成果分析,可以得到較為精細(xì)的異常區(qū)特征:
① 有利(可靠)區(qū)段為異常③區(qū),它位于廣域電磁法GY01測線約1 175~1 400 m段(SY01線1 000~1 225 m段)附近,其規(guī)模較大、連通性較好、深度及地溫適中,兼具鉆井施工成本、潛在出水量、出水溫度優(yōu)勢。該異常區(qū)內(nèi)深部優(yōu)勢儲、導(dǎo)水靶點為B點(點位:測線中部約1 450 m段處,深度約-1 600 m,靶點范圍大小約50~100 m、SW傾向)和A點(點位:測線中部約1 225 m段處,深度約-1 800 m,靶點范圍大小約100~200 m、SW傾向)。在地面用地、用熱、鉆探施工作業(yè)條件符合時,可考慮布設(shè)ZK1、ZK2鉆井。實際設(shè)計施工時,還可根據(jù)地面、鉆探揭露地質(zhì)及地層、抽水水量、水溫等情況,設(shè)計、施工由B點向A點的定向井。
② 較有利區(qū)段為異常②區(qū),它位于廣域電磁法GY01測線約1 675~1 875 m段(SY01測線1 500~1 700 m段)附近,其埋深較淺,推測與異常③區(qū)有一定連通。該異常區(qū)內(nèi)深部優(yōu)勢儲、導(dǎo)水靶點為C點(點位:測線中部約1 750~1 825 m段處,深度-1 200~-1 400 m,靶點范圍大小約50~100 m、主體SW傾向)。在地面用地、用熱、鉆探施工作業(yè)條件符合時,可考慮布設(shè)ZK3鉆井。實際設(shè)計施工時,還可根據(jù)地面、鉆探揭露地質(zhì)及地層、抽水水量、水溫等情況,設(shè)計、施工由②區(qū)C點向異常③區(qū)B、A靶點的定向井。

圖14 研究區(qū)地?zé)豳x存條件及地?zé)徙@井綜合分析圖Fig.14 Geothermal occurrence conditions and comprehensive analysis of geothermal drilling in the study area
由于探測手段、電磁干擾及作業(yè)環(huán)境、地質(zhì)及地球物理條件等的影響,在探測試驗和解譯、綜合分析過程中尚存在較多問題需要進(jìn)一步試驗、探討:
(1) 大功率電性源時間域電磁法解譯剖面淺部細(xì)節(jié)縱橫向分辨率不高,分析其原因是因為發(fā)射波形為TD50-4s雙極性方波,電偶極矩大(2.4 km),發(fā)射功率大,受關(guān)斷時間的影響,二次衰減場第一個時間點為10 ms,淺部存在約500 m的盲區(qū)。
(2) 大功率電性源時間域電磁法測線大偏移距一側(cè)(南側(cè))解譯結(jié)果與廣域電磁法解譯結(jié)果存在較大差異,可能與探測時發(fā)射源、測線的布置方式有關(guān),因為在平行現(xiàn)實際試驗測線的東西兩側(cè),合適的收發(fā)距范圍內(nèi),適宜布設(shè)與勘探線平行的電偶極發(fā)射源的地方均遇湖水,試驗探測時只能將大功率電性源時間域電磁法電偶極發(fā)射源方向布設(shè)為與測線方向垂直,且該發(fā)射場所發(fā)射電磁場主方向與主地質(zhì)構(gòu)造線方向一致,對結(jié)果有影響。建議在具備條件的地方,再采用平行測線布置方式進(jìn)行探測、對比。
(3) 采用大功率電性源時間域電磁法探測一維、二維解譯成果顯示,北部異常的位置、深度存在一定差別,產(chǎn)生的原因可能是靜態(tài)位移的影響。
(4) 本次使用大功率電性源時間域電磁法探測雖然發(fā)射電流值較高,但在武漢地區(qū)強干擾、淺層為厚層低阻環(huán)境下,探測深度似乎仍難超過2 000 m。
(5) 研究區(qū)鄰近城區(qū)和多個變電站,前期電磁環(huán)境探測時實測只有在凌晨1:30-3:30時為基本無強干擾物探探測窗口,本次大功率電性源時間域電磁法探測數(shù)據(jù)采集均為白天進(jìn)行,對探測結(jié)果是否有影響及影響程度大小需要進(jìn)一步分析并引起注意。
(6) 極化探測可影響深度和極化率解譯剖面適用推斷深度值得進(jìn)一步探討。
(7) 研究區(qū)地?zé)豳x存條件和特征尚需進(jìn)一步勘查、驗證。
致謝:感謝武漢地質(zhì)工程勘察院有限公司,中國地質(zhì)大學(xué)(武漢)姚春亮,中南大學(xué)李帝銓、朱云起,課題組成員及有關(guān)領(lǐng)導(dǎo)、專家對本項目的支持、指導(dǎo)!