巖漿巖侵入煤層的地震屬性識別技術

吳 斌，孟凡彬

巖漿巖侵入煤層的地震屬性識別技術

吳 斌1，孟凡彬2

(1. 中國礦業大學 資源與地球科學學院，江蘇 徐州 221116；2. 中國煤炭地質總局地球物理勘探研究院，河北 涿州 072750)

巖漿巖體嚴重影響井田的煤質、儲量等，對采區工作面布置、巷道掘進及回采工作造成嚴重的影響，因而探測巖漿巖侵入煤層的范圍對煤礦安全高效生產具有重大意義。在總結巖漿巖侵入體的地質特征基礎上，建立巖漿侵入煤層地震地質模型，通過正演，對多種屬性進行交互分析，發現主頻能量屬性與振幅類屬性對巖漿巖反應較為敏感，并且隨著巖漿巖厚度的逐漸增加能量類與振幅類屬性逐漸降低。通過對QN煤礦103工作面高密度三維地震資料連井剖面10煤進行屬性分析，發現頻帶寬度、中心頻率、主振幅屬性可以較好地識別巖漿侵入區。通過對勘探區10煤層進行頻帶寬度、中心頻率和主振幅屬性分析識別巖漿侵入區，其識別結果與鉆孔實際揭露吻合較好，證明采用地震屬性識別技術可以識別煤層巖漿巖侵入區范圍。

巖漿巖；煤層；正演模擬；地震屬性；高密度三維地震

巖漿巖體對含煤地層的侵入破壞，嚴重影響井田的煤質、儲量等[1]，使連續完整的煤層被分割[2-7]，對采區工作面布置、巷道掘進及回采工作造成嚴重的影響。因此，利用地震勘探技術開展采區巖漿巖體預測方法的研究，不僅具有一定的理論意義，而且具有很大的實用價值[8]。

國內外對火成巖的研究比較廣泛，識別技術也具有多種類型，一類是以鉆井、測井資料為主的預測方法[9-14]；一類是以重磁電等資料為主的物探預測方法[15-17]。隨著三維地震勘探的逐步普及，利用地震資料預測火成巖成為主要手段之一[18-24]，三維地震資料具有大面積密集采集信息的優勢，利用地震信息可以從平面和立體角度研究地層的構造、巖性的變化。更進一步的是高密度三維地震資料具有高信噪比、高分辨率、高保真度的特點，在研究地層的構造、巖性的變化上必將顯著提升其效果。

本文以QN煤礦103采區10煤層為例，采用理論分析、正演模型模擬和實際高密度三維地震勘探資料應用相結合的研究方法，系統地開展巖漿巖侵入對煤層的影響區域等基礎研究，在此基礎上，利用地震多屬性分析方法和技術識別巖漿侵入煤層的影響區域。

1 巖漿侵入體的地質特征

1.1 巖漿巖侵入時受圍巖性質的影響

1.2 巖漿巖侵入煤儲層的選層性

巖漿巖侵入含煤地層時傳遞的熱量、壓力和揮發物質足以使煤發生質的變化，而圍巖相對來說變化很小。煤的體積縮小為巖漿侵入提供初期通道。侵入后巖漿進一步吞蝕煤層，使煤層厚度發生明顯變化。因此，在煤與其他巖石相比較時，巖漿首先選擇煤層侵入。另外，巖漿侵入煤層主要取決于煤層的厚度和煤質[1]。

1.3 巖漿巖侵入的方式

巖漿巖侵入的產狀又分為巖床、巖盆、巖蓋、巖脈、巖株和巖基[26]。層狀、樹杈狀、透鏡狀、細脈狀和串珠狀等侵入形狀較為常見，侵入方式一般是順地下某一軟弱地帶或斷層一盤侵入含煤地層[8]。

1.4 研究區的巖漿侵入情況

QN煤礦103采區巖漿侵入較為嚴重，侵入區幾乎覆蓋整個采區，主要侵入10煤層。巖漿巖的侵入破壞了煤層的原生狀態，使煤層焦化，減少了采區的資源儲量；造成煤層變薄或完全被吞噬，不可采區增大，降低了煤層的穩定性。煤層被巖漿巖穿插，出現分叉合并，使煤層夾矸增多，煤層結構趨于復雜；同時巖漿巖沿煤層頂板侵入，使煤層頂板遭受破壞，降低了煤層頂板強度，增加了回采難度。

該區巖漿侵入規律性差，生產揭露情況往往與勘查確定范圍出入較大。經生產揭露，10煤層巖漿巖分布范圍較精查勘探的范圍大，同時巖漿侵入的厚度變化大，在103工作面10煤層回采時發現，巖漿呈樹枝狀沿煤層中部和頂板侵入。

2 巖漿侵入煤層地質地震模型正演

2.1 正演模型建立

通過地震與地質資料的對比可見，由于巖漿巖入侵的厚度不同，有些區域巖漿巖的侵入引起的地震反射波的變化，表現為反射波能量減弱，同相軸凌亂，地震數據品質變差。而有些地區巖漿巖侵入之后幾乎沒有引起煤層反射波的明顯變化，這就給地震資料解釋帶來較大的困難[9]。因此，研究巖漿巖的基本地震響應特征，有利于從地震數據上識別巖漿巖，提高地震資料解釋精度。

將原生煤層、天然焦與巖漿巖設計在同一模型中，進行正演模擬，分析正演模型反射波的地震響應。正演模型如圖1所示，模型寬5 000 m，深600 m，地震子波為主頻50 Hz的雷克子波。模型中的物理參數見表1。

圖1 放大的煤層巖漿巖侵入正演模型

2.2 正演數值模擬

依據地震波在黏彈性介質中傳播地震波場吸收衰減理論對地質模型(圖1)進行正演模擬，對模擬結果進行主振幅、中心頻率、頻帶寬度、主頻能量等多屬性分析，其結果(歸一化后)如圖2所示。

采用交會分析技術對屬性與巖漿侵入厚度、巖漿侵入位置進行研究，對正演數值模擬分析可知：

表1 正演模型物性參數

① 主頻能量與振幅類屬性對巖漿巖反應較為敏感；而頻帶寬度、中心頻率、平均瞬時頻率屬性只有在巖漿巖厚度超過煤層厚度一半時才會有明顯的屬性異常。主要表現為頻帶寬度增高，振幅最大值、均方根振幅、中心頻率、主振幅、平均瞬時振幅降低，主頻能量減弱。

② 隨著巖漿巖的厚度加大能量類與振幅類屬性逐漸降低，當巖漿巖全部吞蝕煤層，這些屬性值又有小幅度的增加。

③ 巖漿巖侵入的位置不同，其造成的屬性差異不同。同等厚度的巖漿巖的振幅類屬性與能量類屬性幅值從大到小的規律為：從煤層頂部侵入時、從煤層底部侵入時、從煤層中部侵入。

圖2 煤層巖漿巖侵入模型正演反射波屬性分析

3 研究區巖漿巖地震屬性解釋

根據正演，地震屬性對圈定巖漿巖分布區具有明顯的效果，現就鉆孔揭露資料與高密度三維地震資料進行分析，進一步探討地震屬性預測巖漿巖分布區的可行性。

3.1 地震屬性特征

以17-18-4、2016-12、18-1、2015-3、2016-16等鉆孔(表2)連井高密度三維地震剖面為例，其中17-18-4、2016-12、2016-16三個鉆孔存在巖漿侵入煤層的情況，17-18-4、2016-12鉆孔巖漿巖從煤層頂部侵入，巖漿巖厚度均為0.4 m，2016-16鉆孔處巖漿巖從煤層中部侵入，巖漿巖厚度為1.7 m。18-1、2015-3孔處沒有巖漿侵入。對該連井剖面10煤層反射波進行屬性分析，時窗選擇為10 ms，其屬性結果如圖3所示。

表2 連井剖面鉆孔煤層厚度統計

從圖3可以看出：

① 對于巖漿巖較厚區域(如2016-16鉆孔)，其地震反射波屬性明顯異于其他區域，表現為頻帶寬度、中心頻率增高，主振幅、平均瞬時振幅降低，主頻能量減弱。

② 對于巖漿巖較薄區域(如17-18-4鉆孔)，其頻帶寬度、中心頻率、主振幅屬性異于無巖漿侵入區，具體表現為巖漿侵入區其地震反射波的頻帶寬度加寬、中心頻率增高、主振幅降低。而平均瞬時振幅、平均瞬時頻率、主頻能量屬性無法分辨出巖漿巖較薄區。

可見頻帶寬度、中心頻率、主振幅屬性可以較好地識別巖漿侵入區，與鉆孔實際情況符合程度高，可以為巖漿侵入區的預測提供參考。因此，利用地震屬性識別巖漿巖具有可行性。

3.2 高密度三維地震屬性識別技術

基于前述地震地質模型正演及研究區10煤層反射波連井剖面多屬性分析成果，10煤層巖漿侵入區與無巖漿侵入區在頻帶寬度、中心頻率、主振幅屬性有相對明顯的區分，巖漿侵入造成10煤層地震反射波頻帶寬度加大、中心頻率增大、主振幅降低。

3.2.1 識別結果

結合研究區地質鉆孔巖漿侵入特點，對研究區高密度三維地震勘探資料中的10煤層地震反射波進行主振幅、頻帶寬度、中心頻率屬性進行提取分析，成果如圖4—圖6所示黑色鉆孔為無巖漿侵入鉆孔，洋紅色鉆孔為存在巖漿侵入的鉆孔，巖漿巖侵入區范圍為已有資料劃定。

圖3 QN煤礦103采區連井剖面10煤層反射波屬性分析

圖4 QN煤礦103采區10煤層反射波主振幅屬性平面圖

圖5 QN煤礦103采區10煤層反射波頻帶寬度屬性平面圖

圖6 QN煤礦103采區10煤層反射波中心頻率屬性平面圖

由圖4—圖6可知：研究區10煤層巖漿巖的分布極其復雜，地震屬性差異也不完全相同。地震屬性圖中顏色用以預測巖漿巖厚度，紅色為厚層，厚度大于1 m；黃色為薄層，厚度小于1 m；藍色為巖漿巖未侵入區。

3.2.2 識別精度分析

屬性預測結果與鉆孔實見情況符合程度見表3，頻帶寬度屬性預測巖漿侵入的準確度為67.9%，中心頻率屬性預測巖漿侵入的準確度為67.9%，主振幅屬性預測巖漿侵入的準確度為60.7%。

巖漿巖厚度識別結果與鉆孔實見情況符合程度見表4，頻帶寬度屬性預測巖漿侵入的準確度為53.85%，中心頻率屬性預測巖漿侵入的準確度為84.62%，主振幅屬性預測巖漿侵入的準確度為53.85%。

由于各類地震屬性與地下地質目標并非有一一對應關系，多種地質異常都可能導致某一屬性出現相同變化，因而，屬性預測與鉆孔實見情況很難出現完全吻合。

表3 鉆孔處巖漿巖地震屬性識別結果與實際揭露吻合情況統計

表4 鉆孔處巖漿巖厚度識別結果與實際揭露吻合情況統計

注：第一種情況為厚度小于1 m，處于橙黃色區域；第二種情況為厚度大于1 m，應處于橙紅色區域。

4 結論

a.通過多屬性對比分析，主頻能量與振幅類屬性相較于頻帶寬度、中心頻率、平均瞬時頻率屬性對巖漿巖反應較為敏感，隨著巖漿巖的厚度加大能量類與振幅類屬性逐漸降低，當巖漿巖全部吞噬煤層，這些屬性值又有小幅度的增加。

b.對于巖漿巖較厚區域，地震反射波頻帶寬度、中心頻率增高，主振幅、平均瞬時振幅衰減，主頻能量減弱。對于巖漿巖較薄區域，巖漿侵入區地震反射波的頻帶寬度加寬、中心頻率增高、主振幅降低。

c.通過高密度三維地震資料提取地震屬性用以預測巖漿巖侵入的效果較好，中心頻率預測巖漿巖厚度的效果最佳。由于各類地震屬性與地下地質目標并非有一一對應關系，多種地質異常都可能導致某一屬性出現相同變化，因而，屬性預測與鉆孔實見情況很難出現完全吻合。

d.可以對多種屬性值配以不同權重因子，合成為總數值加以比較分析，或者配上其他地質異常體的屬性研究統合分析來提高預測精度。

[1] 劉運超，張利轉. 巖漿巖侵入地震識別及解釋方法分析[J]. 煤炭技術，2010，29(-3)：1-7-1–173. LIU Yunchao，ZHANG Lizhuan. Mogma-t-ic rock i-n-trusion seismic identification a-nd interpretation method analysis[J]. C-o-a-l Technology，2010，29(3)：171–173.

[2] 楊孟達. 煤礦地質學[M]. 北京：煤炭工業出版社，2006. YANG Mengda. Coal Mine Geology[M]. Beijing：China Coal Industry Publishi-n-g House，2006.

[3] 張文永，徐勝平，蔡學斌. 臥龍湖煤礦巖漿侵入規律及其對煤層、煤質、瓦斯的影響[J]. 安徽地質，2005，15(1)：25–28. ZHANG Wenyong，XU Shengping，CAI X-uebin. Magamatic intrusion law and i-t-s effect on coal s-e--am，coal quality and g-as in the WOLONGHU c-o-a-l mine [J]. Geology of Anhui，2005，15(1)：2-5–28.

[4] 王瑞貞，白旭明，王金寬，等. 朝克烏拉凹陷火成巖覆蓋區地震勘探方法研究[J]. 石油物探，2020，59(3)：382–395. WANG Ruizhen，BAI Xuming，WANG Jinkuan，et al. Seismic exploration methods for igneous rock-overlain regions in the Chaokewula Sag[J]. Geophysical Prospecting for Petroleum，2020，59(3)：382–395.

[5] 楊梓懂. 呂鳳礦區巖漿侵入體對煤層煤質的影響[J]. 能源與環境，2012(2)：67–69. YANG Zidong. Influence of Magma Int-r-usion on Coal Seam and Coal Quality i-n- Lufeng Mining A-rea[J]. Energy and E-nvironment，2012(2)：67–6-9.

[6] 譚笑林，王兵. 云駕嶺礦擴大區巖漿巖侵入特征及對煤層煤質的影響[J]. 山東煤炭科技，2008(6)：136–137. TIAN Xiaolin，WANG Bing. Features of magma r-ock invading and its influence t-o coal seam qua-l-i-ty in Yunjialing coal mine expending area[J]. S-h-andong C-o-al Science and Technology，2008(6)：1-36–1-3-7.

[7] 劉福勝，徐培武，鄭榮華，等. 邯邢煤田巖漿侵入及對煤層煤質的影響[J]. 中國煤田地質，2007，19(5)：22–23. LIU Fusheng，XU Peiwu，ZHENG Ron-g-h-u-a，et al. Magmatic intrusion and its i-m-pact to coal sea-m and coal quality in H-anxing coal--field[J]. Co-a-l Geolog-y of China，2007，19(5)：22–23.

[8] 朱紅娟. 巖漿巖侵入體三維地震解釋方法及應用研究[D]. 西安：西安科技大學，2-0-0-5. ZHU Hongjuan. The study of interpreta-t-ion met-h-ods 3D seismic and its appli-c--ation on magmati-c intrusive mass[D]. Xi’an：Xi’an University of S-cience and Technology，2005.

[9] 李江，智敏，朱書階. 巖漿巖地震波阻抗反演與厚度預測[J]. 物探與化探，2020，44(5)：1-2-33–1238.LI Jiang，ZHI Min，ZHU Shujie. The f-i-c-t-itious -impedance inversion and thic-k-n-ess prediction of m-agmatic rock[J]. G-e-ophysical and Geochemical E-xploration，2020，44(5)：1233–1238.

[10] 謝慶賓，韓德馨，朱筱敏，等. 三塘湖盆地火成巖儲集空間類型及特征[J]. 石油勘探與開發，2-0-0-2，29(1)：84–86. XIE Qingbin，HAN Dexin，ZHU Xiaomin，et al. Reservoir space feature and evolu-tion of the volc-a-ni-c rocks in the Santanghu basin[J]. Petroleum E-xp-l-oration and Development，2002，29(1)：84–86.

[11] 羅靜蘭，邵紅梅，張成立. 火山巖油氣藏研究方法與勘探技術綜述[J]. 石油學報，2003，24(1)：31–38. LUO Jinglan，SHAO Hongmei，ZHANG Ch-e-ngli. Summary of research methods and ex-plora-t-ion te-c-hnologies for volcanic reservo-i-rs[J]. Acta Petrolei Sinica，2003，24(1)：31–38.

[12] 黃布宙，潘保芝. 松遼盆地北部深層火成巖測井響應特征及巖性劃分[J]. 石油物探，2001，40(3)：42–47. HUANG Buzhou，PAN Baozhi. The log chara-c-t-eristics of deep igneous rocks and lithol-o-g-y determination in the northern Songliao B-asin[J]. Geophysi-c-al Prospecting for Petrol-e--u-m，2001，40(3)：42–47.

[13] 郭秀娟，潘保芝，吳海波，等. 徐家圍子斷陷火成巖裂縫性地層的測井特征與識別[J]. 石油物探，2004，43(3)：301–305. GUO Xiujuan，PAN Baozhi，WU Haibo，et al. Features and identification of fractur-e-d ig-neous r-eservoirs on well logging data:A c-a-se study at X-u-jiaweizi Fault Depressi-o-n[J]. Geophysical Prospecting for Petrol-e--um，2004，43(3)：301–305.

[14] 許風光，鄧少貴，范宜仁，等. 火成巖儲層測井評價進展綜述[J]. 勘探地球物理進展，2006，29(4)：239–243. XU Fengguang，DENG Shaogui，FAN Yire-n，et--- al. Review on logging evaluation of igneous rese-r-voir[J]. Progress in Exploration Geoph-y-sics，2006，29(4)：239–243.

[15] 王家林，王一新，萬明浩. 石油重磁解釋[M]. 北京：石油工業出版社，1991. WANG Jialin，WANG Yixin，WAN Minghao. Pet-r-oleum gravity and magnetic interpretation[M]. B-e-ijing：Petroleum Industry Press，1991.

[16] 鄧榮來，李慶浩，宋桂橋，等. 重磁聯合反演及重磁與MT綜合解釋巴彥浩特盆地火成巖[J]. 石油物探，2002，41(2)：222–225. DENG Ronglai，LI Qinghao，SONG Guiqiao，et al. Investigation on distribution of igneous rock in Bayabhaote Basin with joint inversion and integra-t-ed interpretation of gravity, magnetic and MT da-t-a[J]. Geophysical Prospecting for Petrol-e--u-m，2002，41(2)：222–225.

[17] 楊輝，張研，劉曉，等. 松遼盆地北部徐家圍子斷陷航磁異常特征及油氣地質意義[J]. 石油地球物理勘探，2006，41(1)：97–99. YANG Hui，ZHANG Yan，LIU Xiao，et al. Aer-o-magnetic anomalous feature of Xujiaweizi F-ault Depression in north part of Songliao Basin a-nd geologic meaning for oil/gas exploration[J]. Oil Geophysical Pospecting，2006，41(1)：97–99.

[18] 趙國連，張岳橋. 大慶火山巖地震反射特征與綜合預測技術[J]. 石油勘探與開發，2002，29(5)：44–46. ZHAO Guolian，ZHANG Yueqiao. Siesmic reflecti-on character of volcanic reservoir of Daqing and the comprehensive prediction technology[J]. Petrol-eum Exploration and Development，2002，29(5)：44–46.

[19] 李明，鄒才能，劉曉，等. 松遼盆地北部深層火山巖氣藏識別與預測技術[J]. 石油地球物理勘探，2002，37(5)：477–484.LI Ming，ZOU Caineng，LIU Xiao，et al. I-dent-i-f-ying and predicting technology of de-ep v-olcanics gas reservoir in north of Son-gliao bas in[J]. Oil Geophysical Pospecting，2002，37(5)：477–484.

[20] 高印軍，郭春東，余忠，等. 綜合解釋技術在孔南地區火成巖油藏研究中的應用[J]. 石油勘探與開發，2000，27(3)：81–83. GAO Yinjun，GUO Chundong，YU Zhong，e-t al. Applying comprehensive interpretation techni-que to the study of igneous rock reservoir in Ko-n-gnan area[J]. Petrol-eum Exploration and Develo-p-ment，2000，27(3)：81–83.

[21] 崔勇，欒瑞樂，趙澄林. 遼河油田歐利坨子地區火山巖儲集層特征及有利儲集層預測[J]. 石油勘探與開發，2000，27(5)：47–49. CUI Yong，LUAN Ruile，ZHAO Chenglin. Volca-n-ic reservoir characteristics and favorable reservoir prediction in Oulituozi area of Liaohe Oilfield[J]. Petrol-eum Exploration and Develo-p-ment，2000，27(5)：47–49.

[22] 呂希學，肖煥欽，劉惠民，等. 根據地震勘查信息預測火成巖儲集體[J]. 浙江大學學報(工學版)，2003，37(2)：123–127. LYU Xixue，XIAO Huanqin，LIU Huimin，et al. P-rediction of the igneous rock reservoir by using multi-seismic information[J]. Jou-r-nal of Zhejiang University(Engineering S-cience)，2003，37(2)：123–1-27.

[23] 吳海波，董守華，黃亞平，等. 基于地震屬性的煤層火成巖侵入預測[J]. 地球物理學進展，2015，30(3)：1376–1381. WU Haibo，DONG Shouhua，HUANG Yaping，e-t al. Prediction of coal seam igneous intrusion ba-sed on seismic attributes[J]. Progress in Geophysi-c-s，2015，30(3)：1376–1381.

[24] 李仁海，崔若飛，毛欣榮，等. 利用巖性解釋方法圈定巖漿巖侵入煤層范圍[J]. 地球物理學進展，2008，23(1)：242–248. LI Renhai，CUI Ruofei，MAO Xinrong，et al. Determination of magnetic rock in-t-rusion using li-t-hology interpretation met-h-od[J]. Progress in Geophysics，2008，23(1)：242–248.

[25] 金法禮，冀明君，張培礎. 淮北煤田永固井田巖漿巖侵入特征及煤厚預測[J]. 中國礦業大學學報，1998，27(2)：209–212. JIN Fali，JI Mingjun，ZHANG Peichu. M-agmati-c rock intrusive characteristics a-nd coal thickness p-rediction in Yonggu minefield of Huaibei coalfie-l-d[J]. Jo-ur-nal of China University of Mi-ning & Technology，1998，27(2)：209–212.

[26] 樂昌碩. 高等學校教材巖石學[M]. 北京：地質出版社，1984. LE Changshuo. Teaching materials of colleges and Universities Petrology[M]. Beijing：Geological Pu-blishing House，1984.

Seismic attribute recognition technology of magmatic rock intrusive coal seam

WU Bin1, MENG Fanbin2

(1.School of Resources and Geosciences, China University of Mining and Technology, Xuzhou 221116, China; 2.Research Institute of Coal Geophysical Exploration, China National Administration of Coal Geology, Zhuozhou 072750, China)

Magmatic rock seriously affects the coal quality and reserves of the mine fields, has serious impact on the mining face layout, roadway excavation and mining work. Therefore, it is of great significance to detect the scope of magmatic rock intrusion into coal seam for safe and efficient production of coal mine. This paper summarizes the geological characteristics of magmatic intrusions, establishes the geological and seismic model of magmatic intrusion into coal seam. Through forward modeling, it is found that the main frequency energy attributes and amplitude attributes are more sensitive to magmatic rock reaction, and with the gradual increase of magmatic rock thickness, the energy and amplitude attributes gradually decrease. Through the attribute analysis of coal seam 10 in the actual high-density 3D seismic data, it is found that the frequency band width, central frequency and main amplitude attributes can better identify the magmatic intrusion area. Through the analysis of frequency band width, center frequency and main amplitude of coal seam 10 in the exploration area, it is found that the identification of magmatic intrusion area is in good agreement with borehole exposure, which proves that the seismic attribute identification technology can identify the magmatic rock intrusion area of coal seam.

magmatic rocks; coal seam; forward modeling; seismic attributes; high density 3D seismic exploration

請聽作者語音介紹創新技術成果等信息，歡迎與作者進行交流

P631.42

A

10.3969/j.issn.1001-1986.2020.06.009

1001-1986(2020)06-0064-08

2020-10-15；

2020-11-03

國家重點研發計劃課題(2018YFC0807803)

National Key R&D Program of China(2018YFC0807803)

吳斌，1996年生，男，安徽宣城人，碩士研究生，從事煤田地震勘探理論與方法研究. E-mail：2016372561@qq.com

吳斌，孟凡彬. 巖漿巖侵入煤層的地震屬性識別技術[J]. 煤田地質與勘探，2020，48(6)：64–71.

WU Bin，MENG Fanbin. Seismic attribute recognition technology of magmatic rock intrusive coal seam[J]. Coal Geology & Exploration，2020，48(6)：64–71.

(責任編輯 聶愛蘭)