发明名称 水下埋设管线受邻近施工振动影响的安全监测评估方法
摘要 本发明公开了一种水下埋设管线受邻近施工振动影响的安全监测评估方法,具体过程为:在施工振动位置距埋设管线最近的直线上及沿管线方向确定测线布置相应振动监测点,利用海上施工船舶进行水上定位、水下钻孔并埋设振动传感器,在工程施工全过程中通过动态信号采集分析系统,实时采集监测的振动数据,结合土体-混凝土管线的局部精细化模型得出管线安全状态的评估值,评估管线安全。本发明建立邻近水下埋设管线在施工振动下的振动响应的监测体系,安全准确的监测近海工程中邻近水下埋设管线在施工振动下的振动响应,并建立局部精细化模型结合监测数据对邻近水下埋设管线安全评估与预警。
申请公布号 CN103835278B 申请公布日期 2015.09.02
申请号 CN201310703471.4 申请日期 2013.12.19
申请人 河海大学 发明人 苏静波;诸裕良;张栋;张龙;朱瑞虎;朱义欢;马志国
分类号 E02D1/00(2006.01)I 主分类号 E02D1/00(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 朱小兵
主权项 一种水下埋设管线受邻近施工振动影响的安全监测评估方法,所述邻近施工为近海工程施工,其特征在于,包括以下步骤:步骤1:确定监测路径、监测点位置及监测点数量;步骤2:根据工程资料,通过有限元软件建立土体‑混凝土管线的局部精细化模型,并定义其损伤本构关系,最终得出管线的起始损伤速度;土体‑混凝土管线的局部精细化模型损伤本构关系的具体表达式如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>&sigma;</mi><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>d</mi><mo>)</mo></mrow><mover><mi>&sigma;</mi><mo>&OverBar;</mo></mover></mrow>]]></math><img file="FDA0000718832900000011.GIF" wi="264" he="69" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>&sigma;</mi><mo>&OverBar;</mo></mover><mo>=</mo><msubsup><mi>D</mi><mn>0</mn><mi>el</mi></msubsup><mrow><mo>(</mo><mi>&epsiv;</mi><mo>-</mo><msup><mi>&epsiv;</mi><mi>pl</mi></msup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000718832900000012.GIF" wi="327" he="80" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mover><mi>&epsiv;</mi><mover><mo>~</mo><mo>.</mo></mover></mover><mi>pl</mi></msup><mo>=</mo><mi>h</mi><mrow><mo>(</mo><mover><mi>&sigma;</mi><mo>&OverBar;</mo></mover><mo>,</mo><msup><mover><mi>&epsiv;</mi><mo>~</mo></mover><mi>pl</mi></msup><mo>)</mo></mrow><mo>&CenterDot;</mo><msup><mover><mi>&epsiv;</mi><mo>.</mo></mover><mi>pl</mi></msup></mrow>]]></math><img file="FDA0000718832900000013.GIF" wi="364" he="73" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msup><mover><mi>&epsiv;</mi><mo>.</mo></mover><mi>pl</mi></msup><mo>=</mo><mi>&lambda;</mi><mfrac><mrow><mo>&PartialD;</mo><mi>G</mi><mrow><mo>(</mo><mover><mi>&sigma;</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><mover><mi>&sigma;</mi><mo>&OverBar;</mo></mover></mrow></mfrac></mrow>]]></math><img file="FDA0000718832900000014.GIF" wi="299" he="129" /></maths>式中,σ为名义应力;d为总的损伤指标;<img file="FDA0000718832900000015.GIF" wi="50" he="59" />为有效应力;<img file="FDA0000718832900000016.GIF" wi="82" he="78" />为材料的弹性本构矩阵;ε为材料应变;ε<sup>pl</sup>为塑性应变;h为硬化函数;<img file="FDA0000718832900000017.GIF" wi="65" he="75" />为等效塑性应变率;<img file="FDA0000718832900000018.GIF" wi="68" he="68" />为等效塑性应变;<img file="FDA0000718832900000019.GIF" wi="67" he="66" />为塑性应变率;λ为塑性乘子;G为塑性势函数;在模型中引入刚度退化指标,通过以下公式定义总的损伤指标:(1‑d)=(1‑s<sub>t</sub>d<sub>c</sub>)(1‑s<sub>c</sub>d<sub>t</sub>) 0≤s<sub>t</sub>,s<sub>c</sub>≤1<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><msub><mi>s</mi><mi>t</mi></msub><mo>=</mo><mn>1</mn><mo>-</mo><msub><mi>&omega;</mi><mi>t</mi></msub><msup><mi>r</mi><mo>*</mo></msup><mrow><mo>(</mo><msub><mover><mi>&sigma;</mi><mo>&OverBar;</mo></mover><mn>11</mn></msub><mo>)</mo></mrow></mtd><mtd><mn>0</mn><mo>&le;</mo><msub><mi>&omega;</mi><mi>t</mi></msub><mo>&le;</mo><mn>1</mn></mtd></mtr></mtable></mfenced>]]></math><img file="FDA00007188329000000110.GIF" wi="840" he="74" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><msub><mi>s</mi><mi>c</mi></msub><mo>=</mo><mn>1</mn><mo>-</mo><msub><mi>&omega;</mi><mi>c</mi></msub><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>r</mi><mo>*</mo></msup><mrow><mo>(</mo><msub><mover><mi>&sigma;</mi><mo>&OverBar;</mo></mover><mn>11</mn></msub><mo>)</mo></mrow><mo>)</mo></mtd><mtd><mn>0</mn><mo>&le;</mo><msub><mi>&omega;</mi><mi>c</mi></msub><mo>&le;</mo><mn>1</mn></mtd></mtr></mtable></mfenced>]]></math><img file="FDA00007188329000000111.GIF" wi="841" he="76" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msup><mi>r</mi><mo>*</mo></msup><mrow><mo>(</mo><msub><mover><mi>&sigma;</mi><mo>&OverBar;</mo></mover><mn>11</mn></msub><mo>)</mo></mrow><mo>=</mo><mi>H</mi><mrow><mo>(</mo><msub><mover><mi>&sigma;</mi><mo>&OverBar;</mo></mover><mn>11</mn></msub><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>if</mi></mtd><mtd><msub><mover><mi>&sigma;</mi><mo>&OverBar;</mo></mover><mn>11</mn></msub><mo>></mo><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>if</mi></mtd><mtd><msub><mover><mi>&sigma;</mi><mo>&OverBar;</mo></mover><mn>11</mn></msub><mo>&lt;</mo><mn>0</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA00007188329000000112.GIF" wi="861" he="164" /></maths>式中,d<sub>t</sub>为受拉损伤指标;d<sub>c</sub>为受压损伤指标;s<sub>t</sub>和s<sub>c</sub>分别为应力状态的函数,用来描述应力状态改变即拉压互相转换时对刚度退化的影响;ω<sub>t</sub>和ω<sub>c</sub>分别为权重因子;<img file="FDA00007188329000000113.GIF" wi="77" he="68" />为施加冲击速度作用的X方向的有效应力,X方向指在水平面内施工振动位置与该振动监测点的连线方向;<img file="FDA00007188329000000114.GIF" wi="156" he="78" />为自变量<img file="FDA00007188329000000115.GIF" wi="72" he="70" />的Heaviside函数;<img file="FDA00007188329000000116.GIF" wi="158" he="86" />为用于计算应力状态的函数,通过<img file="FDA00007188329000000117.GIF" wi="159" he="77" />来赋值;在土体‑混凝土管线的局部精细化模型边界施加冲击速度,计算冲击速度变化过程中的损伤情况;若总的损伤指标d&gt;0,即管线出现损伤,此时的冲击速度即为管线的起始损伤速度;步骤3:利用海上施工船舶进行水上定位,在预定振动监测位置,进行水下钻孔并埋设振动传感器;步骤4:在工程施工全过程中实时采集监测的振动数据,结合管线的起始损伤速度来计算管线安全状态的评估值θ<sub>min</sub>,评估管线安全,具体如下:401,对第i个监测点,管线的起始损伤速度记为V<sub>i</sub>,根据监测的振动数据所得振动峰值为V<sub>i</sub>',计算管线安全状态的相对值θ<sub>i</sub>,具体公式为:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>&theta;</mi><mi>i</mi></msub><mo>=</mo><mfrac><msub><mi>V</mi><mi>i</mi></msub><msup><msub><mi>V</mi><mi>i</mi></msub><mo>&prime;</mo></msup></mfrac></mrow>]]></math><img file="FDA0000718832900000021.GIF" wi="174" he="144" /></maths>式中,i为自然数;402,管线安全状态评估值选取各监测点的管线安全状态的相对值的最小值:θ<sub>min</sub>=min(θ<sub>1</sub>,θ<sub>2</sub>,…,θ<sub>i</sub>,…)则θ<sub>min</sub>即为管线目前所处的总体安全状态,用以评估近海工程施工振动对邻近水下埋设管线安全影响;若θ<sub>min</sub>小于1,表明管线出现损伤;否则表明管线处于安全状态。
地址 210098 江苏省南京市鼓楼区西康路1号