圍壓條件下砂巖循環(huán)沖擊損傷的力學(xué)與超聲分析

劉少赫1，許金余1,2，王鵬1，張國喜3

(1.空軍工程大學(xué)機場建筑工程系，西安710038； 2. 西北工業(yè)大學(xué)力學(xué)與土木建筑學(xué)院，西安710072；3.空軍航空大學(xué)飛行基礎(chǔ)訓(xùn)練基地后勤部, 長春130022)

摘要：采用帶圍壓裝置的φ100 mm霍普金森壓桿系統(tǒng)，對砂巖試樣進行了單軸沖擊和4 MPa、20 MPa圍壓條件下的循環(huán)沖擊試驗，并對每次沖擊前后的砂巖試樣進行超聲縱波檢測，分析了砂巖在沖擊荷載循環(huán)作用下的應(yīng)力、應(yīng)變特征，定義了砂巖試樣的屈服-彈性比用以描述圍壓下試樣應(yīng)力應(yīng)變曲線的彈塑性特征，采用縱波波速定義了砂巖試樣的沖擊損傷并分析了循環(huán)沖擊試驗中砂巖試樣縱波波速和應(yīng)力應(yīng)變之間的關(guān)系。分析發(fā)現(xiàn)：在圍壓作用下，應(yīng)力-應(yīng)變曲線呈現(xiàn)典型的彈塑性特征。隨著沖擊荷載循環(huán)作用次數(shù)的增加，砂巖試樣屈服應(yīng)力、峰值應(yīng)力降低，屈服應(yīng)變、峰值應(yīng)變增加。隨著循環(huán)沖擊作用次數(shù)的增加，砂巖試樣的應(yīng)力、應(yīng)變特征與縱波波速之間存在良好的相關(guān)關(guān)系。較低圍壓狀態(tài)下累計損傷度明顯高于較高圍壓下砂巖的累積損傷，砂巖循環(huán)沖擊損傷具有明顯的圍壓效應(yīng)。研究結(jié)果對地下工程的建設(shè)和防護有一定的指導(dǎo)意義。

關(guān)鍵詞：巖石動力學(xué);圍壓;循環(huán)沖擊;砂巖;縱波波速;損傷

中圖分類號:TU452文獻標(biāo)志碼：A

基金項目：中國煤炭工業(yè)科技計劃項目(MTKJ2009-264)

收稿日期：2013-08-08修改稿收到日期：2013-12-04

Mechanical and ultrasonic analysis on damage of sandstone under cyclical impact loading with confining pressure

LIUShao-he1,XUJin-yu1,2,WANGPeng1,ZHANGGuo-xi3(1. Department of Airfield and Building Engineering, Air Force Engineering University, Xi’an 710038, Shaanxi, China；2. College of Mechanics and Civil Architecture, Northwest Polytechnic University, Xi’an 710072, Shaanxi, China；3. Base of Flight Training, Aviation University of Air Force, Changchun 130022, Jilin, China)

Abstract:Based on a φ100-mm-diameter split Hopkinson pressure bar with confining pressure apparatus uniaxial impact tests and cyclical impacting tests with confining pressure of 4 MPa and 20 MPa on sandstone specimens were carried out. Meanwhile, the ultrasonic longitudinal wave test was executed before and after each impact. The stress-strain properties of sandstone under cyclical impact were analyzed. The relationship between these properties and longitudinal wave velocity (Vp) in sandstone was analyzed and the impact damage was defined by the Vp. The yield-elasticity-ratio was defined to analyze the elastic-plastic characteristics of sandstones under confining pressure. The results show that, the stress-strain curves under confining pressure present a typical elastic-plastic feature. With the increase of cyclical impact times, yield stress and peak stress of sandstone sample reduce while yield strain and peak strain increase. The accumulative damage under lower confining pressure is apparently higher than that under higher confining pressure. The cyclical impact damage of sandstone shows an obvious confining pressure effect. The results have a guiding significance to the underground construction and protection.

Key words:rock dynamics; confining pressure; cyclical impact; sandstone; longitudinal wave velocity; damage

國民經(jīng)濟和國防建設(shè)的不斷發(fā)展促使巖石工程不斷向特殊地域、特殊環(huán)境延伸。隨著地下工程的深度不斷加大，地質(zhì)條件趨于復(fù)雜，深部巖體力學(xué)問題[1]已經(jīng)成為關(guān)系國家財產(chǎn)和人民生命安全的重大課題。深部巖體所處的“三高一擾動[1]”特殊地質(zhì)力學(xué)環(huán)境導(dǎo)致沖擊地壓、瓦斯爆炸等現(xiàn)象頻繁發(fā)生，由此引起的沖擊動荷載的不斷作用又使巖石發(fā)生損傷甚至破壞；實際上，大量工程巖體的宏觀失效破壞正是由地沖擊荷載循環(huán)作用引起的損傷累積造成的。深地下工程高圍壓環(huán)境下巖石在沖擊荷載循環(huán)作用下的累積損傷效應(yīng)[2]，已經(jīng)成為巖土力學(xué)與工程界的熱門前沿課題。

對此，國內(nèi)外進行了廣泛的研究，Schulze[3]應(yīng)用超聲波監(jiān)測技術(shù)研究了巖鹽變形過程中損傷發(fā)展演化規(guī)律。Meglis等[4]應(yīng)用超聲層析成像現(xiàn)場測試方法研究了隧道開挖誘發(fā)的圍巖損傷問題，得到了圍巖損傷程度和損傷分布規(guī)律。許金余等[5]使用帶圍壓裝置的SHPB系統(tǒng)對多種巖石循環(huán)沖擊能量特性進行了研究。林大能等[6]對大理巖試件在壓力試驗機上進行模擬沖擊加載，測試受沖擊后試件軸向超聲波波速，研究了巖石循環(huán)沖擊荷載作用下?lián)p傷的圍壓效應(yīng)。

本文對砂巖試樣進行了單軸沖擊和4 MPa、20 MPa圍壓情況下的循環(huán)沖擊試驗，并對每次沖擊前后的砂巖試樣進行超聲縱波檢測，分析了巖石在循環(huán)沖擊荷載作用下的屈服點和峰值點應(yīng)力、應(yīng)變的變化特性，采用縱波波速定義了砂巖試樣的沖擊損傷并分析了循環(huán)沖擊試驗中砂巖試樣縱波波速和應(yīng)力應(yīng)變之間的關(guān)系。研究結(jié)果對地下工程的建設(shè)和防護有一定的指導(dǎo)意義。

1砂巖基本物理力學(xué)特性

本文所用的巖石材料取自陜西太白山區(qū)秦嶺巖群。經(jīng)國土資源部西安礦產(chǎn)資源監(jiān)督檢測中心檢驗，巖樣礦物成分組成如表1所示。

表1　砂巖礦物成分組成

巖石是一種特殊的天然材料，內(nèi)部礦物成分和組織結(jié)構(gòu)呈現(xiàn)出顯著的多相性和不均勻性。為得到真實可靠的巖石材料的動態(tài)力學(xué)性能，需要采用大尺寸的試樣，進而也需要相應(yīng)的大直徑Hopkinson壓桿試驗裝置[7]。本文采用直徑為100 mm的SHPB試驗裝置。砂巖根據(jù)ISRM巖石力學(xué)試驗標(biāo)準(zhǔn)[8]制作成φ100×50 mm的圓柱體試樣。此外，為盡量排除巖樣差異性對試驗結(jié)果的影響，在試驗研究之前對其進行密度和超聲檢測篩選，所選用試樣的基本物理力學(xué)特性(平均值)如表2所示。

表2　砂巖基本物理力學(xué)特性

2圍壓條件下砂巖的循環(huán)沖擊試驗及超聲檢測

2.1試驗設(shè)備及試驗設(shè)計

原始SHPB試驗裝置主要由主體設(shè)備、能源系統(tǒng)、測試系統(tǒng)三大部分組成。本文與洛陽立特公司合作，對原始SHPB系統(tǒng)進行圍壓裝置改造，帶圍壓裝置的SHPB系統(tǒng)示意圖如圖1所示。

根據(jù)世界范圍內(nèi)的現(xiàn)場資料的統(tǒng)計[9]可知，絕大多數(shù)地下工程所處環(huán)境地應(yīng)力位于0-30 MPa區(qū)間內(nèi)。為此，本文圍壓設(shè)定為沒有圍壓(單軸狀態(tài)作為對照)、較低圍壓(C=4 MPa)和較高圍壓(C=20 MPa)三個等級。對于每一個圍壓等級采用三個巖石試樣進行試驗，試驗結(jié)果取平均值進行分析。

圖1　帶圍壓裝置的SHPB系統(tǒng)示意圖 Fig.1 SHPB test system with confining pressure device

經(jīng)單軸試驗可知，該批砂巖在無圍壓條件下，當(dāng)試驗入射應(yīng)力波的峰值強度僅為150 MPa(對應(yīng)的打擊桿速度為7.39 m/s)時，試樣已經(jīng)開裂。因此，為了研究圍壓條件下砂巖遭受循環(huán)沖擊荷載打擊的力學(xué)性能變化規(guī)律，確定試驗入射應(yīng)力波的最大峰值強度為345 MPa左右。

對于每一個砂巖試樣，以同一強度的沖擊荷載循環(huán)作用到使試件破壞為止。每一次沖擊前后均對砂巖試樣進行超聲檢測，記錄波形并讀取縱波波速值，每一次沖擊過程均采用超動態(tài)應(yīng)變儀對桿件內(nèi)的應(yīng)力波傳播狀態(tài)進行測試，直至試樣最終破壞。

2.2SHPB試驗基本原理

(1)

式中:A、As分別為壓桿、試件的截面積；ls為試件的初始厚度；C0為壓桿的彈性縱波波速；P1、P2分別為作用在試件兩端的力；E是壓桿的彈性模量；A是壓桿的橫截面面積。

3試驗結(jié)果分析

3.1圍壓條件下砂巖試樣的彈塑性特征

本文由三波法得到單軸狀態(tài)和圍壓狀態(tài)的應(yīng)力應(yīng)變曲線如圖2所示。

圖2　單軸狀態(tài)和圍壓(4 MPa) 狀態(tài)下的動態(tài)應(yīng)力-應(yīng)變曲線 Fig.2 Dynamic stress-strain curves under uniaxial compression or confining pressure

在無圍壓條件下，砂巖試件的破壞屬于脆性破壞，而在圍壓作用下，應(yīng)力-應(yīng)變曲線呈現(xiàn)典型的彈塑性特征，在達到屈服應(yīng)力后，應(yīng)力增長不多，而應(yīng)變持續(xù)增加，即塑性應(yīng)變增大，出現(xiàn)明顯的屈服平臺，巖樣的破壞屬于延性破壞，具有明顯的脆性-延性轉(zhuǎn)化特征[12]。塑性變形的本質(zhì)是材料整體發(fā)生相對滑移或剪切，因而不發(fā)生剪切失穩(wěn)是塑性變形的前提[13]?？芍?，圍壓的加載作用對阻止試件發(fā)生剪切失穩(wěn)的作用相當(dāng)明顯。

圍壓作用下砂巖的應(yīng)力應(yīng)變曲線可以劃分為三個階段，即彈性階段、屈服階段和峰后階段，不同圍壓條件下屈服階段長度(以屈服階段應(yīng)變增量衡量)不同，其與彈性階段長度(以彈性階段應(yīng)變增量衡量)的比值也有所不同。本文以屈服階段應(yīng)變增量與彈性階段應(yīng)變增量之比定義屈服-彈性比R，用以表征砂巖應(yīng)力應(yīng)變塑性特征的顯著程度，其計算表達式如下：

(2)

式中：R為屈服-彈性比，εp、εf分別為峰值應(yīng)變、屈服應(yīng)變。砂巖試樣循環(huán)沖擊作用下，屈服-彈性比隨沖擊作用次數(shù)的變化規(guī)律如圖3所示。

圖3　沖擊荷載循環(huán)作用次數(shù)與屈服-彈性比的關(guān)系 Fig.3 Relationship between circular action times and R

由圖3可知，與較低圍壓條件(C=4 MPa)下相比，較高圍壓條件(C=20 MPa)下砂巖試樣的屈服-彈性比明顯降低，砂巖應(yīng)力應(yīng)變曲線屈服階段變短。隨著沖擊作用次數(shù)增加，兩種圍壓條件下砂巖試樣在最終破壞前屈服-彈性比均增大，砂巖試樣的塑性特征越發(fā)明顯。最后一次沖擊，砂巖試樣發(fā)生斷裂，試樣屈服-彈性比的變化規(guī)律與破壞前不同。屈服-彈性比的變化規(guī)律反映出試樣隨沖擊作用次數(shù)的增加內(nèi)部結(jié)構(gòu)發(fā)生變化，動載作用下塑性特征變得明顯。

3.2圍壓條件下砂巖試樣的應(yīng)力應(yīng)變特征

表3所列為不同圍壓狀態(tài)下砂巖試樣各次沖擊試驗的屈服應(yīng)力應(yīng)變、峰值應(yīng)力應(yīng)變值，同時列出的還有各次沖擊前后超聲檢測獲取的縱波波速值，用以反映各次沖擊后砂巖試樣的損傷狀態(tài)。為更直觀的描述砂巖試樣屈服應(yīng)力應(yīng)變、峰值應(yīng)力應(yīng)變隨沖擊荷載循環(huán)作用次數(shù)的變化規(guī)律，繪制圖4和圖5。

圖4　沖擊荷載循環(huán)作用次數(shù)與屈服應(yīng)力和屈服應(yīng)變關(guān)系 Fig.4 Relationship between circular action times and yield stress, yield strain

如圖4所示，隨著沖擊荷載循環(huán)作用次數(shù)的增加，試樣的屈服應(yīng)力降低，屈服應(yīng)變增加。在沖擊荷載作用下，砂巖試樣進入塑性變形階段的屈服平臺后，試樣內(nèi)部不斷產(chǎn)生塑性變形，從而導(dǎo)致材料性能的劣化，宏觀上表現(xiàn)為在后續(xù)循環(huán)加載時試樣的屈服應(yīng)力降低。同時，高圍壓條件下砂巖試樣的屈服應(yīng)力值較低圍壓條件高而屈服應(yīng)變值低。

表3　不同圍壓狀態(tài)下砂巖試樣循環(huán)沖擊試驗結(jié)果

注：屈服應(yīng)力應(yīng)變值指的是圍壓狀態(tài)下砂巖應(yīng)力應(yīng)變曲線由彈性階段進入屈服階段的拐點處的應(yīng)力應(yīng)變值。

圖5　沖擊荷載循環(huán)作用次數(shù)與峰值應(yīng)力和峰值應(yīng)變關(guān)系 Fig.5 Relationship betweencircular action times and peak stress, peak strain

如圖5所示，隨著沖擊荷載循環(huán)作用次數(shù)的增加，試樣的峰值應(yīng)力應(yīng)變具有與屈服應(yīng)力應(yīng)變相似的變化規(guī)律。圍壓的增加可以極大的提高砂巖的抗壓強度及變形性能。且隨著沖擊荷載循環(huán)作用次數(shù)的增加，圍壓的改善作用越明顯。巖石內(nèi)存在各種方向的裂隙，無圍壓時所有傾角大于內(nèi)摩擦角的裂隙承載能力為零；隨著圍壓的增大，陡傾角裂隙的承載能力迅速增大，延緩了滑移破壞的產(chǎn)生，砂巖承載力提高。

3.3圍壓條件下砂巖試樣循環(huán)沖擊特性超聲分析

砂巖試樣屈服應(yīng)力應(yīng)變、峰值應(yīng)力應(yīng)變隨沖擊荷載循環(huán)作用次數(shù)而變化。沖擊荷載循環(huán)作用次數(shù)不具備實質(zhì)性的物理力學(xué)含義，應(yīng)力應(yīng)變值的變化究其原因在于沖擊荷載循環(huán)作用引起砂巖試樣不同程度的損傷，隨著沖擊荷載作用次數(shù)的增加，砂巖損傷不斷積累直至破壞。超聲波經(jīng)過巖樣將會攜帶大量的巖樣信息[14]。砂巖宏觀的超聲特性能夠較好地反映其內(nèi)部的損傷狀態(tài)。趙洪寶等[15]認為損傷材料發(fā)生損傷后應(yīng)力-應(yīng)變關(guān)系應(yīng)滿足：

σ=E(1-D)εe

(3)

并基于縱波波速與介質(zhì)彈性系數(shù)的關(guān)系定義巖石

損傷因子為:

(4)

式中:Vp、Vf分別為損傷前、后砂巖中的縱波波速。

基于式(4)，取Vp為砂巖試樣每次沖擊后的縱波波速，Vf為未遭受沖擊砂巖試樣的縱波波速，可以構(gòu)建砂巖試樣累積損傷隨沖擊荷載循環(huán)作用次數(shù)的變化關(guān)系，如圖6所示?？梢娒看螞_擊對砂巖試樣的損傷并不均等。與較高圍壓(C=20 MPa)下相比，較低圍壓(C=4 MPa)狀態(tài)下砂巖累計損傷更高。需要說明的是，試樣最后一次沖擊作用破壞后，無法測得試樣的準(zhǔn)確波速值，而將其縱波波速值視為0或者將其損傷值視為1均欠妥，故本文未對其進行討論。

圖6　砂巖累積損傷隨沖擊荷載作用次數(shù)的關(guān)系 Fig.6 Relationship between accumulative damage of sandstone and circular action times

觀察表3 可知，作為損傷的表征量，縱波波速與屈服應(yīng)力應(yīng)變、峰值應(yīng)力應(yīng)變隨作用次數(shù)的增加具有相似的變化規(guī)律。為研究不同圍壓條件下，砂巖試樣力學(xué)特性與各次沖擊作用前的超聲縱波波速間的相關(guān)關(guān)系，本文分析了屈服應(yīng)力應(yīng)變、峰值應(yīng)力應(yīng)變與各次沖擊作用前縱波波速的相關(guān)關(guān)系，如圖7、圖8所示。

圖7　砂巖屈服應(yīng)力應(yīng)變隨縱波波速的關(guān)系 Fig.7 Yield stress and strain versus longitudinal wave velocity

圖8　砂巖峰值應(yīng)力應(yīng)變隨縱波波速的關(guān)系 Fig.8 Peak stress and strain versus longitudinal wave velocity

分析圖7、圖8 可知，砂巖試樣縱波波速與屈服應(yīng)力、屈服應(yīng)變之間具有較強的線性相關(guān)關(guān)系。其中，縱波波速與屈服應(yīng)力為正相關(guān)關(guān)系，與屈服應(yīng)變?yōu)樨撓嚓P(guān)關(guān)系。砂巖試樣峰值應(yīng)力、峰值應(yīng)變而講，隨縱波波速的變化規(guī)律總體上呈現(xiàn)單調(diào)性，但并不是簡單的線性關(guān)系。

4結(jié)論

本文以秦嶺太白山區(qū)砂巖為試驗對象進行了單軸沖擊和4 MPa、20 MPa圍壓條件下的循環(huán)沖擊試驗，并對每次沖擊前后的砂巖試樣進行超聲縱波檢測，分析結(jié)果發(fā)現(xiàn)：

(1)在無圍壓條件下，砂巖試件的破壞屬于脆性破壞；而在圍壓作用下，應(yīng)力-應(yīng)變曲線出現(xiàn)較明顯的屈服平臺，呈現(xiàn)出典型的彈塑性特征。隨著沖擊作用次數(shù)增加，砂巖試樣在最終破壞前屈服-彈性比均增大，砂巖試樣的塑性特征越發(fā)明顯。

(2)隨著沖擊荷載循環(huán)作用次數(shù)的增加，砂巖試樣屈服應(yīng)力、峰值應(yīng)力降低，屈服應(yīng)變、峰值應(yīng)變增加。較低圍壓(C=4 MPa)狀態(tài)下應(yīng)力和應(yīng)變隨沖擊作用次數(shù)變化更明顯，較高圍壓(C=20 MPa)狀態(tài)下則較小。

(3)隨著循環(huán)沖擊作用次數(shù)的增加，砂巖試樣的應(yīng)力、應(yīng)變特征與縱波波速之間存在良好的相關(guān)關(guān)系。基于砂巖試樣超聲波速定義了沖擊累積損傷。每次沖擊對砂巖試樣的損傷并不均等，較低圍壓狀態(tài)下砂巖試樣累計損傷明顯高于較高圍壓下砂巖的累積損傷，砂巖循環(huán)沖擊損傷具有明顯的圍壓效應(yīng)。

參考文獻

[1]何滿潮, 錢七虎. 深部巖體力學(xué)研究進展[C]. 中國巖石力學(xué)與工程學(xué)會編. 第九屆全國巖石力學(xué)與工程學(xué)術(shù)大會論文集. 北京: 科學(xué)出版社, 2006: 49-62.

[2]許金余, 范建設(shè), 呂曉聰. 圍壓條件下巖石的動態(tài)力學(xué)特性[M]. 西安: 西北工業(yè)大學(xué)出版社, 2012. 01.

[3]Sehulze O, Popp T, Kem H. Development of damage and Permeability in deforming rock salt[J]. Engineering Geology, 2001, 61(2):163-180.

[4]Meglis I L, Chow T M, Martin C D. Assessing in situ micro crack damage using ultrasonic velocity tomography[J]. International Journal of Rock Mechanics & Mining Sciences, 2005, 42(l): 25-34.

[5]許金余, 呂曉聰, 張軍, 等. 圍壓條件下巖石循環(huán)沖擊損傷的能量特性研究[J]. 巖石力學(xué)與工程學(xué)報, 2010, 29(S2): 4159-4165.

XU Jin-yu, LU Xiao-cong, ZHANG Jun, et al. Research on energy properties of rock cyclical impact damage under confining pressure[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(S2): 4159-4165.

[6]林大能, 陳壽如, 劉優(yōu)平. 巖石循環(huán)沖擊損傷演化圍壓效應(yīng)的實驗研究[J]. 礦冶工程, 2005, 25(4): 12-15.

LIN Da-neng, CHEN Shou-ru, LIU You-ping. Experimental study on confining pressure effect of rock damage under cyclical impact loading[J]. Mining and Metaliurgical Engineering, 2005, 25(4): 12-15.

[7]陳德興, 胡時勝, 張守保, 等. 大尺寸Hopkinson壓桿及其應(yīng)用[J].實驗力學(xué), 2005, 20(3): 398-402.

CHEN De-xing, HU Shi-sheng, ZHANG Shou-bao. Larger dimension hopkinson pressure bar and it’s application[J]. Journal of Experimental Mechanics, 2005, 20(3): 398-402.

[8]International Society for Rock Mechanics. Suggested methods for determining the uniaxial compressive strength and deformability of rock materials[S].[S. l.]: [s. n.], 1979.

[9]Brown E T，Hoek E. Trends in relationships between measured rock in situ stress and depth[J]. Int J Rock Mech Min Sci Geornech Abstr,1978,15：21l-215.

[10]Ravichandran G, Subhash G. Critical appraisal of limiting strain rates for compression testing ceramics in a split Hopkinson pressure bar[J]. Journal of American Ceramic Society, 1994,77(1): 263-267.

[11]王禮立. 應(yīng)力波基礎(chǔ)[M].北京: 國防工業(yè)出版社, 2005: 52-60.

[12]Wang Peng, Xu Jin-yu, Liu Shi. Staged moduli: A quantitative method to analyze the complete compressive sress-strain response for thermally damaged rock[J]. Rock Mech Rock Eng, 2004.

[13]薛志剛, 胡時勝. 水泥砂漿在圍壓下的動態(tài)力學(xué)性能[J]. 工程力學(xué), 2008, 25(12): 184-188.

XUE Zhi-gang, HU Shi-sheng. Dynamic behavior of cement mortar under confinement[J]. Engineering Mechanics, 2008, 25(12): 184-188.

[14]王鵬, 許金余, 劉石, 等. 熱損傷砂巖力學(xué)與超聲時頻特性研究[J]. 巖石力學(xué)與工程學(xué)報, 2014, 33(9): 1897-1904.

WANG Peng, XU Jin-yu, LIU Shi, et al. Mechanical properties and ultrasonic time-frequency characteristics of thermally damaged sandstone[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(9): 1897-1904.

[15]趙洪寶, 尹光志, 諶倫建. 溫度對砂巖損傷影響試驗研究[J]. 巖石力學(xué)與工程學(xué)報，2009，28(s1): 2784-2788.

ZHAO Hong-bao, YIN Guang-zhi, CHEN Lun-jian. Experimental study on effect of temperature on sandstone damage[J]. Chinese Journal of Rock Mechanics and Engineering, 2009, 28(s1): 2784-2788.

第一作者趙麗娟女，博士，教授，博士生導(dǎo)師,1964年4月生

通信作者田震男，博士生，1987年1月生