基于繞射信息提取技術(shù)的斷裂識(shí)別方法及應(yīng)用

肖 曦 周 鵬 張益明 霍守東 王志紅 陳建華

(①中海油研究總院有限責(zé)任公司，北京 100028；②中國(guó)科學(xué)院地質(zhì)與地球物理研究所，北京 100029；③中國(guó)海洋石油國(guó)際有限公司，北京 100028)

0 引言

鄂爾多斯盆地是中國(guó)第二大沉積盆地，油氣資源豐富[1]。該盆地上古生界發(fā)育致密砂巖氣藏，已發(fā)現(xiàn)地質(zhì)儲(chǔ)量超過(guò)5×1012m3，勘探潛力巨大[2]。A區(qū)塊位于鄂爾多斯盆地東緣，現(xiàn)有探井均有氣層發(fā)現(xiàn)，氣層鉆遇率高，但不同井的氣層測(cè)試產(chǎn)量差異很大，其中一半以上的井測(cè)試產(chǎn)量較低，不具有商業(yè)開(kāi)發(fā)價(jià)值。分析井孔的成藏特征表明，距斷層、裂縫近的井往往氣層厚度更大，即斷層、裂縫影響油氣富集，具有后期調(diào)整、再分配氣層分布的作用。因此，斷裂體系研究對(duì)致密氣儲(chǔ)層的井位優(yōu)選具有重要意義。

前人廣泛研究了刻畫斷裂體系的方法。常用的基于疊后屬性的方法包括相干體屬性[3-4]、邊緣檢測(cè)[5]、曲率屬性[6]、方差體[7]等斷裂識(shí)別方法，通過(guò)研究道間不連續(xù)變化以及構(gòu)造變形預(yù)測(cè)斷層和構(gòu)造突變區(qū)，從而預(yù)測(cè)斷裂發(fā)育帶。此外，近年來(lái)人們深入研究了基于機(jī)器學(xué)習(xí)的斷裂識(shí)別技術(shù)。機(jī)器學(xué)習(xí)是實(shí)現(xiàn)人工智能的重要方法之一，主要目的是從數(shù)據(jù)中直接獲取知識(shí)。在過(guò)去20年中，發(fā)展了數(shù)十種機(jī)器學(xué)習(xí)算法，如遺傳算法、深度學(xué)習(xí)、人工神經(jīng)網(wǎng)絡(luò)等[8]。根據(jù)地震資料，人們應(yīng)用深度學(xué)習(xí)方法預(yù)測(cè)斷層和裂縫。Cai等[9]利用卷積神經(jīng)網(wǎng)絡(luò)并采用多種疊后地震屬性解釋斷層。段艷廷等[10]基于三維半密度卷積神經(jīng)網(wǎng)絡(luò)檢測(cè)斷裂，該法具有一定的遷移能力和較好的泛化能力。Di等[11]提出一種基于多層感知器(MLP)技術(shù)的地震斷層自動(dòng)檢測(cè)方法。

A區(qū)塊地表為黃土塬區(qū)，區(qū)域黃土厚度變化大，復(fù)雜地表的強(qiáng)烈吸收造成地震波能量急劇衰減，地震資料分辨率低[12]。該區(qū)致密氣儲(chǔ)層非均質(zhì)性強(qiáng)、厚度小，斷層、裂縫體系通常斷距小，小尺度裂縫(斷距較小的裂縫)豐富，常規(guī)數(shù)據(jù)難以滿足薄儲(chǔ)層中斷層、裂縫體系識(shí)別的需求。此外，基于疊后屬性的方法主要以地震反射數(shù)據(jù)為主[3-7]，分辨率受菲涅耳帶半徑和鏡面反射波長(zhǎng)約束，對(duì)小尺度地質(zhì)異常體探測(cè)能力較弱[13]。鑒于深度學(xué)習(xí)的大數(shù)據(jù)特性，業(yè)界主要將人工合成地震記錄及斷層標(biāo)簽作為訓(xùn)練數(shù)據(jù)學(xué)習(xí)斷層特征，最終識(shí)別效果受合成地震記錄的精度、斷層標(biāo)簽的準(zhǔn)確性以及網(wǎng)絡(luò)參數(shù)的合理性等多因素影響，目前大部分方法仍停留在實(shí)驗(yàn)室階段，在油田實(shí)際工作中應(yīng)用較少。

通過(guò)繞射信息識(shí)別小尺度地質(zhì)體是近年的研究熱點(diǎn)[14-17]。地震繞射波攜帶了大量的與地下突變點(diǎn)(斷層、溶洞、斷裂系統(tǒng)等)有關(guān)的信息。由于繞射波信息能量弱，常淹沒(méi)在能量較強(qiáng)的反射波中，因此難以準(zhǔn)確識(shí)別。通過(guò)分離繞射波與反射波以探測(cè)小尺度斷裂分布。為此，人們提出了不同的繞射波分離成像方法，包括疊加、濾波、Radon變換—反Radon變換、聚焦—反聚焦、平面波分解等方法[14]。陳明政等[15]應(yīng)用繞射波成像技術(shù)，對(duì)塔河油田的碳酸鹽巖小尺度縫洞儲(chǔ)集體(具弱反射特征)的識(shí)別效果較好。李繼光[16]利用繞射波成像技術(shù)提高了河道砂體識(shí)別能力。然而目前鮮有利用繞射波信息識(shí)別致密氣儲(chǔ)層中斷裂體系的研究實(shí)例。

鑒于該區(qū)致密氣儲(chǔ)層厚度小、斷裂體系斷距小、小尺度裂縫豐富，常規(guī)數(shù)據(jù)及常規(guī)基于地震反射信息的斷裂識(shí)別方法無(wú)法高精度預(yù)測(cè)儲(chǔ)層中斷裂體系分布，筆者首先利用基于L0范數(shù)稀疏反演的壓縮感知地震資料處理方法提高數(shù)據(jù)分辨率，為后續(xù)薄儲(chǔ)層中斷裂體系識(shí)別創(chuàng)造良好的資料條件。在此基礎(chǔ)上探索基于中值濾波技術(shù)的繞射信息提取方法，在保幅偏移地震數(shù)據(jù)體上提取繞射信息，得到清晰的斷點(diǎn)和斷層信息；根據(jù)繞射波地震數(shù)據(jù)振幅及相位在空間上的不連續(xù)性，通過(guò)螞蟻體追蹤刻畫斷層和裂縫分布。結(jié)果表明：處理后的地震數(shù)據(jù)分辨率得到有效提高，并且斷點(diǎn)清晰，更易識(shí)別斷層；與常規(guī)基于疊后屬性的方法相比，基于繞射信息提取技術(shù)的斷裂識(shí)別方法能夠識(shí)別更小尺度的斷裂分布區(qū)。上述結(jié)果為探索該區(qū)氣層發(fā)育的主控因素奠定了基礎(chǔ)，為促進(jìn)天然氣的增儲(chǔ)上產(chǎn)提供指導(dǎo)。

1 斷裂體系識(shí)別關(guān)鍵技術(shù)

1.1 基于壓縮感知的寬頻數(shù)據(jù)處理

壓縮感知(Compressed Sensing, CS)是應(yīng)用數(shù)學(xué)和信號(hào)處理領(lǐng)域近十年來(lái)的最新研究成果，也稱為基于稀疏性的信號(hào)采樣理論，在遠(yuǎn)小于Nyquist采樣率的條件下，利用隨機(jī)采樣獲取信號(hào)的離散樣本，把高維空間的信號(hào)通過(guò)觀測(cè)矩陣投影到低維空間，進(jìn)而通過(guò)非線性重建算法完備重建信號(hào)[18]。CS需要滿足三個(gè)條件：①信號(hào)是稀疏的或可壓縮的；②觀測(cè)矩陣是隨機(jī)的，和信號(hào)本身互不相關(guān)；③通過(guò)一定的重建算法完成信號(hào)重建。根據(jù)地震褶積模型，地震數(shù)據(jù)s可以描述為地震子波w與地下反射系數(shù)r的褶積和隨機(jī)噪聲n之和

s=w*r+n

(1)

上式在頻率域可以表示為

S=WR+N

(2)

式中S、W、R、N分別為地震記錄、地震子波、地下反射系數(shù)和隨機(jī)噪聲的Fourier變換。

假設(shè)地震反射系數(shù)是稀疏的，通過(guò)地震數(shù)據(jù)恢復(fù)寬帶反射系數(shù)滿足CS理論[19]。通過(guò)重建算法尋求式(2)的最稀疏解，最直接的方法就是在L0范數(shù)下求解最優(yōu)化問(wèn)題，使非零元素的個(gè)數(shù)最少。L0范數(shù)求解是個(gè)非確定性多項(xiàng)式難題(NP難題，僅有局部最優(yōu)解)，通過(guò)將匹配追蹤(MP)方法得到的初始值作為輸入，可以得到較好的結(jié)果。MP是一種信號(hào)稀疏分解技術(shù)，通過(guò)將地震信號(hào)分解為與局部結(jié)構(gòu)特征最相似的基元函數(shù)的線性組合，在此基元函數(shù)集合上稀疏分解，能夠反演地震反射系數(shù)[20]。該過(guò)程使基元函數(shù)的線性組合逼近地震數(shù)據(jù)，但當(dāng)反射系數(shù)個(gè)數(shù)偏多時(shí)，算法收斂性不高。將MP得到的反射系數(shù)作為L(zhǎng)0范數(shù)稀疏反演的初始值，可有效降低稀疏反演的迭代次數(shù)，使得到的反射系數(shù)更好地反映地層結(jié)構(gòu)，斷點(diǎn)清晰，更易識(shí)別斷層。

本次研究利用MP獲取初始反射系數(shù)，采用基于L0范數(shù)的稀疏反演方法探索提高分辨率的方法。其目標(biāo)方程為

(3)

式中：γ為正則化因子；C為非零反射系數(shù)個(gè)數(shù)。

通過(guò)L0范數(shù)約束優(yōu)化問(wèn)題，求取稀疏表達(dá)的反射系數(shù)，在提高分辨率的同時(shí)，識(shí)別薄層界面、透鏡體邊界和地層尖滅位置等薄層結(jié)構(gòu)。在L0范數(shù)稀疏反演反射系數(shù)基礎(chǔ)上，進(jìn)一步通過(guò)Morlet小波拓頻替換[21]拓寬頻帶。

基于CS的提高分辨率的具體方法[22]為：①利用MP方法獲取初始反射稀疏模型；②以MP算法得到的結(jié)果作為初始反射系數(shù)，進(jìn)行基于CS算法的疊后稀疏反演；③對(duì)已知井進(jìn)行合成地震記錄標(biāo)定，獲取地震數(shù)據(jù)的初始子波；④在提取的子波約束下，通過(guò)Morlet小波連續(xù)變換，嚴(yán)格控制變換子波的頻率域峰值和相位與原子波完全一致，進(jìn)而通過(guò)連續(xù)小波逆變換提高地震數(shù)據(jù)頻帶寬度。

1.2 基于繞射信息提取的斷裂識(shí)別

與常規(guī)反射信息相比，地震數(shù)據(jù)中的繞射波與小尺度地質(zhì)體有關(guān)，且具有較高的分辨率，是精細(xì)描述小尺度地質(zhì)異常體的重要信息。前人的研究表明，繞射成像能夠檢測(cè)小尺度斷裂信息，最小尺度達(dá)5～10m[23]。疊后數(shù)據(jù)中依然存在有效繞射信息，能夠很好地反映特殊構(gòu)造(溶洞、斷層、斷裂系統(tǒng)等)的空間位置及規(guī)模等重要信息[24]。

反射信息分布在局部光滑波前上，振幅變化較緩慢。由繞射響應(yīng)函數(shù)衰減特征可知，繞射振幅、波形沿繞射走時(shí)曲線變化很快。經(jīng)偏移處理后，繞射能量殘存于端點(diǎn)成像位置的一定區(qū)域內(nèi)。

反射成像結(jié)果的波形、振幅穩(wěn)定，空間相關(guān)性相對(duì)較好。中值濾波技術(shù)已被用于預(yù)測(cè)反射信息(圖1)：當(dāng)中值濾波窗口完全位于反射同相軸時(shí)，可完全預(yù)測(cè)反射信息；當(dāng)中值濾波窗口一端位于反射軸、另一端位于端點(diǎn)成像干涉帶內(nèi)時(shí)，在端點(diǎn)附近繞射振幅和相位的變化均較快，通常為非線性變化。因此，經(jīng)中值濾波后，難以預(yù)測(cè)繞射信息，即越接近端點(diǎn)，不能預(yù)測(cè)的繞射殘留能量越多。

圖1 中值濾波與反射信息預(yù)測(cè)示意圖[24]

由上述分析可知，通過(guò)估計(jì)地震數(shù)據(jù)的傾角，沿著傾向進(jìn)行中值濾波，可以預(yù)測(cè)反射信息。由于繞射產(chǎn)生的振幅、相位相對(duì)不穩(wěn)定區(qū)域接近斷點(diǎn)，存在殘留能量，因此能夠分離非均質(zhì)地質(zhì)體產(chǎn)生的繞射波場(chǎng)，但分離的波場(chǎng)也存在部分噪聲。由于由斷裂系統(tǒng)產(chǎn)生的邊緣繞射通常在空間上具有一定的相關(guān)性，將Radon變換應(yīng)用于繞射成像結(jié)果，充分利用繞射波場(chǎng)的空間相關(guān)性進(jìn)行繞射特征增強(qiáng)處理，可以得到直觀、可靠的繞射波分析結(jié)果[13]，該過(guò)程類似于基于Radon變換的數(shù)據(jù)插值處理。具體方法為：①確定時(shí)窗，選擇時(shí)窗內(nèi)地震數(shù)據(jù)；②對(duì)時(shí)窗內(nèi)地震數(shù)據(jù)Radon變換；③在Radon域選擇合適閾值，保留在Radon域較為聚集的能量點(diǎn)(即繞射波具有一定的空間相關(guān)性的點(diǎn))，進(jìn)而確定主要繞射點(diǎn)的主要延伸方向；④通過(guò)Radon反變換增強(qiáng)繞射特征。

1)榆林市礦區(qū)樣本水化學(xué)類型主要為HCO3-Ca以及HCO3-Ca·Mg，其中超過(guò)《地下水質(zhì)量標(biāo)準(zhǔn)》的第Ⅲ類標(biāo)準(zhǔn)離子有8個(gè)，分別為TDS、TH、Cl-、SO2-4、NO-2、NO-3、F-以及CODCr。其中NO-3的超標(biāo)水樣數(shù)量最多，為15個(gè)，超標(biāo)率為35.7%，最大的超標(biāo)倍數(shù)接近5倍，其次為NO-2、SO2-4、TDS和總硬度。SO2-4超標(biāo)數(shù)為2個(gè)，超標(biāo)率為4.8%；NO-2超標(biāo)數(shù)為3個(gè)，超標(biāo)率為7.2%；其余的指標(biāo)超標(biāo)水量均為1個(gè)。

2 模型試算

為了說(shuō)明繞射波對(duì)斷裂的識(shí)別能力，建立不同尺度裂縫的全波場(chǎng)二維模型(圖2)進(jìn)行試算。通過(guò)Kirchhoff算法合成零炮檢距地震記錄，并通過(guò)前述方法提取反射波和繞射波剖面。圖2為小斷裂正演模擬結(jié)果。由圖可見(jiàn)：在零炮檢距地震記錄(圖2a)上，垂直于斷層方向存在繞射波；當(dāng)斷層斷距小于λ/16(λ為波長(zhǎng))時(shí)，反射波分辨率已經(jīng)明顯不足(圖2b)；繞射波對(duì)斷點(diǎn)有較好的響應(yīng)(圖2c)。假設(shè)地層速度為4500m/s，地震波主頻達(dá)到50Hz，則λ=90m，那么在較理想的情況下，繞射波分辨率(λ/16)可以達(dá)到5.625m。

圖2 小斷裂正演模擬結(jié)果

3 應(yīng)用效果分析

A區(qū)塊位于鄂爾多斯盆地東緣，主要受加里東期、印支期、燕山期和喜馬拉雅期等古構(gòu)造應(yīng)力場(chǎng)影響[25]?；趹?yīng)力場(chǎng)的區(qū)域斷層、裂縫分析表明，該區(qū)主要裂縫方向?yàn)榻黃N向、NW—SE向和NE—SW向。該區(qū)及周邊地層自下而上為下奧陶統(tǒng)、中石炭統(tǒng)本溪組，下二疊統(tǒng)太原組、山西組、下石盒子組，上二疊統(tǒng)上石盒子組。主力產(chǎn)氣層段為本溪組、太原組。應(yīng)用上述斷裂體系識(shí)別關(guān)鍵技術(shù)研究A區(qū)塊儲(chǔ)層中斷裂體系，探索該區(qū)氣層發(fā)育的主控因素，以促進(jìn)天然氣的增儲(chǔ)上產(chǎn)。

圖3為原始數(shù)據(jù)和寬頻數(shù)據(jù)的疊加剖面及頻譜。由圖可見(jiàn)：原始數(shù)據(jù)主頻約為30Hz，最高頻率約為70Hz(圖3a下)，原始資料斷點(diǎn)不清晰(圖3a上黑色橢圓和黑色箭頭指示區(qū)域)；寬頻數(shù)據(jù)(圖3b上)主頻約為60Hz，最高頻率約為120Hz(圖3b下)，分辨能力得到提高，斷點(diǎn)更清晰。

將寬頻處理數(shù)據(jù)濾波，并與原始數(shù)據(jù)對(duì)比，測(cè)試有效頻帶范圍。圖4為原始數(shù)據(jù)與寬頻數(shù)據(jù)。由圖可見(jiàn)：①原始數(shù)據(jù)剖面中存在兩個(gè)斷裂區(qū)域(黑色橢圓區(qū)域)，斷點(diǎn)清晰可辨(圖4a)；②將寬頻數(shù)據(jù)濾除30Hz以下成分后，斷裂分布帶(黑色橢圓區(qū)域)與圖4a形態(tài)接近，但分辨率明顯提高(圖4b)；③將寬頻數(shù)據(jù)分別濾除60Hz以下成分(圖4c)和90Hz以下成分(圖4d)后，斷裂形態(tài)與圖4a較接近，斷點(diǎn)位置幾乎相同；④將寬頻數(shù)據(jù)分別濾除120Hz以下成分(圖4e)和150Hz以下成分(圖4f)后，斷點(diǎn)位置(黑色橢圓區(qū)域)與圖4a差異較大，并且出現(xiàn)“斷面波”(黑色箭頭指示區(qū)域)，推測(cè)是高頻噪聲。因此有效頻帶范圍與圖3b下基本一致，最高不超過(guò)120Hz，有效頻帶范圍為100Hz以內(nèi)。圖3和圖4表明，寬頻數(shù)據(jù)濾波處理有效提高了數(shù)據(jù)的時(shí)間和空間分辨率，為后續(xù)致密砂巖斷裂體系識(shí)別研究奠定了堅(jiān)實(shí)的基礎(chǔ)。

圖3 原始數(shù)據(jù)(a)和寬頻數(shù)據(jù)(b)的疊加剖面(上)及頻譜(下)白色橢圓區(qū)域指示煤層和砂巖互層

圖4 原始數(shù)據(jù)與寬頻數(shù)據(jù)剖面

通過(guò)繞射波提取方法，從反射地震數(shù)據(jù)(圖5a)中提取繞射波屬性(圖5b)，有效壓制了反射界面，繞射波發(fā)育帶能夠清晰地反映斷點(diǎn)信息。

圖5 基于CS的寬頻地震剖面(a)與繞射波屬性剖面(b)

根據(jù)繞射波振幅及相位在空間上的不連續(xù)性，應(yīng)用螞蟻體追蹤技術(shù)刻畫小尺度斷裂體系。螞蟻體追蹤算法通過(guò)在地震數(shù)據(jù)中散播大量的“螞蟻”，滿足預(yù)設(shè)斷裂條件的“螞蟻”將釋放某種信號(hào)，召集其他區(qū)域的螞蟻集中在斷裂處進(jìn)行追蹤，直到完成對(duì)斷裂的追蹤和識(shí)別[26]。對(duì)繞射波屬性運(yùn)用該方法，能夠突顯小尺度斷裂的分布特征。

為驗(yàn)證斷裂體系識(shí)別關(guān)鍵技術(shù)在小尺度斷裂識(shí)別中的應(yīng)用效果，將基于原始數(shù)據(jù)的相干屬性與基于寬頻數(shù)據(jù)的繞射波屬性進(jìn)行對(duì)比。將原始數(shù)據(jù)和相干屬性疊合顯示(圖6a)，通過(guò)相干屬性可識(shí)別尺度為幾十米到幾百米的大斷裂，可見(jiàn)大斷裂主要發(fā)育在下石盒子組和奧陶系，太原組、本溪組并不發(fā)育。分析鉆井?dāng)?shù)據(jù)可知，該區(qū)氣層主要發(fā)育在太原組和本溪組，因此大斷裂分布與氣層分布相關(guān)性不大，即大斷裂不是氣層富集的關(guān)鍵因素。進(jìn)一步將基于CS處理的寬頻數(shù)據(jù)與基于繞射波數(shù)據(jù)提取的螞蟻體屬性疊合顯示(圖6b)，可見(jiàn)繞射波數(shù)據(jù)刻畫的中小尺度斷裂在下石盒子組、山西組、太原組、本溪組和奧陶系均有分布，其中太原組和本溪組的斷裂細(xì)節(jié)(藍(lán)色橢圓區(qū)域)比圖6a的大斷裂更豐富，說(shuō)明繞射波數(shù)據(jù)在消除大尺度斷裂的波場(chǎng)特征后，能夠突出中小尺度斷裂系統(tǒng)特征。繞射波數(shù)據(jù)刻畫的中小尺度斷裂整體分布規(guī)律與氣層分布規(guī)律一致，表明中小尺度斷裂有利于氣層成藏。

圖6 原始數(shù)據(jù)與相干屬性疊合剖面(a)、寬頻數(shù)據(jù)與基于繞射波數(shù)據(jù)提取的螞蟻體屬性疊合剖面(b)

圖7、圖8分別為本溪組儲(chǔ)層頂面相干平面圖、基于繞射波提取的本溪組儲(chǔ)層頂面螞蟻體屬性平面圖。由圖可見(jiàn)：相干屬性能夠識(shí)別大尺度不連續(xù)性，較清晰地刻畫了斷層及河道邊界(圖7)；基于繞射信息提取的斷裂屬性直觀地識(shí)別了更小尺度的斷裂發(fā)育帶，斷裂的主要發(fā)育方向分別為NNE和NNW向(圖8)。綜上所述，推斷燕山期鄂爾多斯盆地構(gòu)造應(yīng)力場(chǎng)大致由盆地四周向盆內(nèi)擠壓，工區(qū)內(nèi)應(yīng)力場(chǎng)方向呈NW—SE向；喜馬拉雅期盆地最大主壓應(yīng)力場(chǎng)方向呈NNE—SSW向。因此，繞射波刻畫的斷裂主要與上述兩期構(gòu)造運(yùn)動(dòng)有關(guān)。

測(cè)井解釋成果表明：A-3井、A-7井的本溪組氣層厚度分別為1.7、1.9m，壓裂后均有產(chǎn)能；A-2井指示厚度為0.6m氣層和0.5m差氣層，壓裂后均無(wú)產(chǎn)能；其他井在本溪組無(wú)氣層。由于A-3井和A-7井均處于繞射斷裂發(fā)育帶，A-2井未處于斷裂發(fā)育帶(圖8)，因此推測(cè)產(chǎn)能與油氣充注和疏導(dǎo)情況有關(guān)，即氣層和斷裂同時(shí)存在時(shí)，往往有產(chǎn)能；在相干屬性中沒(méi)有這種規(guī)律(圖7)。上述認(rèn)識(shí)為井位部署提供了重要參考依據(jù)。

圖8 基于繞射波提取的本溪組儲(chǔ)層頂面螞蟻體屬性平面圖

4 結(jié)論

(1)應(yīng)用基于L0范數(shù)稀疏反演的CS寬頻地震資料處理方法提高數(shù)據(jù)分辨率，鑒于L0范數(shù)的稀疏特性，處理后剖面斷點(diǎn)清晰，斷層更易識(shí)別。

(2)與常規(guī)疊后屬性相比，基于繞射信息提取技術(shù)的斷裂識(shí)別方法能夠識(shí)別更小尺度的斷裂，在識(shí)別研究區(qū)致密砂巖薄儲(chǔ)層斷層、裂縫體系時(shí)具有獨(dú)特優(yōu)勢(shì)。

(3)對(duì)比鉆井氣層解釋結(jié)果、產(chǎn)能資料與繞射斷裂識(shí)別結(jié)果發(fā)現(xiàn)，斷裂和氣層同時(shí)存在時(shí)產(chǎn)能較高，該認(rèn)識(shí)可為井位部署提供重要依據(jù)。