DOI：10.3969/j.issn.1001-2206.2024.06.016

摘" " 要：打孔盜油導(dǎo)致了國(guó)有資產(chǎn)流失，同時(shí)給生態(tài)環(huán)境和人身財(cái)產(chǎn)安全帶來嚴(yán)重威脅，因此需要對(duì)打孔盜油事件發(fā)生后的可能危害進(jìn)行充分的研究與分析，合理評(píng)估風(fēng)險(xiǎn)水平。以獨(dú)烏成品油管道某處打孔盜油問題為例，針對(duì)事故發(fā)生后的潛在風(fēng)險(xiǎn)事件樹，詳細(xì)分析了氣象數(shù)據(jù)、沿途人口分布、點(diǎn)火源及點(diǎn)火概率等相關(guān)參數(shù)，給出了后果估計(jì)，對(duì)成品油管道打孔盜油的個(gè)人風(fēng)險(xiǎn)和社會(huì)風(fēng)險(xiǎn)進(jìn)行了定量評(píng)估研究。結(jié)果表明，打孔盜油事件會(huì)明顯增大成品油管道對(duì)周圍人員的風(fēng)險(xiǎn)，致使個(gè)人風(fēng)險(xiǎn)和社會(huì)風(fēng)險(xiǎn)均超出相關(guān)標(biāo)準(zhǔn)要求的限值。

關(guān)鍵詞：打孔盜油；定量風(fēng)險(xiǎn)評(píng)估；成品油管道；個(gè)人風(fēng)險(xiǎn)；社會(huì)風(fēng)險(xiǎn)

Abstract：Stealing oil by drilling holes in the product pipeline can lead to the loss of state-owned assets， and bring serious threats to the ecological environment， personal safety， and property security. Therefore， it is necessary to carry out sufficient research and analysis on the possible damage after the occurrence of stealing oil by drilling holes and to reasonably assess the risk level. Taking stealing oil by drilling holes occurred in the Dushanzi-Urumqi product oil pipeline as an example， this paper analyzes the meteorological data， population distribution along the route， ignition source， and ignition probability in detail with the consideration of the potential risk event tree after the accident. Then， this paper estimates the consequences and quantitatively evaluates the individual risk and social risk brought by this behavior. The results show that stealing oil by drilling holes can significantly increase the risk faced by personnel around the pipeline， raising the individual risk and social risk beyond the value limited by related standards.

Keywords：stealing oil by drilling holes; quantitative risk assessment; product oil pipeline; individual risk; social risk

我國(guó)油氣管道總長(zhǎng)度已突破千萬公里，但油品盜竊問題一直存在[1]。打孔盜油不僅竊取了國(guó)有資產(chǎn)，而且盜取油品時(shí)極易發(fā)生泄漏，造成環(huán)境污染，對(duì)生態(tài)文明造成長(zhǎng)期破壞。此外，泄露的油品具有易燃易爆特性，有可能引起火災(zāi)、爆炸等災(zāi)害，對(duì)周圍建筑與人民群眾造成巨大的生命財(cái)產(chǎn)威脅[2]。另一方面，長(zhǎng)輸管道一旦遭遇打孔盜油事故，或?qū)⒚媾R全線停輸，給石油企業(yè)和國(guó)民經(jīng)濟(jì)、生活的正常運(yùn)行帶來阻礙[3]。研究表明，我國(guó)石油公司每年因油品被盜的損失約30多億元，僅大慶油田每年就損失10億元 [4]。管道油品盜竊是一個(gè)嚴(yán)重的全球性問題，在過去幾年中一直呈上升趨勢(shì)[5]。

打孔盜油問題具有散落性、隱蔽性、多發(fā)性、季節(jié)性、流竄性等特點(diǎn)，且犯罪手段呈日益專業(yè)化的趨勢(shì)[3，6]。當(dāng)前國(guó)內(nèi)外對(duì)該問題開展了一些研究，并在綜合治理方面提出了一些較為先進(jìn)的技術(shù)手段[7-8]，例如韓永波等[9]開發(fā)了一種基于BP神經(jīng)網(wǎng)絡(luò)模型的打孔盜油事件預(yù)測(cè)模型，Sun等[4]分析了我國(guó)原油盜竊的多樣性和石油防盜的難點(diǎn)，并進(jìn)行了需求分析，設(shè)計(jì)了一種基于物聯(lián)網(wǎng)應(yīng)用的智能防盜系統(tǒng)，能及時(shí)制止盜油行為并提前防止不法分子竊取原油。盡管國(guó)內(nèi)外各石油公司已經(jīng)采取了多種措施和技術(shù)來預(yù)防打孔盜油事件的發(fā)生，但始終無法完全消除該問題。因此，需要對(duì)打孔盜油事件發(fā)生后的危害進(jìn)行充分的研究與分析，合理評(píng)估風(fēng)險(xiǎn)水平，進(jìn)而給出各類處理預(yù)案。

本文以獨(dú)烏成品油管道某處打孔盜油問題為例，對(duì)成品油管道打孔盜油的個(gè)人風(fēng)險(xiǎn)和社會(huì)風(fēng)險(xiǎn)進(jìn)行定量評(píng)估研究。

1" " 風(fēng)險(xiǎn)標(biāo)準(zhǔn)和風(fēng)險(xiǎn)事件樹

風(fēng)險(xiǎn)是某一特定危險(xiǎn)事件在特定時(shí)期或給定情況下發(fā)生的概率，而定量風(fēng)險(xiǎn)評(píng)估則是指把風(fēng)險(xiǎn)大小以量化的形式體現(xiàn)出來的風(fēng)險(xiǎn)評(píng)價(jià)方法。因此，首先需要明確具有潛在危害的裝置或行為。

1.1" " 風(fēng)險(xiǎn)標(biāo)準(zhǔn)

風(fēng)險(xiǎn)標(biāo)準(zhǔn)是評(píng)判風(fēng)險(xiǎn)是否可接受的尺碼，或者對(duì)其重要性做出一些其他判斷。依照設(shè)定的標(biāo)準(zhǔn)比較風(fēng)險(xiǎn)級(jí)別并對(duì)其主要風(fēng)險(xiǎn)因子進(jìn)行識(shí)別，這一過程即為風(fēng)險(xiǎn)評(píng)價(jià)[10]。本文以特定場(chǎng)所個(gè)人風(fēng)險(xiǎn)和社會(huì)（群體）風(fēng)險(xiǎn)作為風(fēng)險(xiǎn)標(biāo)準(zhǔn)進(jìn)行研究。

特定場(chǎng)所個(gè)人風(fēng)險(xiǎn)表示界區(qū)外某一個(gè)體持續(xù)出現(xiàn)在某一特定場(chǎng)所所遭遇的某種危險(xiǎn)發(fā)生的頻率，通常以年度個(gè)人死亡率表示。對(duì)該風(fēng)險(xiǎn)進(jìn)行評(píng)估時(shí)假設(shè)暴露個(gè)體無法逃生、在室外且無保護(hù)。

社會(huì)（群體）風(fēng)險(xiǎn)表示群體在危險(xiǎn)區(qū)域承受某種程度傷害的頻發(fā)程度，可以采用F-N曲線和潛在損失生命來表征[11]。F-N曲線是用來明確表示累積頻率（F）和死亡人數(shù)（N）之間關(guān)系的曲線圖，通常使用對(duì)數(shù)坐標(biāo)圖。社會(huì)死亡風(fēng)險(xiǎn)描述的是可能受到災(zāi)難性事件影響的死亡人數(shù)。對(duì)位于不可接受區(qū)的風(fēng)險(xiǎn)，需要提出可行的風(fēng)險(xiǎn)減緩措施。在整個(gè)研究過程中通過反復(fù)進(jìn)行定量風(fēng)險(xiǎn)評(píng)估，可以評(píng)判這些措施的有效性。

1.2" " 風(fēng)險(xiǎn)事件樹

目標(biāo)管道輸送成品油介質(zhì)，泄漏出來的輕組分與空氣混合能形成爆炸性混合物，遇熱源或明火時(shí)有燃燒爆炸的危險(xiǎn)；成品油泄漏可能會(huì)形成池火或者流淌火[11]；另外成品油泄漏會(huì)污染周圍的土壤和水源，因此該處打孔盜油的潛在風(fēng)險(xiǎn)事件樹如圖1所示。

池火或流淌火是指液體或者液化碳?xì)浠衔镌诳拷c(diǎn)火源處泄漏，可燃?xì)怏w空氣混合物中的蒸發(fā)氣體會(huì)在液池之上燃燒。

對(duì)于閃火事件，最常見的情況是可燃蒸氣云團(tuán)在擴(kuò)散時(shí)遇到一個(gè)強(qiáng)點(diǎn)火源，在云團(tuán)的邊緣處被點(diǎn)燃，火焰前緣會(huì)反向迅速回燃，直到整個(gè)氣云團(tuán)被燃燒殆盡。當(dāng)空氣中的可燃云團(tuán)濃度處于可燃下限和可燃上限濃度之間時(shí)，任何一處的擴(kuò)散云團(tuán)都能被點(diǎn)燃。

蒸氣云爆炸是指碳?xì)浠衔镎魵庠诳扇紳舛确秶鷥?nèi)，在受限或者部分受限空間內(nèi)被點(diǎn)燃而發(fā)生的爆炸。

2" " 評(píng)估計(jì)算參數(shù)與指標(biāo)參數(shù)

打孔盜油事件發(fā)生后的可能危害與該事件中長(zhǎng)輸管道的工作參數(shù)及周圍環(huán)境條件等諸多參數(shù)有關(guān)[12]，要進(jìn)行定量風(fēng)險(xiǎn)評(píng)估需要對(duì)這些參數(shù)進(jìn)行綜合分析。本文以獨(dú)烏成品油管道某處打孔盜油問題為例進(jìn)行分析。

2.1" " 打孔盜油事故

本文研究的事例中，打孔支管直徑為1.5 in（1 in=25.4 mm），獨(dú)烏成品油管道基本參數(shù)見表1。該打孔盜油處周圍有一排毛坯土房，其中一戶有人居住，見圖2。自接管以來，共發(fā)生兩起打孔盜油事件，發(fā)生頻率為6年一次，折合年發(fā)生概率為0.17。

2.2" " 大氣環(huán)境參數(shù)

大氣環(huán)境通常通過實(shí)時(shí)變化的多個(gè)參數(shù)進(jìn)行描述，本文在詳細(xì)分析了該事故位置附近的環(huán)境數(shù)據(jù)后，給出了圖3所示的打孔盜油處所在地區(qū)的風(fēng)頻數(shù)據(jù)圖。

為了進(jìn)一步對(duì)事故風(fēng)險(xiǎn)進(jìn)行定量化評(píng)估，本文選取了事故所在地區(qū)兩組典型的氣象條件進(jìn)行研究，表2給出了白天及晚上大氣穩(wěn)定度和風(fēng)速的氣象數(shù)據(jù)，其中，白天和晚間平均大氣溫度分別為15.3 ℃和2.0 ℃，大氣壓力為1.029 bar（1 bar=100 kPa），太陽輻射為0.5 kW/m2。

2.3" " 點(diǎn)火源及點(diǎn)火概率

由于點(diǎn)火源的存在，管道泄漏后可能會(huì)引發(fā)閃火和池火風(fēng)險(xiǎn)。本例中周圍存在的確定點(diǎn)火源主要是來往的交通車輛和周圍的人口活動(dòng)，事故周圍的點(diǎn)火源分布位置如圖4所示。

根據(jù)燃燒后果的不同，對(duì)立即點(diǎn)火和延遲點(diǎn)火兩種模型進(jìn)行分析，其中閃火只可能是易燃蒸氣混合物經(jīng)過延遲點(diǎn)火而產(chǎn)生的，而立即點(diǎn)火只能引發(fā)池火事故[12]。

2.3.1" " 立即點(diǎn)火概率

立即點(diǎn)火概率由SAFETI軟件自帶的點(diǎn)火模型選取，其來源于荷蘭應(yīng)用科學(xué)研究院（TNO）的點(diǎn)火模型，目前在行業(yè)具有權(quán)威性。

2.3.2" " 延遲點(diǎn)火

延遲點(diǎn)火概率既可以使用SAFETI軟件中預(yù)定義的點(diǎn)火源信息，也可在軟件里直接定義。對(duì)每一個(gè)泄漏事件，計(jì)算出來的延遲點(diǎn)火概率值會(huì)隨著天氣類別、風(fēng)速和風(fēng)向、點(diǎn)火源的數(shù)量、強(qiáng)度和與擴(kuò)散可燃蒸氣云團(tuán)的接觸時(shí)間的不同而不同。本文考慮了點(diǎn)火源特性、泄漏物組分及點(diǎn)火源處于蒸氣云團(tuán)內(nèi)的概率。在任意時(shí)刻t∈（0， T］發(fā)生點(diǎn)火的概率p（ t ）由下式給出[13]。

式中：ω表示點(diǎn)火效率，與點(diǎn)火源特性有關(guān)，可根據(jù)點(diǎn)火源在某一段時(shí)間內(nèi)的點(diǎn)火概率計(jì)算得出。

本文中不同火源特性的點(diǎn)火概率計(jì)算可以參考文獻(xiàn)[13]進(jìn)行選取。

2.4" " 評(píng)估計(jì)算指標(biāo)參數(shù)

根據(jù)圖1所示的事件樹可知，成品油管道發(fā)生打孔盜油事故后形成的后果類型包括閃火、爆炸和池火3種。

閃火事故通常持續(xù)時(shí)間短，不存在超壓，但對(duì)室外人員致死率為100%，因此需判定造成閃火的氣體云團(tuán)所覆蓋的范圍。

池火產(chǎn)生的熱輻射強(qiáng)度指標(biāo)參數(shù)在表3中給出，其中熱輻射強(qiáng)度由入射通量I表示，其表達(dá)式為：I=Qtc/（ 4πR2 ）。式中：Q為總熱輻射通量，tc為熱傳導(dǎo)系數(shù)，R為目標(biāo)點(diǎn)到液池中心距離。

本文在分析蒸氣云爆炸產(chǎn)生的沖擊波對(duì)建筑物的影響時(shí)，按照文獻(xiàn)[13]（表H.3）的規(guī)定判斷超壓對(duì)建筑物的影響。

3" " 后果風(fēng)險(xiǎn)評(píng)估

3.1" " 評(píng)估計(jì)算方案流程

在全面獲取第1.1節(jié)所述打孔盜油事件中該長(zhǎng)輸管道的工作參數(shù)以及周圍環(huán)境條件等諸多參數(shù)并完成數(shù)據(jù)收集后，對(duì)該打孔盜油事件進(jìn)行定量風(fēng)險(xiǎn)評(píng)估；對(duì)潛在危險(xiǎn)進(jìn)行辨識(shí)后，通過失效頻率分析與失效后果分析，對(duì)潛在風(fēng)險(xiǎn)事件樹進(jìn)行計(jì)算；最后，詳細(xì)研究該事件帶來的個(gè)人風(fēng)險(xiǎn)和社會(huì)風(fēng)險(xiǎn)，評(píng)估方案流程如圖5所示。

對(duì)于個(gè)人風(fēng)險(xiǎn)等高線的評(píng)估計(jì)算流程為：第一步，獲取打孔盜油事故的年發(fā)生概率，并對(duì)收集到的環(huán)境參數(shù)進(jìn)行計(jì)算，得到圖3所示的風(fēng)頻數(shù)據(jù)圖；第二步，根據(jù)潛在風(fēng)險(xiǎn)事件樹中的事故類型設(shè)置點(diǎn)火概率；第三步，計(jì)算在特定環(huán)境參數(shù)與特定后果事件下網(wǎng)格單元上的個(gè)體死亡概率，然后對(duì)所有的后果類型以及所有環(huán)境參數(shù)重復(fù)以上計(jì)算過程；第四步，取以上計(jì)算得到的概率之和，并得到個(gè)人風(fēng)險(xiǎn)等高線。

對(duì)于社會(huì)風(fēng)險(xiǎn)等高線的評(píng)估計(jì)算流程，前兩步與個(gè)人風(fēng)險(xiǎn)的計(jì)算一致。第三步則先計(jì)算在特定環(huán)境參數(shù)與特定后果事件下網(wǎng)格單元內(nèi)的人口死亡百分比和死亡人數(shù)（人口死亡百分比和網(wǎng)格單元內(nèi)人數(shù)的乘積），得到死亡總?cè)藬?shù)（所有網(wǎng)格單元內(nèi)死亡人數(shù)之和）和發(fā)生該特定后果的累積頻率（各條件的概率/頻率的乘積），然后對(duì)所有的后果類型以及所有環(huán)境參數(shù)重復(fù)以上計(jì)算過程；第四步，繪制死亡總?cè)藬?shù)和累積頻率的對(duì)應(yīng)關(guān)系圖，得到打孔盜油事件后的社會(huì)風(fēng)險(xiǎn)F-N曲線。

3.2" " 池火熱輻射風(fēng)險(xiǎn)計(jì)算

打孔盜油處不同天氣條件下形成的池火熱輻射影響距離如圖6所示。

3.3" " 閃火風(fēng)險(xiǎn)計(jì)算

成品油形成的閃火影響范圍與液池影響范圍基本一致，經(jīng)計(jì)算得到打孔盜油形成的液池影響范圍如圖7所示。

不同天氣條件下液池影響范圍如表4所示。

3.4" " 蒸氣云爆炸和環(huán)境污染

該打孔盜油處周圍相對(duì)開闊，不存在高擁塞度區(qū)域，因此不會(huì)引發(fā)較大超壓的爆炸，爆炸風(fēng)險(xiǎn)較小。

此外，成品油泄漏后如未點(diǎn)燃則會(huì)對(duì)范圍內(nèi)的土壤和水源等環(huán)境造成污染，該打孔盜油處地處戈壁區(qū)域，環(huán)境污染風(fēng)險(xiǎn)較小。

3.5" " 個(gè)人風(fēng)險(xiǎn)等高線

該段管道未發(fā)生打孔盜油事件和發(fā)生該事故后的個(gè)人風(fēng)險(xiǎn)等高線影響范圍見圖8和圖9。圖中各顏色曲線表示產(chǎn)生的個(gè)人風(fēng)險(xiǎn)等高線，各等高線具體對(duì)應(yīng)數(shù)值見圖例。

從上述計(jì)算結(jié)果可以看出：其一，由圖8和圖9的結(jié)果可知，在未發(fā)生打孔盜油事件前該管道產(chǎn)生的最大個(gè)人風(fēng)險(xiǎn)等高線為1×10-6，在發(fā)生打孔盜油事件后該管道產(chǎn)生的最大個(gè)人風(fēng)險(xiǎn)等高線上升為1×10-2，打孔盜油事件明顯增大了該段管道的個(gè)人風(fēng)險(xiǎn)。其二，在未發(fā)生打孔盜油事件前管道周圍毛坯土房住戶的個(gè)人風(fēng)險(xiǎn)值為7.77×10-7，滿足SY/T 6859—2020對(duì)個(gè)人風(fēng)險(xiǎn)的限值要求；在發(fā)生打孔盜油事件后管道周圍毛坯土房住戶的個(gè)人風(fēng)險(xiǎn)值上升為1.1×10-2，超出SY/T 6859—2020對(duì)個(gè)人風(fēng)險(xiǎn)的限值要求。

3.6" " 社會(huì)風(fēng)險(xiǎn)評(píng)估

該段管道未發(fā)生打孔盜油事件的社會(huì)風(fēng)險(xiǎn)計(jì)算結(jié)果（F-N曲線）見圖10，發(fā)生打孔盜油事件后的社會(huì)風(fēng)險(xiǎn)計(jì)算結(jié)果（F-N曲線）見圖11。

圖中黃色線表示風(fēng)險(xiǎn)可接受線，綠色線表示產(chǎn)生的社會(huì)風(fēng)險(xiǎn)曲線。從計(jì)算結(jié)果可知，在未發(fā)生打孔盜油事件前，該段管道的F-N曲線均位于風(fēng)險(xiǎn)可接受區(qū)域，滿足行業(yè)標(biāo)準(zhǔn)SY/T 6859—2020中規(guī)定的社會(huì)風(fēng)險(xiǎn)基準(zhǔn)；在發(fā)生打孔盜油事件后，該段管道的F-N曲線均位于不可接受區(qū)域，不滿足行業(yè)標(biāo)準(zhǔn)SY/T 6859—2020中規(guī)定的社會(huì)風(fēng)險(xiǎn)基準(zhǔn)。

4" " 結(jié)論

本文以獨(dú)烏成品油管道某處打孔盜油問題為例，針對(duì)事故發(fā)生后的潛在風(fēng)險(xiǎn)事件樹，詳細(xì)分析了氣象數(shù)據(jù)、沿途人口分布、點(diǎn)火源及點(diǎn)火概率等相關(guān)參數(shù)，給出了后果估計(jì)，對(duì)成品油管道打孔盜油的個(gè)人風(fēng)險(xiǎn)和社會(huì)風(fēng)險(xiǎn)進(jìn)行定量評(píng)估研究。通過本文的研究，可以得到如下結(jié)論。

1）打孔盜油事件明顯增大了成品油管道的個(gè)人風(fēng)險(xiǎn)，發(fā)生打孔盜油事件后管道周圍毛坯土房住戶的個(gè)人風(fēng)險(xiǎn)值不滿足SY/T 6859—2020對(duì)個(gè)人風(fēng)險(xiǎn)的限值要求。

2）打孔盜油事件明顯增大了成品油管道的社會(huì)風(fēng)險(xiǎn)，發(fā)生打孔盜油事件后該段管道的F-N曲線均位于不可接受區(qū)域，不滿足SY/T 6859—2020中規(guī)定的社會(huì)風(fēng)險(xiǎn)基準(zhǔn)。

綜上所述，打孔盜油事件會(huì)明顯增大成品油管道對(duì)周圍人員的風(fēng)險(xiǎn)，致使個(gè)人風(fēng)險(xiǎn)和社會(huì)風(fēng)險(xiǎn)均超出相關(guān)標(biāo)準(zhǔn)要求的限值。

參考文獻(xiàn)

[1]" 王春枝，孫兆光. “12·19”大案懸疑：盜油何時(shí)方能休[J]. 中國(guó)石油石化，2004（12）：41-43.

[2]" 賈玉高. 長(zhǎng)輸管道遭受打孔盜油的危害及對(duì)策探討[J]. 安全、健康環(huán)境，2006（6）：44-45.

[3]" 陳振友. 輸油管線打孔盜油作案的特點(diǎn)分析及防范措施[J]. 石油庫與加油站，2010，19（3）：12-14，47.

[4]" SUN J F，ZHANG Z Y，SUN X L. The intelligent crude oil anti-theft system based on iot under different scenarios[J]. Procedia Computer Science，2016，96：1 581-1 588.

[5]" COLLINS H W，ELEKWACHI W. Crude oil theft in the Niger Delta： the oil companies and host communities conundrum[J]. International Journal of Research and Scientific Innovation，2020，7：22-32.

[6]" 常番. “打孔盜油”的綜合治理研究——以河北地區(qū)為研究視角[D]. 北京：中國(guó)人民公安大學(xué)，2020.

[7]" 齊新春. 加強(qiáng)綜合治理防范打孔盜油——淺談防范打孔盜油案件的有效途徑[J]. 化工管理，2013（ 8 ）：20，22.

[8]" 王欣. 長(zhǎng)輸管道打孔盜油案件的偵查與防范策略[J]. 四川警察學(xué)院學(xué)報(bào)，2016，28（6）：31-36.

[9]" 韓永波，張鵬宇，田娜. 基于大數(shù)據(jù)和神經(jīng)網(wǎng)絡(luò)的油氣管道打孔盜油（氣）事件預(yù)測(cè)方法研究[J]. 中外能源，2020，25（8）：41-44.

[10] 尚騰飛，吳達(dá)嶺，齊學(xué)振. 失效頻率修正對(duì)風(fēng)險(xiǎn)結(jié)果的影響分析[J]. 廣州化工，2022，50（11）：128-130，154.

[11] 魏超南，孫舒，趙俊丹，等. 關(guān)于某?；肺锪鞴镜亩匡L(fēng)險(xiǎn)評(píng)估案例與思考[J]. 安全、健康和環(huán)境，2017，17（11）：34-37.

[12] 楊光輝，宋漢成. 高后果區(qū)高壓燃?xì)夤艿朗鹿屎蠊绊懸蛩匮芯縖J]. 科技通報(bào)，2023，39（3）：89-95，102.

[13] AQ/T 3046—2013，化工企業(yè)定量風(fēng)險(xiǎn)評(píng)價(jià)導(dǎo)則[S].

作者簡(jiǎn)介：馬建文（1993—），男，甘肅莊浪人，2018年畢業(yè)于西安石油大學(xué)材料物理專業(yè)，現(xiàn)從事生產(chǎn)運(yùn)行管理方面的工作。

Email：majw1993@126.com

收稿日期：2024-08-23