復雜煤巖體結構動力失穩多參量預報方法研究*

曹建濤，來興平，崔　峰，單鵬飛，楊毅然

(1.西安科技大學 能源學院，陜西 西安 710054；2.教育部 西部礦井開采及災害防治重點實驗室，陜西 西安 710054)

?

復雜煤巖體結構動力失穩多參量預報方法研究*

曹建濤1，2，來興平1，2，崔峰1，2，單鵬飛1，2，楊毅然1，2

(1.西安科技大學 能源學院，陜西 西安 710054；2.教育部 西部礦井開采及災害防治重點實驗室，陜西 西安 710054)

針對煤巖體在開采擾動影響下由于動力載荷的作用將發生復雜的運動演化過程，形成新的更為復雜的煤巖體結構。煤巖體結構所受的應力水平變化特異，煤巖體結構在高應力強卸荷作用下結構發生嚴重畸變，出現煤巖體結構動力失穩、沖擊地壓、大面積片幫和大體積冒頂等動力失穩現象，且煤巖體結構動力失穩始終隨著采動煤巖層的時空演化而發生著動態運動，危險源產生的位置、時間具有不確定性，存在自由面的煤層、采場、巷道、掘進工作面等圍巖結構空間都有發生動力失穩的可能，因此導致了動力失穩危險源的搜索范圍廣泛、定位難度較大，采用單一手段難以滿足監測預報的要求。本論文基于能量積聚和釋放的圍巖結構動力失穩前兆信息研究，利用應力、應變、聲發射、鉆孔窺視等聲、光、電監測手段，通過多參量信息的綜合分析，提出復雜煤巖體結構動力失穩多參量動力失穩預測方法。

復雜條件；煤巖體結構；開采擾動載荷；動力失穩

2.KeyLaboratoryofWesternMineExplorationandHazardPrevention，MinistryofEducation，Xi’an710054，China)

0　引　言

復雜煤巖體結構動力失穩現象發生頻繁，動力災害時有發生。隨采深增加與開采強度的加大，動力災害發生頻次上升，大洪溝礦(采深350 m)累計發生10多次動力失穩破壞現象，烏東煤礦發生動力災害，造成1死2傷；堿溝煤礦和小紅溝煤礦也有類似災害發生，河南義馬煤業集團千秋煤礦發生動力失穩災害，動力作用致使巷道嚴重損壞并冒落，近400 m巷道損毀嚴重，部分巷道完全合攏，作業工人被冒落巷道封堵或掩埋在事發地點，最后造成10人遇難的重大傷亡事故，陜西省榆林市神木縣采空區動力失穩引發了3.0級塌陷地震，造成嚴重的社會影響。

繆協興[1-2]等利用斷裂力學原理，建立了煤巖體滑移裂紋擴展的動力失穩模型。章夢濤[3-5]等提出以動力失穩過程判別準則和能量不穩定判別準則，通過建立數學模型分析煤巖動力失穩問題。齊慶新、劉天泉[6-7]等研究了煤巖動力失穩的黏滑失穩機理，提出了煤巖體結構破壞的“三因素”準則，齊慶新等還研究了煤層動力失穩與巖爆之間的聯系與區別。譚云亮、宋揚等研究了煤巖動力失穩的非吸納性特性。唐春安、潘岳[8-11]等采用突變理論，分析了斷層誘發煤巖結構體動力失穩問題，提出了煤巖體系統失穩彈性能釋放量表達式和破裂臨界條件。唐春安等也采用分析軟件(RFPA)對巖石破壞過程的煤巖動力失穩、巖爆的孕育過程進行了較好的分析。

對于沖擊地壓監測預警、預測方法及設備已有較多成果，例如應力監測[12-13]、微震監測[14-16]、鉆屑量監測[17]、聲發射監測[18-19]、電磁輻射監測[20-22]、震動波CT 監測[23-24]等，但是，這些預測、預警方法多是單套系統進行監測預測或預警，多參量的預測或預警方法現場實踐較少。

基于能量積聚和釋放的圍巖結構動力失穩前兆信息研究，利用應力應變、聲發射、鉆孔窺視等監測手段，利用多參量信息的綜合分析，可以實現復雜煤巖體結構動力失穩多參量動力失穩預測方法。針對應力、應變、聲發射、鉆孔窺視等判定指標需要綜合判定復雜煤巖體結構動力失穩，建立復雜煤巖體結構動力失穩危險性綜合評價函數。

1　復雜煤巖體結構結構特征

圖1　采動作用下動力學破壞情況Fig.1　Deformation and damage of roadway induced by EDZ

論文針對烏東煤礦急傾斜(45°～87°)特厚煤層。采深達400 m.采用水平分段綜放開采，工作面上方采空區不均衡垮落和沉降，局部區域形成載荷集中和不均勻受力，沿垂直剖面煤柱形成“塔形”分布結構，最深部水平分段煤柱載荷產生應力疊加，極易導致煤柱斷裂。由于急斜煤層傾角大，煤巖體在開采擾動作用下產生非對稱性“撬動”作用力，煤巖體結構內部不斷發生著運動演化，煤巖體結構經歷“損傷-變形-斷裂-滑移-失穩”復雜的結構力學變化過程，復雜煤巖體結構發生十分明顯的“剪切-滑移”相對運動，最后由于結構整體強度的劣化發生煤巖體坍塌進而誘發動力失穩災害，對工作面形成動力學破壞。煤巖體能量突然釋放使巷道出現嚴重底臌幫鼓與冒頂，最大幫鼓達1.0 m，最大底鼓達1.2 m；巷道內風筒與錨網撕裂和皮帶側翻(圖1)等。

2　復雜煤巖體結構動態失穩致災時空預測原理

基于能量改變的煤巖體結構動態調控技術原理，以復雜煤巖結構動力失穩防治為目標，通過先進測試方法、綜合理論研究、室內實驗和現場工程試驗等方法與手段，揭示復雜煤巖結構動力失穩規律，建立動壓預報技術，判定動壓防治重點區域，優化采掘布局，減弱或消除動力學失穩危險，制訂適合復雜煤巖結構動力失穩危險預警與管理制度，實現安全開采。

以地面調查、鉆探和超前探查為基礎，建立煤巖體結構空間信息。以室內試驗和現場監測方法，研究分析應力、變形、位移、聲發射等指標參數，研究煤巖體開挖狀態下結構面性狀的變化特征及其演化規律，完成復雜煤巖體動力學破壞與時間密切相關的強度實驗，確定煤巖動力學失穩傾向性及其破壞時聲發射參數、實驗室尺度、不同應力水平下煤巖體內部損傷孕育與發展過程，揭示煤巖體結構特征與力學行為及其對災害調控作用。將復雜煤巖體結構動力失穩多元前兆信息指標進行概率統計分析，應用隸屬函數進行綜合評判，為安全生產管理提供定量依據，建立集監測、預報與管理為一體的高度信息化和智能化安全監控體系。

3　開采擾動動力顯現監測

由于煤巖體結構動力失穩始終隨著采動煤巖層的時空演化而發生著動態運動，危險源產生的位置、時間具有不確定性，存在自由面的煤層、采場、巷道、掘進工作面等圍巖結構空間都有發生動力失穩的可能，因此導致了動力失穩危險源的搜索范圍廣泛、定位難度較大，采用單一手段難以滿足監測預報的要求。本論文基于能量積聚和釋放的圍巖結構動力失穩前兆信息研究，利用應力應變、聲發射、鉆孔窺視等聲、光、電監測手段，通過多參量信息的綜合分析，提出復雜煤巖體結構動力失穩多參量預測方法。

3.1應力監測

工作面順槽內每20 m一組安裝錨桿測力計，在巷道兩幫及頂部各安裝一塊，測力計初始壓力設定為3 MPa.每隔1 h監測1次，并記錄數據，根據數據分析，確定預報值。

圖2　應力監測Fig.2　Stress monitoring

底板巖體的應力大于煤體應力，而煤體部位的應力大于頂板部位的應力值。底板的應力變化是開始時應力上升而后應力趨于平穩，到后期應力發生了急速下降。煤體和頂板部位的應力變化不是十分明顯，基本趨于穩定，煤體部位的應力基本維持在6 MPa左右、頂板煤體的應力基本維持在2 MPa左右，應力變化比較明顯，應力值突然增加到近似原來應力值的2倍，如圖2所示。

3.2變形監測

按照巷道變形監測方案，總結分析巷道寬度、高度受上分層綜采面推進及掘進頭掘進影響的規律，為及時準確提供巷道嚴重變形進行預測預報如圖3，圖4所示。測定的壓力值與變形量與額定壓力值與變形的差值超過正常值1倍以上時，說明監測區域可能具有動力失穩危險性，應立即采取相應動力失穩危險解危措施。

圖3　兩幫相對移近量監測Fig.3　Monitoring of convergence of two sides

圖4　頂板下沉量監測Fig.4　Monitoring of subsidence of roof

3.3結構動力失穩聲發射信號預報

煤巖體損傷破裂聲發射(Acoustic emission)AE信號分析為現場圍巖穩定性預報提供十分有效的信息參數。通過聲發射與煤巖體動態破裂與失穩聯合監測，揭示了動力作用下煤巖破裂與失穩AE大事件和能率特征及規律(圖5)。

圖5　煤巖體破裂與失穩聲發射(AE)特征Fig.5　AE characteristics of instability and failure in coal and rock mass (a)AE大事件　(b)AE能率

如圖5所示，復雜煤巖體破裂失穩分為5個階段：初始壓密階段(ST-1)、破裂擴展階段(ST-2)、加速破裂階段(ST-3)、整體失穩階段(ST-4)和破壞后能量急劇釋放階段(ST-5)，破壞加速累積導致煤巖體整體破壞，AE能率和總事件都達到峰值，之后儲能大量釋放，AE能率和總事件驟然釋放與重構擴展，承載力顯著下降。能率和總事件值分別為4 300和770個/s以上時發出預報信號。

3.4基于鉆孔攝像結構演化內部變形預報

煤巖體局部化破裂孕育演化與失穩致災，實質上是煤巖體局部結構與應力環境作用的動力學演化過程。鉆孔內部煤體出現顯著變形、破壞與坍塌，壓應力和剪應力致使煤巖體產生結構畸變與動態失穩。在鉆孔不同位置和不同深度破壞形式迥異，表明煤巖體內部不同位置局部應力集中程度或應力卸荷度也完全不同，極易在局部區域破壞并誘發動力學失穩。圖6描述了鉆孔內部不同深部破裂與變形結果。當鉆孔內部煤體出現顯著變形、破裂與坍塌，在裂隙邊界點上兩點間產生相對位移，沿鉆孔縱向和橫向發生擴張或移動時發出預報信號。

4　復雜煤巖體結構動力失穩綜合指標預報

復雜煤巖體結構動力失穩現象發生的前兆信息主要表現在，利用應力、應變、聲發射、鉆孔窺視4個判定指標監測“力-震-能”的變化過程，將4個不同量綱判定指標采用隸屬函數進行危險性綜合評判，給出復雜煤巖體結構動力失穩一種簡單的綜合判定方法。

每個指標的評判函數為

fi=a1x1+a2x2+a3x3+a4x4.

(1)

式中影響系數與傾向等級的數值采用現場采集數據集應用專家評判法得出系數具體數值。ai(i=1,2,3)為每個指數的影響系數，分別為0.595，0.601，0.450，0.4；xi(x=1,2,3)為每個因素評判結果，設每個參數對應動力失穩傾向等級Ⅰ，Ⅱ，Ⅲ，Ⅳ 4個類別的分別取值為：0.168，0.333，0.667，1，見表1.

采用二次拋物型隸屬函數，動力失穩傾向性對災害危險性的隸屬函數為

(2)

圖6　煤巖體內部不同深部破裂與變形光學觀測結果Fig.6　Crack and deformation of in-situ coal and rock mass at various depth (a)孔口　(b)1.0 m　(c)2.0 m　(d)3.0 m　(e)4.0 m　(f)5.0 m

類別ⅠⅡⅢⅣ動力失穩傾向無弱中強指標壓力指數/MPa<2[2,5][5,12]>12變形指數/mm<5[5,15][15,50]>50聲發射指數<2500[2500,3000][3000,5000]>5000鉆孔變形指數/mm<1[1,2][2,5]>5

復雜煤巖體結構動力失穩危險性綜合評價函數為

U=c1μ1+c2μ2+c3μ3+c4μ4.

(3)

式中c1，c2，c3，c4為權重，通過原始數據綜合分析和現場實踐經驗總結確定各項指標權重分別為0.25，0.25，0.25和0.25.

表2　復雜煤巖體結構動力失穩危險性分類、名稱

由公式(1)～(3)進行計算，再根據表2可對復雜煤巖體結構動力失穩危險性進行判定。

根據現場應力、應變、聲發射、鉆孔窺視4監測結果進行綜合判定復雜煤巖體結構動力失穩現象。

根據評判函fi=a1x1+a2x2+a3x3+a4x4，得出判定指標對應的函數值，見表3.

表3　判定指標對應評判函數值

根據動力失穩傾向性對災害危險性隸屬函數式(2)得出災害危險性隸屬函數值，見表4.

表4　判定指標對應災害危險性隸屬函數值

根據復雜煤巖體結構動力失穩危險性綜合評價函數

U=c1μ1+c2μ2+c3μ3+c4μ4=0.25×0.878+0.25×0.833+0.25×0.727+0.25×0.973=0.854.

得出復雜煤巖體結構動力失穩危險性綜合評價函數值為0.854，對照表2可得出動力失穩危險性等級為強級。成功預報了動力災害的現象，避免了事故發生。

5　結　論

動力失穩作為一種復雜的圍巖結構動力失穩災害，使得對動力失穩發生機理的研究非常復雜和困難。通過對復雜煤巖體結構形態分析研究，研究復雜煤巖體運動演化、動力作用，利用應力應變監測、聲發射信號分析、鉆孔變形等指標進行多參量綜合預報。

1)采用水平分段綜放開采，工作面上方采空區不均衡垮落和沉降，局部區域形成載荷集中和不均勻受力，最深部水平分段煤柱載荷產生應力疊加，極易導致煤柱斷裂。由于急斜煤層傾角大，煤巖體在開采擾動作用下產生非對稱性“撬動”作用，致使集中載荷突然卸荷，導致煤巖體結構整體失穩，對工作面形成動力學破壞；

2)煤巖體結構內部不斷發生著運動演化，煤巖體結構經歷“損傷-變形-斷裂-滑移-失穩”復雜的結構力學變化過程，由于結構整體強度的劣化發生煤巖體坍塌進而誘發動力失穩災害；

3)基于能量積聚和釋放的圍巖結構動力失穩前兆信息研究，利用應力、應變、聲發射、鉆孔窺視等多元測控技術，利用多參量信息的綜合分析，提出針對復雜煤巖體結構動力失穩的多參量預測方法。

References

[1]繆協興，安里千，翟明華，等.巖(煤)壁中滑移裂紋擴展的沖擊礦壓模型[J].中國礦業大學學報，1999，28(2)：113-117.

MIAO Xie-xing，AN Li-qian，ZHAI Ming-hua，et al.Model of rockburst for extension of slip fracture in palisades[J].Journal of China University of Mining & Technology，1999，28(2)：113-117.

[2]繆協興，孫海，吳志剛.徐州東部軟巖礦區沖擊礦壓機理分析[J].巖石力學與工程學報，1999，18(4)：428-431.

MIAO Xie-xing，SUN Hai，WU Zhi-gang.Mechanism analysis of rockburst in softrock mines in eastern of Xuzhou[J].Chinese Journal of Rock Mechanics and Engineering,1999，18(4)：428-431.

[3]章夢濤.沖擊地壓失穩理論與數值模擬計算[J].巖石力學與工程學報，1987，6(3)：197-204.

ZHANG Meng-tao.Instability theory and mathematical model for coal/rock bursts[J].Chinese Journal of Rock Mechanics and Engineering,1987，6(3)：197-204.

[4]章夢濤，徐曾和，潘一山，等.沖擊地壓和突出的統一失穩理論[J].煤炭學報，1991，16(4)：48-53.

ZHANG Meng-tao，XU Zeng-he，PAN Yi-shan，et al.A united instability theory on coal(rock)burst and outburst[J].Journal of China Coal Society，1991，16(4)：48-53.

[5]章夢濤.我國沖擊地壓預測和防治[J].遼寧工程技術大學學報，2001，20(4)：434-435.

ZHANG Meng-tao.Forecast and prevention of rock burst in China[J].Journal of Liaoning Technical University，2001，20(4)：434-435.

[6]齊慶新，劉天全，史元偉.沖擊地壓的摩擦滑動失穩機理[J].礦山壓力與頂板管理，1995，3(4)：174-177.

QI Qing-xin，LIU Tian-quan，SHI Yuan-wei.Mechanism of friction sliding destability of rock burst[J].Ground Pressure and Strata Control，1995，3(4)：174-177.

[7]齊慶新，史元偉，劉天全.沖擊地壓粘滑失穩機理的實驗研究[J].煤炭學報，1997，22(2)：144-148.

QI Qing-xin，SHI Yuan-wei，LIU Tian-quan.Mechanism of instability caused by viscous sliding in rock burst[J].Journal of China Coal Society，1997，22(2)：144-148.

[8]唐春安.脆性材料破壞過程分析的數值實驗方法[J].力學與實踐，1999，21(2)：21-24.

TANG Chun-an.Numeical tests of progressive failure in brittle materials[J].Mechanics and Engineering，1999，21(2)：21-24.

[9]劉紅元，劉建新，唐春安.采動影響下覆巖垮落過程的數值模擬[J].巖土工程學報，2001，23(2)：201-204.

LIU Hong-yuan，LIU Jian-xin，TANG Chun-an.Numerical simulation of failure process of overburden rock strata caused by mining excavation[J].Chinese Journal of Geotechnical Engineering，2001，23(2)：201-204.

[10]付東波，齊慶新，秦海濤，等.采動應力監測系統的設計[J].煤礦開采，2009，14(6)：13-16.

FU Dong-bo，QI Qing-xin，QIN Hai-tao，et al.Design for mining stress monitoring system[J].Coal Mining Technology，2009，14(6)：13-16．

[11]王平，姜福興，王存文，等.沖擊地壓的應力增量預報方法[J].煤炭學報，2010，35(增1)：5-9.

WANG Ping，JIANG Fu-xing，WANG Cun-wen，et al.The stress incremental forecasting method of rock burst[J].Journal of China Coal Society，2010，35(Suppl.1)：5-9．

[12]陸菜平，竇林名，王耀峰，等.堅硬頂板誘發煤體沖擊破壞的微震效應[J].地球物理學報，2010，53(2)：450-456.

LU Cai-ping，DOU Lin-ming，WANG Yao-feng，et al.Microseismic effect of coal materials rockburst failure induced by hard roof[J].Chinese Journal of Geophysics，2010，53(2)：450-456.

[13]姜福興，楊淑華，成云海，等.煤礦沖擊地壓的微地震監測研究[J].地球物理學報，2006，49(5)：1 511-1 516.

JIANG Fu-xing，YANG Shu-hua，CHENG Yun-hai，et al.A study on microseismic monitoring of rock burst in coal mine[J].Chinese Journal of Geophysics，2006，49(5)：1 511-1 516.

[14]夏永學，康立軍，齊慶新，等.基于微震監測的5個指標及其在沖擊地壓預測中的應用[J].煤炭學報，2010，35(12)：2 011-2 015.

XIA Yong-xue，KANG Li-jun，QI Qing-xin，et al.Fiveindex of microseismic and their application in rock burst forecastion[J].Journal of China Coal Society，2010，35(12)：2011-2015.

[15]Gu S T，Wang C Q，Jiang B Y，et al.Field test of rock burst danger based on drilling pulverized coal parameters[J].Disaster Adv，2012，5(4)：237-240.

[16]齊慶新，李首濱，王淑坤.地音監測技術及其在礦壓監測中的應用研究[J].煤炭學報，1994，19(3)：221-232.

QI Qing-xin，LI Shou-bin，WANG Shu-kun.Application of AE technique in monitoring ground pressure[J].Journal of China Coal Society，1994，19(3)：221-232.

[17]賀虎，竇林名，鞏思園，等.沖擊礦壓的聲發射監測技術研究[J].巖土力學，2011，32(4)：1 262-1 268.

HE Hu，DOU Lin-ming，GONG Si-yuan，et al.Study of acoustic emission monitoring technology for rockburst[J].Rock and Soil Mechanics，2011，32(4)：1 262-1 268.

[18]張偉，來興平，漆濤，等.急斜特厚煤層動力學失穩現場聲發射監測[J].西安科技大學學報，2011，31(3)：253-255.

ZHANG Wei，LAI Xing-ping，QI Tao，etal.Acoustic emission monitoring for dynamic destabilization of heavy and steep coal seam[J].Journal of Xi’an University of Science & Technology，2011，31(3)：253-255.

[19]鄒磊，來興平，王寧波，等.復雜地質條件下急斜厚煤層巷道變形特征[J].西安科技大學學報，2010，30(2)：145-149.

ZOU Lei，LAI Xing-ping，WANG Ning-bo，et al.Roadway deformation characteristics of steep-thick coal seam under complicated geological conditions[J].Journal of Xi’an University of Science & Technology，2010，30(2)：145-149.

[20]竇林名，何學秋，王恩元.沖擊礦壓預測的電磁輻射技術及應用[J].煤炭學報，2004，29(4)：396-399.

DOU Lin-ming，HE Xue-qiu，WANG En-yuan.Electro-magnetic emission technique of monitoring rock burst and its application[J].Journal of China Coal Society，2004，29(4)：396-399.

[21]王書文，毛德兵，杜濤濤，等.基于地震CT 技術的沖擊危險性評價模型[J].煤炭學報，2012，37(增1)：1-6.

WANG Shu-wen，MAO De-bing，DU Tao-tao，et al.Rockburst hazard evaluation model based on seismic CT technology[J].Journal of China Coal Society，2012，37(Suppl.1)：1-6.

[22]竇林名，蔡武，鞏思園，等.沖擊危險性動態預測的震動CT技術研究[J].煤炭學報，2014，39(2)：238-244.

DOU Lin-ming，CAI Wu，GONG Si-yuan，et al.Dynamic risk assessment of rock burst based on the technology of seismic computed tomography detection[J].Journal of China Coal Society，2014，39(2)：238-244.

[23]姜福興，楊淑華，XUN Luo.微地震監測揭示的采場圍巖空間破裂形態[J].煤炭學報，2003，28(4)：357-360.

JIANG Fu-xing，YANG Shu-hua，XUN Luo.Spatial fracturing progresses of surrounding rock masses in longwall face monitored by microseismic monitoring techniques[J].Journal of China Coal Society，2003，28(4)：357-360.

[24]李春睿，康立軍，齊慶新，等.深部巷道圍巖分區破裂與沖擊地壓關系初探[J].煤炭學報，2010，35(2)：185-189.

LI Chun-rui，KANG Li-jun，QI Qing-xin，et al.Probe intorelationship between zonal fracturing and rock burst in deep tunnel[J].Journal of China Coal Society，2010，35(2)：185-189.

Multi-parameter forecasting method of dynamic destabilization of complex coal and rock mass structure

CAO Jian-tao1，2，LAI Xing-ping1，2，CUI Feng1，2，SHAN Peng-fei1，2，YANG Yi-ran1，2

(1.CollegeofEnergyScienceandEngineering，Xi’anUniversityofScienceandTechnology，Xi’an710054，China；

Due to the dynamic loading caused by exploitation，the structure of coal and rock mass has become more complex.Great deformation generated under high stress and strong unloading particularly when encountered large variation of stress level.Phenomenons such as dynamical destabilization of structure，rockburst，wall caving and roof falling etc maybe occur on this account.The dynamic change of dynamic destabilization in structure of coal and rock mass coupled with temporal and spatial evolution of exploitation area，and the location and formation time is uncertain.Locations such as coal seam，stope，roadway and working face etc are very likely to be dangerous zones，resulting that the dynamical destabilization source hard to locate.Therefore，it is unable to meet the requirement of monitoring and prediction with single-parameter forecasting method.This paper employs monitoring measures such as stress-strain，acoustic emission(AE)，borehole imaging and photoelectric，based on the study of message produced before dynamic destabilization by energy accumulation and release.By means of synthetical analysis of multi-parameter message，we have proposed a multi-parameter forecasting method of dynamic destabilization of complex coal and rock mass structure.

complex condition；structure of coal and rock mass；mining disturbance loading；dynamic destabilization

10.13800/j.cnki.xakjdxxb.2016.0301

1672-9315(2016)03-0301-07

2016-01-10責任編輯：劉潔

國家自然科學基金(煤炭聯合基金)重點項目(U13612030);陜西省重點科技創新團隊計劃(2013KCT-16);博士啟動基金(2014QDJ069);博士培育基金(2014051)

曹建濤(1981-),男,山西大同人,博士,講師,E-mail：463583050@qq.com

TD 823

A