发明名称 一种基于板形检测的冷轧带材跑偏自动识别方法
摘要 本发明公开了一种基于板形检测的冷轧带材跑偏自动识别方法,包括计算冷轧机出口带材边部前张应力变化量、判别带材是否跑偏以及跑偏方向和计算确定带材跑偏量三个步骤,利用安装于轧机出口接触式板形仪所检测到的前张应力反馈值,通过分析带材两侧边部前张应力反馈值变化量与相应板形执行机构反馈值变化量之间的对应关系,识别出得出带材的跑偏方向,并通过带材的宽度板形和板形仪通道之间的关系,计算出带材跑偏量,可以供人工手动控制或计算机自动控制跑偏,因此,本发明根据冷轧机出口带材边部前张应力反馈值及其相应板形执行机构反馈值的变化量,即可识别出带材的跑偏方向及跑偏量,从而不用另外增加设备,又比较实用可靠。
申请公布号 CN104275353B 申请公布日期 2016.08.24
申请号 CN201410487793.4 申请日期 2014.09.22
申请人 宁波宝新不锈钢有限公司 发明人 何汝迎;顾廷权;徐锋;乔爱民;袁意林
分类号 B21B39/14(2006.01)I;B21B38/02(2006.01)I 主分类号 B21B39/14(2006.01)I
代理机构 宁波诚源专利事务所有限公司 33102 代理人 袁忠卫
主权项 一种基于板形检测的冷轧带材跑偏自动识别方法,其特征在于:该识别方法根据冷轧机出口板形仪检测冷轧机出口带材边部前张应力反馈值及其相应板形执行机构反馈值,确定带材的跑偏方向及跑偏量,具体包括以下步骤:(1)计算当前监测周期相对于前一个监测周期带材边部前张应力以及各板形执行机构的变化量,相应计算模型为:当前监测周期相对于前一个监测周期带材操作侧边部前张应力变化量Δσ<sub>w</sub>(kT)的计算式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&Delta;&sigma;</mi><mi>w</mi></msub><mrow><mo>(</mo><mi>k</mi><mi>T</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&sigma;</mi><mi>w</mi></msub><mrow><mo>(</mo><mi>k</mi><mi>T</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&sigma;</mi><mi>w</mi></msub><mo>&lsqb;</mo><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>T</mi><mo>&rsqb;</mo><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>&Delta;S</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mi>T</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mfrac><mrow><mo>&part;</mo><msub><mi>&sigma;</mi><mi>w</mi></msub></mrow><mrow><mo>&part;</mo><msub><mi>S</mi><mi>i</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000901104000000011.GIF" wi="1070" he="143" /></maths>式中,k为当前监测周期序号,T为给定带材跑偏监测周期时间间隔,i=1、2、Λ、m,i表示为覆盖于板形仪宽向相应检测单元上各条带材的序号,m为自然数表示冷轧机板形执行机构数量,σ<sub>w</sub>(kT)为冷轧机出口板形仪检测的当前监测周期的带材操作侧边部前张应力反馈值,σ<sub>w</sub>[(k‑1)T]为冷轧机出口板形仪检测的前一个监测周期的带材操作侧边部前张应力反馈值,ΔS<sub>i</sub>(kT)当前监测周期相对于前一个监测周期冷轧机第i板形执行机构反馈值变化量,<img file="FDA0000901104000000012.GIF" wi="94" he="110" />为冷轧机第i板形执行机构对带材操作侧边部前张应力的影响系数,Δσ<sub>w</sub>(kT)为带材操作侧边部前张应力变化量,当前监测周期相对于前一个监测周期带材传动侧边部前张应力变化量Δσ<sub>d</sub>(kT)的计算式为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&Delta;&sigma;</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>k</mi><mi>T</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&sigma;</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>k</mi><mi>T</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&sigma;</mi><mi>d</mi></msub><mo>&lsqb;</mo><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>T</mi><mo>&rsqb;</mo><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>&Delta;S</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mi>T</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mfrac><mrow><mo>&part;</mo><msub><mi>&sigma;</mi><mi>d</mi></msub></mrow><mrow><mo>&part;</mo><msub><mi>S</mi><mi>i</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000901104000000013.GIF" wi="1069" he="140" /></maths>式中,σ<sub>d</sub>(kT)为冷轧机出口板形仪检测的当前监测周期的带材传动侧边部前张应力反馈值,σ<sub>d</sub>[(k‑1)T]为冷轧机出口板形仪检测的前一个监测周期的带材传动侧边部前张应力反馈值,<img file="FDA0000901104000000014.GIF" wi="95" he="119" />为冷轧机第i板形执行机构对带材传动侧边部前张应力的影响系数,Δσ<sub>d</sub>(kT)为带材传动侧边部前张应力变化量,当前监测周期相对于前一个监测周期冷轧机第i板形执行机构反馈值变化量ΔS<sub>i</sub>(kT)的计算式为:ΔS<sub>i</sub>(kT)=S<sub>i</sub>(kT)‑S<sub>i</sub>[(k‑1)T] i=1,2,Λ,m式中,S<sub>i</sub>(kT)当前监测周期冷轧机第i板形执行机构反馈值,S<sub>i</sub>[(k‑1)T]为前一个监测周期冷轧机第i板形执行机构反馈值;(2)判别带材是否跑偏以及跑偏方向,具体判别方式为:1)当带材边部所在板形仪检测通道的带材理论覆盖率x<设定的带材覆盖率x<sub>0</sub>时,此时,若带材操作侧边部前张应力变化量Δσ<sub>w</sub>(kT)≥K<sub>o</sub>·σ<sub>0</sub>,且带材传动侧边部前张应力变化量Δσ<sub>d</sub>(kT)≤K<sub>n</sub>·σ<sub>0</sub>,则带材向操作侧跑偏,若带材操作侧边部前张应力变化量Δσ<sub>w</sub>(kT)≤K<sub>n</sub>·σ<sub>0</sub>,且带材传动侧边部前张应力变化量Δσ<sub>d</sub>(kT)≥K<sub>o</sub>·σ<sub>0</sub>,则带材向传动侧跑偏,否则,带材视为没有跑偏,其中,x<sub>0</sub>为设定的带材覆盖率,K<sub>o</sub>表示带材边部前张应力增大速率,且大于0,K<sub>n</sub>表示带材边部前张应力减小速率且小于0,σ<sub>0</sub>为冷轧机出口带材平均前张应力,2)当带材边部所在板形仪检测通道的带材理论覆盖率x≥设定的带材覆盖率x<sub>0</sub>时,此时,若带材传动侧边部前张应力变化量Δσ<sub>d</sub>(kT)≤K<sub>n</sub>·σ<sub>0</sub>,则带材向操作侧跑偏,若带材操作侧边部前张应力变化量Δσ<sub>w</sub>(kT)≤K<sub>n</sub>·σ<sub>0</sub>,则带材向传动侧跑偏,否则,带材视为没有跑偏;(3)计算确定带材跑偏量,假设带材边部前张应力变化量是由带材跑偏以及冷轧机各板形执行机构引起的,故,<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>&sigma;</mi><mi>w</mi></msub><mrow><mo>(</mo><mi>k</mi><mi>T</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&sigma;</mi><mi>w</mi></msub><mo>&lsqb;</mo><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>T</mi><mo>&rsqb;</mo><mo>=</mo><mfrac><mrow><msub><mi>&Delta;x</mi><mi>w</mi></msub><mrow><mo>(</mo><mi>k</mi><mi>T</mi><mo>)</mo></mrow></mrow><mi>x</mi></mfrac><mo>&CenterDot;</mo><msub><mi>&sigma;</mi><mi>w</mi></msub><mo>&lsqb;</mo><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>T</mi><mo>&rsqb;</mo><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><mfrac><mrow><mo>&part;</mo><msub><mi>&sigma;</mi><mi>w</mi></msub></mrow><mrow><mo>&part;</mo><msub><mi>S</mi><mi>i</mi></msub></mrow></mfrac><mo>&CenterDot;</mo><msub><mi>&Delta;S</mi><mi>i</mi></msub><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&sigma;</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>k</mi><mi>T</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&sigma;</mi><mi>d</mi></msub><mo>&lsqb;</mo><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>T</mi><mo>&rsqb;</mo><mo>=</mo><mo>-</mo><mfrac><mrow><msub><mi>&Delta;x</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>k</mi><mi>T</mi><mo>)</mo></mrow></mrow><mi>x</mi></mfrac><mo>&CenterDot;</mo><msub><mi>&sigma;</mi><mi>d</mi></msub><mo>&lsqb;</mo><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>T</mi><mo>&rsqb;</mo><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><mfrac><mrow><mo>&part;</mo><msub><mi>&sigma;</mi><mi>d</mi></msub></mrow><mrow><mo>&part;</mo><msub><mi>S</mi><mi>i</mi></msub></mrow></mfrac><mo>&CenterDot;</mo><msub><mi>&Delta;S</mi><mi>i</mi></msub><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000901104000000021.GIF" wi="1381" he="299" /></maths>其中,Δx<sub>w</sub>(kT)和Δx<sub>d</sub>(kT)分别为带材跑偏引起的当前监测周期带材操作侧和传动侧边部所在板形仪检测通道覆盖率变化量,若Δx<sub>w</sub>(kT)&gt;0且Δx<sub>d</sub>(kT)&gt;0,则带材向操作侧跑偏,若Δx<sub>w</sub>(kT)&lt;0且Δx<sub>d</sub>(kT)&lt;0,则带材向传动侧跑偏,根据带材边部所在板形仪检测通道的带材理论覆盖率x的大小,计算带材跑偏量:当带材边部所在板形仪检测通道的带材理论覆盖率x<设定的带材覆盖率x<sub>0</sub>时,当前监测周期带材跑偏量的计算式为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>&Delta;x</mi><mi>w</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow><mo>=</mo><mi>x</mi><mo>&CenterDot;</mo><mrow><mo>{</mo><mrow><msub><mi>&sigma;</mi><mi>w</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow><mo>-</mo><msub><mi>&sigma;</mi><mi>w</mi></msub><mrow><mo>&lsqb;</mo><mrow><mrow><mo>(</mo><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>&rsqb;</mo></mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><mrow><mfrac><mrow><mo>&part;</mo><msub><mi>&sigma;</mi><mi>w</mi></msub></mrow><mrow><mo>&part;</mo><msub><mi>S</mi><mi>i</mi></msub></mrow></mfrac><mo>&CenterDot;</mo><msub><mi>&Delta;S</mi><mi>i</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow><mo>}</mo></mrow><mo>/</mo><msub><mi>&sigma;</mi><mi>w</mi></msub><mrow><mo>&lsqb;</mo><mrow><mrow><mo>(</mo><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>&rsqb;</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&Delta;x</mi><mi>d</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>x</mi><mo>&CenterDot;</mo><mrow><mo>{</mo><mrow><msub><mi>&sigma;</mi><mi>d</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow><mo>-</mo><msub><mi>&sigma;</mi><mi>d</mi></msub><mrow><mo>&lsqb;</mo><mrow><mrow><mo>(</mo><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>&rsqb;</mo></mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><mrow><mfrac><mrow><mo>&part;</mo><msub><mi>&sigma;</mi><mi>d</mi></msub></mrow><mrow><mo>&part;</mo><msub><mi>S</mi><mi>i</mi></msub></mrow></mfrac><mo>&CenterDot;</mo><msub><mi>&Delta;S</mi><mi>i</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow><mo>}</mo></mrow><mo>/</mo><msub><mi>&sigma;</mi><mi>d</mi></msub><mrow><mo>&lsqb;</mo><mrow><mrow><mo>(</mo><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>&rsqb;</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&Delta;y</mi><mi>w</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>l</mi><mi>w</mi></msub><mo>&CenterDot;</mo><msub><mi>&Delta;x</mi><mi>w</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&Delta;y</mi><mi>d</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>l</mi><mi>d</mi></msub><mo>&CenterDot;</mo><msub><mi>&Delta;x</mi><mi>d</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&Delta;y</mi><mi>a</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow><mo>=</mo><mrow><mo>&lsqb;</mo><mrow><msub><mi>l</mi><mi>w</mi></msub><mo>&CenterDot;</mo><msub><mi>&Delta;x</mi><mi>w</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow><mo>+</mo><msub><mi>l</mi><mi>d</mi></msub><mo>&CenterDot;</mo><msub><mi>&Delta;x</mi><mi>d</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow></mrow><mo>&rsqb;</mo></mrow><mo>/</mo><mn>2</mn></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000901104000000031.GIF" wi="1539" he="535" /></maths>式中,l<sub>w</sub>和l<sub>d</sub>分别为操作侧和传动侧带材边部所在板形仪检测通道的宽度,Δy<sub>w</sub>(kT)、Δy<sub>d</sub>(kT)和Δy<sub>a</sub>(kT)分别为当前监测周期带材操作侧跑偏量、传动侧跑偏量和平均跑偏量,当带材边部所在板形仪检测通道的带材理论覆盖率x≥设定的带材覆盖率x<sub>0</sub>时,若带材向操作侧跑偏,则当前监测周期带材跑偏量的计算式为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>&Delta;x</mi><mi>w</mi></msub><mrow><mo>(</mo><mi>k</mi><mi>T</mi><mo>)</mo></mrow><mo>=</mo><mi>x</mi><mo>&CenterDot;</mo><mo>{</mo><msub><mi>&sigma;</mi><mi>w</mi></msub><mrow><mo>(</mo><mi>k</mi><mi>T</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&sigma;</mi><mi>w</mi></msub><mo>&lsqb;</mo><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>T</mi><mo>&rsqb;</mo><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><mfrac><mrow><mo>&part;</mo><msub><mi>&sigma;</mi><mi>w</mi></msub></mrow><mrow><mo>&part;</mo><msub><mi>S</mi><mi>i</mi></msub></mrow></mfrac><mo>&CenterDot;</mo><msub><mi>&Delta;S</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>}</mo><mo>/</mo><msub><mi>&sigma;</mi><mi>w</mi></msub><mo>&lsqb;</mo><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>T</mi><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&Delta;y</mi><mi>a</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>&Delta;y</mi><mi>w</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>l</mi><mi>w</mi></msub><mo>&CenterDot;</mo><msub><mi>&Delta;x</mi><mi>w</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000901104000000032.GIF" wi="1421" he="231" /></maths>若带材向传动侧跑偏,则当前监测周期带材跑偏量的计算式为,<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mo>{</mo><mrow><mtable><mtr><mtd><mrow><msub><mi>&Delta;x</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>k</mi><mi>T</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>x</mi><mo>&CenterDot;</mo><mo>{</mo><msub><mi>&sigma;</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>k</mi><mi>T</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&sigma;</mi><mi>d</mi></msub><mo>&lsqb;</mo><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>T</mi><mo>&rsqb;</mo><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><mfrac><mrow><mo>&part;</mo><msub><mi>&sigma;</mi><mi>d</mi></msub></mrow><mrow><mo>&part;</mo><msub><mi>S</mi><mi>i</mi></msub></mrow></mfrac><mo>&CenterDot;</mo><msub><mi>&Delta;S</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>}</mo><mo>/</mo><msub><mi>&sigma;</mi><mi>d</mi></msub><mo>&lsqb;</mo><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>T</mi><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&Delta;y</mi><mi>a</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>&Delta;y</mi><mi>d</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>l</mi><mi>d</mi></msub><mo>&CenterDot;</mo><msub><mi>&Delta;x</mi><mi>d</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mi>T</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>.</mo></mrow></mrow>]]></math><img file="FDA0000901104000000033.GIF" wi="1476" he="246" /></maths>
地址 315807 浙江省宁波市经济技术开发区(北仑霞浦)