nsga2的源代码

比较实用的常用低速接口verilog代码，包含了Uart，SPI, I2C等，可以为自己的开发设计节约时间，可以直接调用

使用matlab编写的元胞自动机模拟交通路况，双车道，有超车模型，参数可改

三维孔洞储层建模及其地震波场正演模拟，理论讲解很透彻，分析思路清晰1290地球物理学进展26卷预测结果,即从具有确定性资料的控制点(如井点)解释.同时,利用该过程中产生的6冂井的时深关出发,推测出点间(如井间)确定的、唯一的储层参系,用三角剖分网格建立了速度模型并经过井点校数随机模拟是从一个随机函数z(v)中抽取多个可正,实现了对工区构造框架的时深转换通过以波阻能实现,即人工合成反映Z()空间分布的可供选择抗为协变量的孔隙度属性模拟,借助于三维可视化的、等概率的高分辨率实现技术,我们可以大致看到孔洞储集体的形态、分布、对于该工区来说,三维地震资料分辨率较高,对规模及连通性(图1).根据孔隙度发育情况,我们将孔洞储层已经有一定的反映(常表现为低频、不连续强储层分为孔洞欠发育(致密)、孔洞较发育(较致密)振幅反射).通过岩石物理分析又发现孔洞储层低速、和孔洞发育(较疏松)三种类型并分别设计了各自相低密,常规波阻抗反演能够刈其几何形态、空间接触应的弹性参数,同时以模拟出的孔洞形态约束弹性关系定量表征.因此本文将波阻抗数据体作为协变波正演模拟时孔洞体的空间分布量地震属性纵向等值法),采用确定性的协克里金2波动方程正演模拟原理插值算法,对孔洞储层的物性参数进行了三维建模反演所起到的作用,是通过归一化的测井曲线对碳酸盐岩岩溶风化壳孔洞型油气藏属于一种典原始地震数据进行校正,使数据在空间上得到了有型的、复杂的非均质范畴,可以视为由准均匀介质中效的平衡,从而使孔洞反映的更清楚;二是在地质建呈不规则分布的、大小和形状各异的低速体共同组模过程中通过宏观控制,充分利用空间变量的相关成的非层状储集体.在地震剖面上看到的储集体的性,克服低频模型的不足提高属性模拟的分辨率波应是这些低速体的散射(绕射)波.若利用常规波2【V动方程正演模拟方法所使用的均匀介质中的声波方N()程或弹性波方程,难以得到具有复杂非均匀性的孔cline洞型油气藏的地震波场响应2.因此,本文采用非均匀横向各向同性弹性介质中的弹性波波动方程进l()行正演模拟计算,取z轴为垂直对称轴,它可以表示为如下的一阶方程组:a0awta(λ+2naU)+AW(1)7(1(λ1+21)2+λax图1地质模型孔洞储集体俯视图Fig 1 Top view of the cavity reservoiμ(ain geological model其中:(U(x,z,t),W(x,z,t)是速度向量;B(x,z)是密度ρ=g(x,z)的倒数,或者叫疏度;r建模过程中最大的难点是求取准确的速度场,τ(x,z,t),za=x(x,z;t),n=rn(x,z,t)是应木文首先收集整理了工区内6口探井、评价井的钻力张量.A,P/和A1,p分别为水平和垂直方向上井分层数据及多种测井曲线(电阻率、声波时差、密的拉梅系数;为一新的弹性常数可见,在横向各度、自然伽玛等),对其进行了归一化和环境校正,并向同性弹性介质中,独立的弹性常数有五个,它们是制作了合成记录.通过与井旁地震道对比准确标定密度、在垂直方向上的纵、横波速度及纵、横波的各了前中生界侵蚀顶面(15)石炭系双峰灰岩顶面向异性系数,即:(T5。),中下奥陶统顶面(T7)、下奥陶统蓬莱坝组顶面(T7),它们都是区域性的波峰反射在此基础+2uL, Vsi-A pL上采用25m*25m测线密度对该区块6.25km2的3U地震数据体T5°、T5、T74、T78层位进行了精细√λ1+2,CSV闵小刚,等:三维孔洞储层建模其地震波场演模拟1291在具体的有限差分解法上,除了规则网格外,非均匀介质模型的弹性波动方程正演模拟特别是种较为先进的交错网格(图2)最早由 madariaga提当每一个波长中的网格点数多于10个时, Levander出, Virtex在模拟各向同性介质屮SH波和P(1988)2的结果显示,网格色散与网格各向异性均Sⅴ波时也使用了这种差分网格,其精度为o(△2十可忽略不计△x2),在不增加计算工作量和存储容间的前提下,假设U,W分别为介质在x,z两个方向上的速和常规差分网格相比局部精度提高了4倍,且收敛度分量的离散量,R,T,H分别为rxr=和τx的离速度也较快. Levander2又将这种差分网格的精度散量,Lo,M,L1,M1和M2分别为y,kM,A⊥?P1和提高到o(△t2+△x). Crase2则发展了精度可达任g的离散量在(1)式中,各导数项均用中心差分来意阶的高阶交错网格法,但其计算量和内存要求比代替,在如图2所示的一个交错网格中,U,B在节低阶有限差分法大幅度增加.本文使用的是 Virieux点1处计算;W,B在节点2处计算;R,T,M,L,(1986)1的交错网格差分公式,其差分精度为和M1,L1在节点3处计算;H,M2在节点4处计算(△2+△x2), Ikelle l t和 Yung$ K(199)21说这样()式离散为4明该算法可以糈确、稳定地应用于任总复杂变化的=U+B,(R年,-R…)十B,(I1+-H),wH,n-v++B++△(r}一rn)+By2(T+-+),ry=对++(n+2M4)+,△m+-RW(2)+T+,;+(L1+2M1)△tW+U2)+M2△t鲁←z以将震源函数直接赋在rx和n的节点上来模拟震源,即Soure,t)=R(t)t_(source_x, soure_x, t)=R(t.此外,在震源没有激发之前地下介质内部所质点都是静止的,包括质点振动速度为零和所受应力为零.因此,初始条件为图2一个交错树格Fig. 2 A st0,r(x,z,t)=0(t≤0)(3)对于自山表面边界条件,本文采用了模型空间其中,上标k为时间t的离散量,下标i,分别为x的上部加空气的条件,然后再采用吸收边界条件把和z的离散量.△,△x,△z分别为t,x,z的步长空气上边界的弹性波吸收掉,对于空气的下界面,则鉴于 Ricker子波对地震波的分辨率较其它子作正常的分界面来认识,从而获得和实际应用中波函数高,因此,震源选用 Ricker子波,其形式为所采用的地表放炮、地表接收达到一致的效果.R(t)=[1-2Lmf(t-to)] Jexp[-(rf(E to))2]有限差分法在求解波动方程时,会产生不期望式中f表示子波的主频,t为子波持续时间,t为f的数值频散或称网格频散,导致数值模拟结果分辨的函数,在模拟地下激发的地震波时,有限差分交错辛降低2所谓数值频散实质上是一种因离散化求网格中的正应力x和x=是在同一节点上赋值的,解波动方程而产生的伪波动,这种频散既不同于波而vr和vn在节点处的数值并没有参与计算,因此可动方程本身引起的频散,也不同于因波传播的速度1292地球物理学进展26卷频率和角度变化而引起的频散,它是有限差分方法果我们在这里仅分为三种类型:孔涧欠发育(致求解波动方程时所固有的本质特征,无法避免.为了密)、孔洞较发育(较致密)和孔洞发育〔较疏松).消除这种数值频散,前人进行了大量研究,他们的结论是基木一致的,即为了消除数值频散,在使用二阶表1地层框架内各层物性参数有限差分方法时,每个功率对应的波长至少必须使Table I The properties ot each layer用11个网格点,面四阶有限差分则可用二阶差分网in stratigraphie framewor格点数的一半.木文采用的稳定性条件,即计算稳定p(nu/s) v(m/s) (kg/m3)的离散参数区域为151:r4G界面2500三叠系)以Lmd2m≤1(2m-1)fT50界面下石炭系顶)~T46l730≤Ld2m≤T56界面(2m-1)!(巴楚组顶)~T5023102350其中,T74界囿(下奥陶系顶)~T56±8002470T78界面(蓬芠坝组颠)灬T7460002650界面以3702此外,在做波动方程的模型计算时,由于只能在对于试验工区的每条线,其长度均为1625m个有限区域进行,而弹性波在其计算边界上能量为了侏证该区域内均为满叠、孔洞的绕射波收敛以衰诚并不为零,从而产生很强的边界反射,这是模型及边界吸收较为干净,我们在模型的左边延长了计算时所不希望的,需要做人工吸收戌衰减处理,计1200m,右边延长了1575m(延长部分的地层接触算吸收边界的方法有许多种,一般情况下网格周围系并不代表真实情况),即模型总长度为1、4km,的耗散采用质点的速度和应力值乘上一个小于1深度范围为4000~6500m每条线均采用了同样的的因子来平滑的衰减;另种可能性是在网格周围观测系统,具体为:采用零相位对称雷克子波作为震使用低Q值来实现吸收作用,但是后者的吸收效果源(主频40玎z),单边放炮(共20炮,每炮128个检不如前者的吸收效果好,因此本文采用的是第一种波器接收)炮间距160m,检波器间距20m,8次叠方法,具体实珧时釆用了〔 eran等的吸收边界条加,最小偏移距0m,最大偏移距2540m,记录长件实现边界吸收1.6s,Δt=2ms,第一炮的坐标位置为(-1200,0)exp[-a2(I-i)2],1≤i≤1.基于差分稳定条件,取模型中最小介质速度2500m/s其中,I为给定的吸收边界带总节点数;i为吸收边为参考,得到的计算参数为:网格剖分间隔3m界内的节点编号;a为衰减系数,其值的选定与1的3m,时间延拓步长为0.27ms,每个波长(62.5m)大小密切相关,且对吸收效果的影响很大本文中Ⅰ内有20.8个网格.我们一共对33条线进行了正演取为40(即围绕计算区域,再向外设置宽度为40个模拟,图3展示了较为典型的 inline2585线(位于研网格的条形吸收区域)a=0.305/40,i取从0~40究工区的中心部位,地层接触关系以及孔洞体的分节点号.在条形吸收区域中的每个网格结点处,对全布相对比较复杂),从中可以大致看出二维正演模拟部的5个波场量(U,W,R,T,H),在每计算一个时的普遍情况与孔洞体波场响应特征的一般规律问步长后,都做少量的波场减表2展示∫该条线上各孔洞体的几何及物性参3模型计算数,其中④号属于欠发育(充填致密物)类型,①③⑥号属于较发育(充填较致密物)类型,②⑤号属于发在正演之前,我们统计了工区的速度、密度资育(充填较疏松物)类型.此外,建模过程中,我们还料,为了简化模型,并使得孔洞体的地震响应特征更考虑了线与线之间地层起伏渐变、孔洞大小渐变孔具有针对性,我们采用了背景为常速介质、蜜度参数洞物性参数渐变的过程,即所有建模因素都渐变由( arner公式计算的思路(表1).对于孔洞储集的而不是突变,最终保证了三维地震数据体的连体,根据钴井揭示和前面提到的孔隙度属性模拟结续性4期闵小刚,等;三维孔洞储层建模及其地震波场正演模拟1293表2各孔洞储集体的几何及其物性参数最大振幅,且绕射曲率与反射曲率相同,表明二者具Table 2 The geometry and propcrty parameters有不同的传播速度;每个绕射波可分为左右、上下f each cavity reservoir正、反向绕射分支,正向绕射分攴的相位与反射P孔润体尚度宽度vVP波相同,反向绕射分支的相位反转180°,与反射P(m)(m)(m/s)(m/s)(kg/m3)充填物屮心距界面(m)波的相反17396500029002503较致密105弹性波正演模拟生成的炮域合成记录被导人10113480027822470较疏松6FOCUS软件进行常规处理,包括速度谱拾取、动82784500029002500较致密85校、切除、增益、滤波、叠加、偏移和变面积、变密度显①575520030222530致密示等.由于在观测系统中只设计了8次覆盖,为了增⑤18115480027822470较疏松104加速度谱拾取精度,本文采用了由相邻的7个CDP2714850029002500较致密86道集混合构成一个超道集的办法,隔10个CDP拾图4是该模型在590ms时的波场快照,其波场取一个速度文件,并在拾取前先作常规NMO校正清晰,网格频散小,边界吸收较干净这表明,在求解切除,使得原始道集记录能量更强、信噪比史高二维弹性波动方程时,将差分解法和交错网格技术图5、图6分别是TK610井、TK623井所在位置处相结合,通过较好地使用吸收边界和稳定性条件可CDP道集记录及其速度谱,从图中可见各个反射界以显著削弱数值频散,有效地提扃计算精度.同时面的同相轴清晰可辨,对应的能量团集中,而在合成在保证一定的精度前提下,可以采用铰大的空间网记录上T7界面下孔洞所在位置处都有一明显的格间距,提高计算效率.从图巾还可看出,孔洞绕射同相轴,能量团也比较集中,由于TK610升比波和反射波在绕射点处相切,在切点处绕射波具有TK623井孔洞储集体更为发育(尽管二者振幅相1200-80004008001200160020002400280040004505500图3主测线2585地质模型Fig 3 Geclcgical model ol inline258-12004004008001200160045000.10.3图4主测线2585在590ms时波场快照Fig. Snapshot of wave field at 590ms in inlinc25851294地球物理学进展26卷Sg224-230CDP49 SE QNO250030003500400045000.240.60.60.80.8TE1.01. 01.21.2141.4图5TK610井所在位置处CDP道集记录及其速度谱ig. 5 The CDP gather and velocity spectrum at well TK610Sgl58-1640.2ONO250030003500400045000.0.40.60.60.8081.0:1.01212623(2565图6TK623井所在位置处CDP道集记录及共速度谱Fig 6 The CDP gather and velocity spectrum at well TK623当,在地质模型设计时均认为是充填较疏松物,但相消),使得T7界面断断续续,并在该界面下出现TK610井比TK623井在目的层段的厚度要大,横些“短反射”通过仔细分析,我们发现“短反射”中向展布范围也更宽测试产能更高),在合成记录上较强者出现的时间,与孔洞位置相对应.从该模型的孔洞对应的同相轴振幅更强、波形更连续,速度谱上偏移剖面上(图b)可以看到,所有的孔洞体均得到能量团也更强、更集中比铰好的偏移成像,并表现为负正负三个相位的图7是处理完后的叠加和偏移剖面.从叠加剖波形.但鉴于反射波地震勘探的纵向分辨率(大于面上(图a)可以比较清楚的看到孔洞体顶、底的两1/4波长),所有能检测出的孔洞或孔洞组合在叠加组强反射,但是二者之间出现具有绕射特征的弱波剖面上都叠合在T74界面下第一个波峰轴上,在偏代替了成层的背景,这些绕射波的相互下涉(相长、移剖而上都体现在T7界面下第一个黑椭圆体上,4期闵小刚,等:三维孔洞储层建模及其地震波场正演模拟1295601001401802202603003401001401802202603003400.00.00.20.2040.40.60.60.81014露9.926c+089.926e1081022e+091022e+09图?主测线2585对应的叠加剖而(a)和偏移剖面(b)Fig. 7 The stacking section (a) and migration section (b) of inline2585至」其下的“串珠”是孔洞的假象(孔洞组合与围岩(b)之间的多次波及绕射波经偏移归位后形成较强短反L2560L2580L2600射).由于T74界面反射波与沿纵横向有一定分布的孔洞(比较明显的是①、⑤号)的绕射波叠加,使得30001300040孔洞所在位置处T74界面反射波能量变弱,而孔洞底部与下覆围岩之间的正极性反射由于受T7界200600面反射波的负值性续至波叠加,也变得较弱.此外,B40080在构造高点上(④⑥号孔洞体所在位置,④号更为明显),由于孔洞引起的绕射与隆起引起的回转波的相6003600100互丁涉,T56和T7界面不连续,甚至在其间出现空白反射,而实际资料也有这种情况.这说明,对于塔3800800200和油田碳酸盐岩孔洞储集体这类特殊的油气储层来400040001400400说,在解释时遇到层位问断时,不能轻易地开断层,而应该综合考虑构造、孔洞绕射等地震波场特征.这图8联络测线2795实际剖面(a)和正演剖面(b)对比也是塔河油出勘探开发实践中发现“表层弱反射、内Fig 8 Comparison of the actual section (a) and幕强反射”地震特征对应有利储层的一个佐证forward modeling section (b)in crossline2795依据以上思路与工作流程,我们得到了33条沿主测线方向的二维偏移剖面,在并成三维体之前,为而正演模拟釆取的是8次叠加、道间距10m),正演了尽量消除线与线之问因地层起伏造成的不闭合,剖面较好的反映了实际情况.这不仅体现在层位的我们采取先把33条线的速度文件并成三维体,整体形态、分布比较相似(由于速度取了平均,各层的厚平滑后两用每条线对应的、平滑后的速度对其原始度不一致,但不影响我们的主要的,即对孔洞体地共中心点数据进行动校、叠加、偏移的办法,得到33震响应特征的分析),更重要的是,我们所设计的孔条新的二维偏移剖面,再并成一新的三维体,此外,泂体,其位置、形状规模、振幅强弱均与实际地震资由于正演模拟数据体线间距为50m,道间距为料具有相当好的对应关系,这表明我们在止演模拟10m,其空间采样率比实际资料低,本文编制了相和处理时的设计思想和参数选取原则是合理的,这应的算法在频率域对其进行插值,使线间距加密到结果也为我们进行后续工作提供了比较好的数据25m图8是联络测线2795在时间域的实际剖面源由于实质上是2.5维,不是基于面元的真三维,(a)和正演剖面(b)对比,排除二者在采集时的一些所以沿联络测线的剖面上同相轴有抖动现象,这是差异(如实际三维采集资料为24次叠加道间距25m,不可避免的)1296地球物理学进展26卷4结论与建议2]谢桂生,刘洪,赵连功,伪谱法地震波正演模拟的多线程并行计算[冂.地球物理学进展,2005,20(1);17~23.本文从三维角度,建立了与实际资料比较吻合Xie G S, Liu H, Zhao L G. Parallel Algorit hm based on the的孔洞储层模型,并进行了弹性波正演模拟,总结了multithread Technique for pseudospectal modeling of seismic地震响应规律,主要结论如下:wave[J]. Progress in Geophysics(in Chinese), 2005, 20(1)1)结合地震资料建立储层地质模型能够有效[3]刘财,张智,邯志刚,等.线性粘弹体屮地震渡场伪谱法模拟地降低储层模型的不确定性,提高建模精度.同时利技术[门].地球物理学进展。:005,20(3),640~644,用协克里金技术,用波阻抗反演的确定性信息约束Liu C, Zhang Z, Shao Z G, et aL. Pseudo-spectral forward储层的平面非均质性,可以实现孔隙度属性的确定modeling nf seismic wave in linear viscoelasic solid [J]P1性建模),2005,20(3):640~644.4」张智,刘财,邵志刚,伪谱法在常Q粘弹介质地震彼场模拟(2)在求解二维弹性波动方程时,将差分解法和中的应用效果[].地球物理学进展,2005,20(4):945交错网格技术相结合,通过较好地使用吸收边界和949,稳定性条件可以显著削弱数值频散,有效地提高计Zhang Z, Liu C, Shao G. The application of pseudo-spectral算精度.同时,在保证一定的精度前提下,可以采用forward modeling of seismic wave field in constant Q较大的空间网格间距,提高计算效率该方法具有广viscoelastic medium [J]. Progress in Geophysics, 2005,20(4)945~949泛的适用性5]盖良国,马在出,曹景忠,等.一阶弹性波方程交错网格高阶(3)孔洞储集体在偏移剖面上表现为负-正-负差分解法[冂].地球物理学报,200,43(3):411-~419三个相位的波形,但只能确定奥陶系风化面下第Dong LG, Maz T, Cao j Z, et al. A staggered-grid high个负相位是孔洞的发育位置,其下的“串珠”是孔洞order difference method of one-order elastic wave equation]的很象.风化面反射波与沿纵横向有一定分布的孔Chinese J. Geophys. (in Chinese),2000,43(3):411-419洞体的绕射波叠加,使得孔洞所在位置处风化面反[6]董艮国,马在田,曹景忠,一阶弹性波方程交错网格高阶差分解法稳定性研究[门].地球物理学报,200,43(6):856~射能量变弱,而孔洞底部与下覆闱岩之间的正极性反射由于受风化面透射波的负值性续至波叠加,也Dong L G, Ma Z T, Cao J Z. a study on stability of the变得较弱.该结论对于实际地震资料处理、解释以及staggcred-grid high-order difference method of first-order储层预测烃类检测具有普遍的指导意义elastic wave equation. Chinese J. Gcophys. in Chinese)2000。43(6);856~864本文不足之处主要有三点「7]萤良国.复杂地表条件下地震波传播数值模找1.勘探地球(1)在三维孔隙度建模时采用的是常规阻抗信物理进展,2005,28(3);187~194息(约束稀坑脉冲反演),其纵向分辨率不够(只能分Dong L G. Numerical simulation of seismic wave propagation辨1/4波长以上的孔洞储集体),在后续工作中将尝under complex near surface conditions [J]. Progress in试使用地质统计学反演的阻抗体来约束建模以大幅Exploration Geophysics(in Chinese), 2005, 28(3):187--194提高纵向分辨能力[8奚先,姚姚,二维随机介质及波动方程正演模拟[.石油地球物理劫探,2001,36(5):546-552(2)在弹性波正演模拟时采用的是2.5维思想XiX, Yao Y. 2D random media and wavc cquation forward口前正在研制全三维算法有望更逼真的还原孔洞储modeling [J]. Oil Geophysical Prospecting in Chinese集体的真实地下情况001,35(5);546~5523)考虑到缝的各向异性更为复杂,本文尚未涉9]奚先,姚姚,二维粘弹性随机介质中的波场特征分析[刀地及,对于碳酸盐岩中这类油气运移的重要通道,将在球物埋学进展,2004,19(3):608~615今后的工作中进一步研究Xi x, Yao Y. The analysis of the wave field characteristics in2-D viscoelastic random medium LJ. Progress in Geophysics参考文献( References):hinese),2004,19(3):608~[10]奚先,姚姚,二维横各向同性弹性随机介质中的波场特征1]刘文岭.大庆宋芳屯油田芳2区块地震与地质资料综合储层J.地球物理学进展,2004,19(4):924~932地质建模研究(博土论文儿D1.北京:中国地质大学,2002Xi x, Yao Y. The wave field characteristics of 2-DLiu W I. A Study on Reservoir Geological Modeling withclo].ESeismic and Well-log Data in Fang 2 Area of DaqingGeophysics(in Chinese), 2004,19(4):924-932ongfangtun Oil Field (doctor dissertation)(in Chi[111吴永国,贺振华,黄德济.串珠状溶涧模型介质波动方程正Beijing: CUG, 2002.演与偏移[.地球物理学进展,2008,23(2);539~5444期闵小刚,等:三维孔洞储层建模及其地震波场正演模拟1297Wu Y G, He Z H, Huang d J. Wave equation forward[19]肖玉茹,何峰煜,孙义梅,等,古洞穴型碳酸盐岩储层特征modeling and migration for heads-shaped corroded cave model研究一以塔河油田奥陶系古洞穴为例匚门。石油与天然气地EJ]. Progress in Geophysics(in Chinese), 2008, 23(2): 539质,200324(1):71~80.Xiao YR, He f Y, Sun Y M, et al. Reservoir charactetistics12]股文,印兴耀,吴国忧.高特度频率域弹性波方程有限差分of paleocave carhonates-a casc study of Ordovician paleocave方法及波场模拟[」.地球物理学报,2006,49(2):561in tahe oilfield, Tarim basin UJ]. Oil Gas Geology(inChinese),2003:24(1):71-80Yinw, YinXi,WuGC. The method of finite difference of[20]姚蟋,唐文榜.深层碳酸盐岩岩溶风化壳洞缝型油气藏可检high precision elastic wave equations in the frequcncy domain测性的理论研究[门.石油地球物理勘探,2003,38(6):623and wave-field simulation [J. Chinese J, Geophys.629Chinese),2006,49(2):561~568.Yao Y, Tang W B. Theoretical study of detectable cavern[13]马贵,土尚旭,宋建勇.频率域波动方程正演中的多网格Fractured reservoir in weathered Karst of dccp carbonatite迭代箅法[门].石油地球物理勘探,2010,45(1):15[J]. Oil Geophysical Prospecting(in Chinese), 2003,38(6):Ma ZG, Wang S X. Sun J Y. Multigrid iterative algorithm in623~629,domain wave equation forward modeling [J]. Oil [21] Levander A R. Fourth-order finite difference P-SvGeophysical Prospecting(in Chinese ) 2010, 45(1): 1-5seismograms []. Geophysics, 1988, 53(11): 25-36.[14]张金海,王卫民,赵连锋,等.傅里叶有限差分法三维波动[22] Crase e. Iligh- order( space and timc) finite-difference方程正演模拟[.地球物理学报,2007,50(6):1854A, In: 60th SEG Annual1862C].1990:987~991.Zhang j H, Wang W M, Zhao L F, et aL. Modeling 3-D [23] IkelleL T, Yung SK, Daube F. 2-D random media withscalar waves using the Fourier finite-difference method.ellipsoidal autocorrelation function [J]. Geophysics, 199350(6):1854[24]奚先.随机介质模型的构造及其波场模拟(博土论文)[D][15] Qin Z, Lt武汉:中国地质mproved NPML absorbing boundary condition in elastic waveXix. Construction and scismic wave field modeling ofmodeling [J]. Applied Geophysics, 2009. 6(2): 113-121random medium model doctor dissertation ) in Chinese)[16][D].Wwave equation [J. Geophysics, 1986, 51(1): 54-61[25]吴国忱,王华忠.波场模拟中的数值频散分析与校正策略[17] Virieux J. P-Sv wave propagation in heterogeneous mediaLJ.地球物理学进展,2005,20(1):58-65velocity-strcss finite-difference methud LJ]. GeophysicsWu GC, Wang H Z. Analysis of numerical dispersion in1986,V51;889~901.wave-field simulation [J]. Progress in GreaphysiEs ( in18] Igel H, Riollet B. Mora P. Accuracy of staggered 3-D finiteChinese),2005,20(1):58~65difference grids for anisotropie wave propagation [J]. 62th [26] Cerjan C, Kosloff D, Kosloff R, et al. A nonreflectingAnn, Internat, Mtg, Soc. ExpL. Geophys, Expboundary condition for discrete acoustic and elastic- wav1992,1244~1246.equation []. Gcophysics, 1985, 50(4): 705-708.

本代码用于仿真车辆在两条单向三车道道路上移动的情形，考虑到了加速度，超车，背景流量等

滑模控制在工程中的应用，学习和了解滑模的实际应用，学以致用

这是一本非常经典的控制方面的书，值得控制专业的人看看。3.2闭环系绕的方块图(58)333开环传递函数和前向传熄函数59)3,3,4闭环传递函数…""r5935扰动作用下的闭环系统……(59)336画方块图的步骤……………………………………………………………(60)3.3.7方块图的简化……………61)34.状态空间棋型………633.4.1现代控制理论P曾中■■■曾晋■冒■t曾4冒留634.2现代控制理论与传统控制理论的比较…………·■會t■曾■1曾■曾··曾?(6334.3状态"■量■■…………………(63)3.4.4状态变量加■鲁■p■b…(63)3.4.5状态向量…………""…"”TT………………(64)34.6状态空间∵……………………………………"…r………………………(64)34.7状态空间方程…………………………………………………………(64)34.8传递函数与状态空间方程之间的关系……""………{67)34.9传递矩阵………………………………………(69)35动态系统的状态空间表达式………〔69)3.5.1线性徽分方程作用函数中不包含导数项的n阶系统的状态空间表达式693.5,2线性徽分方程作用两数中包含导数项的n阶桑统的状态空阔表达式………………(71)3.6机减系统…幽幽■■……………………………(73)36.1质量…………………………………(74)36.2力中k■bd·tb■■。■hd昌h日b日■…(75)36.3说明………………………………………………………(75)3.64忧械系统(75)3.7电气系统(79)3.7.1LRG电路……………………(80)3.72复阻抗1『甲■■(803.73状态空闯表示…(81)3.74串联元件的传递函数…"…"……………"…"(823.75无负载效应串联元件的传递函数■4■p■p电动■■p●申自■电伽d■……………(83)3.8液位系统………………………(84)3.8.1液位系统的液阻和液容…,…,……,(8438.2液位系统B··d看+!F4日卓中中··■甲宁·日·早……………………(86)38.3相丘有影响的液位系统即B*自日·日·日■·■●·■■■。■。pdψbψhpdb晶■■■dbh■备■“■■■备(87)3.9热力系统………:(88)3.9.1热阻和热容893.92热力系统………………………………………………”…………(89)3.10非线性数学模型的线性化电t中自自■P■■唱■。■91)3.10.1非线性数学模型的线性近似………………………………………(91)Ⅷ3.10.2液压何服系统的线性化·.··.s:中命··。·咖■■要申中中血血··血唱自咖92)例题和解答………………………………………………(95)习题………………………………………………………………………(117)第4章瞬态响应分析4.1引言…■■·鲁■·■■……………………(124)4.1.]典型试验信号.(124)4.【.2瞬态响庇和稳态响应(1244.13绝对稳定性、相对稳定性和稳态误差h■卜■■嚼P自↓h■■4■q『■申曾q鲁『吾P1【會鲁P中·PP曾早目·日·擊中124)本章要点(125)4.2一阶系统(125)4.2.1阶系统的单位阶跃响应·吾B备日吾号目备B■日吾日中吾日『目·即■电中會■·日电中中P■·中目p目■■(125)4.2,2一阶系统的单位斜坡啊应…■■■命■■■b■自命■■■……………………(l27)4.23一阶系统的单位脉冲响应“(128)4.24线性定常系统的重要特性……………………………………………(128)4.3二阶系统画■■■mpq甲b■4.3.1直流伺服马达…………………………(」31)4.3.2同服系统■■■↓■■■L■昌L■↓L■↓L■↓■■↓晶↓■郾↓■↓昌■唱■■↓■■看…(131)4,3,3负载对伺服马达动特性的影响■"44L.··°·4··p.4日·■(134)4.3.4二阶系统的阶跃响应……(134)4,3,5瞬态响应指标的定义(1374.3.6关于瞬态响应指标蚋几点说明……4(139)4.3.7阶系统及其瞬态响应指标:hdh==如如d吾■■亠如139)4.3.8带速度反愤的何服系统…↓亠↓L↓■↓■■LL■甲喜自■卩■司甲4●db■4,3,9二阶系统的酥冲响应4.4用 MATLAB进行瞬态响应分析…………………………………………(1474.1引言…………………………………………………………………(147)44.2线性系统的 MATLAB表示■即■…………………(147)4.4.3传递函数系统单位阶跃响应的求法…………(148)4.4.4在图形屏幕上书写文本…………………………咖■●■·■■●■■■■bdp■■■如k晶■晶司■■■冒h152)4.4.5脉冲响应■幽■·●■■口血■甲■甲……·(15244.6求脉冲响应然另一种方法■b山…………(155)4.4.7斜坡酝效■昌L■■■■■■■■■L■■■■冒?上……(157)4.4在状态空间中定义的系统的单位斜玻响应………………………………………(158)4.49对初始条件的响应(传递函数方法)■■■■■■■山■晶■b■晋b山4.4.10对初始条件的响应(状态空间方法情况1)…………(162)4.4.1l对初始条件的响应(状态空间方法,情况2)…………………………(16345用 MATLAB解题举例_警冒晋日P冒中吾冒q晋個…(165)4.5.1机械振动系统(165)4,5,2计算机仿真(连续时同方法}…t………………………………(168)M4.53计算机仿真《离散时间方法)·…例题和解答(73}习题………(192}第5章控制系统的甚本控制作用和南应L■………………(196)5.1引言(196)5.2基本控制作用p司■(196)52.11业控制器的分类“…+4L昌■■4b亠』L■↓“■■L■L■■■品昌↓■k↓昌■↓■昌522自动控制器,执行机构和传感器(测量元件)………■·■■■自自·血■(197)523自操作控制器……………195.2.4衩态或继电器型控跏作用(198)5.2,5比例控制作用……………………………(199)5.26只分控制作用……………(2005.27比例积分控制作用5,28比例徵分控制作用品日4日吾4品日4………(201)529比例-加积分加微分控制作用……………………!…,"!(2015.21传越器(测量元件〕紂系统性能的影啊……………………(202)53积分和微分控制作用对系统性能的影响…………(203)5.3.1积分控制作用(2035.3.2液位控制系统的积分控制(204)53,3对转矩扰动的响应(比钢控制)……n534利用MAAB求响应………………………………(2陌)535对转矩扰动的响应(比侧加积分控制)…………………………*(207)5.36分控制作用2095.37带惯性负裁系统的比例控制………………0095.38具有惯性负载系统的比例--徽分控制■■↓■晋h■■↓■昏b■h■■+■晋■■晋■昏·画;b(209)539二阶系统的比例加微分控制………………………………………………………(210)5.4高阶系统(211)4.1套阶系统的瞬态响应分析………………………………………………………(2ll54.2闭坏主导极点■中▲■■■↓·;■↓■k昌■■·■』■■■■q昌昌■■即■↓■口是即自↓■……………………(213)54.3复平面内的稳定性分析…(213)55劳斯稳定判据……_p·…………………(215)5.5.1劳斯稳定判据简介…(215)5.5.2特殊慒祝………………………s…………………………(217)5.5.3相对稳定性分析…………(219)5.5.4劳斯稳定判据在控制系统分析中的应用……219)56气动控制器…即即上命聊(225.6.1气动系统和液压系统之间的比较(2205,6、2气动乐统·+4·(21}5,6,3压力系统的气阻和气容……………………s……(22l1)64压归系统(222)5.6.5气劲喷嘴-挡板放大器Lb可·+Lqb命(223)566气矿接续骺………(224)56.7气动比例控制器(力-距离型(225)568气动比例控制器(力平衡型)………甲中幽自司b电血■即…………(228)569气动执行……會■『冒冒口■斷『T■■P■骨■冒會TtP■…………………(229)5.6.10获得徽分控制作用时基本原理………………………………(230)56.11获得气动比例抓积分挽制作用的方法(232)56.12获得气动比倒如积分如微分控制作用的方法上■雪會曾·『中平曾會會會(234)5、7液压控制器(2355.7.1液压系统……(235)5.7.2液压系统的优缺点……………”…(235)5,73说明……-………………………(25)5.74液压积分控制番…(235)575液压比例控制器……………(237)5.76级冲器■L■■■I(238)57.7获得液压比倒加积分控制作用的方法………………(2405.78获得液压比例加黴分控制作用的方法(24158电子控制器(242)5.8.1运算放大器吩■ P Pk J即pbmm……(242)5.8.2反相放大器……………+……………………(243)58,3非反相放大器…………………………………(243)5.8.4求传递函嫩的阻抗法……2465.8.5利用运算放大器构成的超前或滞后网络………■鼻■■■昏昌■■■L■晷■(246)5.86采用运算放大器的PD挽制器……h口口■b口■■…(24859正弦响应中的相位超前和相位滞后……………(2495.10单位反馈控制系统中的稳态误差……(253)5.10.1控制系统的分类………………………………………………………(2535.10.2稳态误差……………………………(253)513静态位置误差常数K…254)5.104静态速度误差常数E……………………(255)5.10.5静态加速度误差常数K……………………(256)5,10,5小结…………(257)5.10.7开环控制系统与闭环控制系统中稳态误差的比较………………………(258)例题和解答(259)习题(285)第6章根轨迹分析294)6.Ⅰ引言中4L昌qb中↓昌↓↓d■↓■昌■■「昌品■■冒(294)6.1.]根教迹法294)6.1.2章要点■■冒冒↓矗■■■晋P■1P■·甲■↓昏冒甲甲11晋曾甲【■■·■■■『平甲口暑甲■冒(295)62根轨迹图■t血●…(295)6.2.1辐角和幅值系统q昌qp中▲L晶mdb4L(295}6,2.2示例………………………………297)6.3根轨迹作图的一般规则■·■■冒血■上■■"■口■1■『曾曾··P平俨中(305)6.3,1作根轨迹图的一般规死·…(306)6.32关于根轨迹图的说明30963,36(g)的极点与Hs)的等点的抵消-……(3106.3,4典型的极零点分布及其相应的根轨迹…■■■↓■■■■L■(311)6.3.5小结……(31264用 MATLAB作根凯迹图………(32)6.5特殊情况··b44如吾b4+吾=如6.5.1变量参数不以乘法因子形式出现时的作根护迹的方法(3226.5.2正反馈系统的板轨迹(326)6.6控制系统的根轨迹分析…30)66.】根轨迹与定常增益轨迹的亚交性……………………(306.6.2条件稳定系统晶■↓↓晶■■噌晷郾↓■■司鲁■即■■即330)663非最小相位系统…(332)具有传递延迟的系统的根轨迹∴……(332)68根轨迹族曲线……即司●中自q中■血↓自咖d·口司D即自■……(336)例题和解答…■q聊d■●■wφh哂看■■■■甲ψ画讠■晋■■■h■晶h■↓■冒h(340)题………-……4(370)第7章控制系统设计的根轨迹法…………………………………………(375)7.l引TP『4中自曾中l“■■■中■■…………:…41(375)7.1.1性能指标………(375)7.2系统的校正■凸■■■■■看■■曾■番↓■警b■h■冒鲁语■■昏鲁十画■P■375)7,1,3串联校正和反馈或并联)校正……………4(375)7.1.4校正装置………………………………………………"……(3767.15设计步骤……·(3777.16本章要点…………………………………………………(377)了2初步设计研究…如吾b4.+4备·"日b4日………………(37772.1控制系统设计的根轨迹法………………·(37872.2增加概点的影响…日日即·中(378)7.2.3增加零点的影响……………………………………"…………(378)7,3超前校正………………………………………………………(379)7.3.1超前网络(3797,32基于根靴迹法的超前校正技术…P自中■■■中(3807.3.3说明………■■↓b福h■■h■昏↓■4h■■4■冒PP■(3857.3.4校正与未校正系统阶联响应的比较……………(3857,4滞后校止………………………(387}74.1●采用运算放大器的电子滞后饺正装置…………,……………(387)74.2应用根轨迹法进行滞后校止……↓■ IJd·d●◆■q·甲·甲■晶■吾■■『·(387)7.4.3用根轨迹法进行带后校正设计的龙骤·即4·bdb■■■+日P中中自·b■……………(38875滞后超前校正……………T冒『·冒4(3957.5.1利用运算放大器构成的电子滞后-超前校正装置7.52基于根轨迹法的淠后超前校下方云……………………………(396例题和解答(405)习题…(433)第8章频率响应分析…看■■■■…(438)8.1小言女_备吾“···*".·""""··P+“8.1.1系统对正弦输入信号的稳态输出备如A如4吾日.甲甲.·昌………∴…(438)81.2用图形表示频率响应性■d■■↓■……,(439)8.13本章要点…………………(439)8.2伯德图甲山日古B4+日··日·日·自咱D中‘■日十F(40)8.2.1伯德图或对数坐标图↓■}■画▲■■晶■■■『P■『"4··▲■备■■4■昏(443)822G(})H(c)的基本因子……………………41(4408.2.3增益K……"………………………(440)积分和徽分因子(〕441)82.5一阶因子(1+7)…一↓■■↓■■∵……(442826二阶因子[L+2(/mn)+/ex)2]+1………………………*………827谐振频率u和谐振峰值Mb■b■督■■冒P冒q甲…·(448)82.8绘制伯德图的一般步骤士晋晋晋!p……………·(449)82.9最小相位泵统和非最小相位系统rp*4+·〓如山血■日■P自噜中·■■·日■唱···(451)8.2.10传递延迟………〔452)8.2.11系统类型与对数幅值曲线之间的关系…(4548,2,12静态位骱误兼常数的确定-……,,……(454)82.13静态速度误差常数的确定…………………4(455)8.2.14静态加速度误差常数的确定……………………………………(456)83用 MATLAR作伯德图……↓4晶昌■■备■司晷阜……,……(457)8.3.1在一定的類率点上增益变成无穷大时对伯德的影响…(464)8.3,2求状态空闸中的系统的伯德图…………"………¨……………(466)84极坐标图v·自司·■■·■■·■日■唱唱申申卓鲁日tb■……(468)8,4.1积分和微分因子(如)3…(468)842…阶因子(1+户T)↓pmu·卓■→·日·q4如■十■吾■■■■唱P…,(46943阶因子[1+26(/an)+(/mn)2]平(470)844传递延迟472)8.45极坐标图的一般形决(474)85用MAB作奈魁斯特图………:…(476)8.5,1注底…479)8.5.2定义在状杰空间的系统的奈魁斯特图画法……………………………(481)8.6对数幅-相图■P·昏P昏P曾■P■晋■■『q■晋■■■晋■昌晶■d……………………s(484)8.7奈魁斯特稳定判据會中号甲看P■(486)8.7.1预备知识………………………(487)87.2映射定理…甲P甲申……(490)8.73映射定理在闭环系统稳定性分析中的应用……,…:4(490)8.74奈魁斯特稳定判据早P日音卡吾日14‘44b4b自(491)875关于奈魁斯特稳定判据的几点说明(492)8.76G(s)H(s)含有位于轴上射极点和(或〕零点的殊情祝……………………………(493)88稳定性分析…『■『餐……………1(4958.8.1条件定系统■■晶■口■■■……(499)8.8.2多凹路系统↑·4中可中中!广曾?冒曾■『■?『■十D■自■血日自咖司即●p●中甲电甲m看·口电(499s.8.3应用于逆极坐标图上的奈魁斯特稳定判据(501)884利用改变的奈魁斯特轨迹分析相对稳定性■_會『會會■■■■會■■個■會■(504)8.9相对稳定性…:·+··=·+“+·4+“+日·中曾十◆·■日『『■■■…………………(506)891遁过保角变换进行相对稳定性分析………(5068.92相位裕量和增益裕量(588.93关于相位裕量和增益裕量的几点说明…………………………………………(51089.4谐振峰值幗值M和谐振蜂值版率n…512)895杯崔二阶系统中阶联瞬态啊应与频率啊应之间的关系5138.96一般系统中的阶跃瞬态响与频率响应之间的关系…………(515)897截止频率和带宽…………………………………………………………(516898剪切率…………………,……(517)8.10闭环频率响应……………………………(518)810.1单位反馈系统的闭坏频率响应4『F:-T自+卡…(518)8.10.2等幅值轨迹(M遴)………(519)8.103等相角轨迹(N)………………,……∵,……………(520)g.104尼枓尔斯图……………s522)810.5非单位反馈系统的闭不频率响应最↓bL(525)8.10.6增益的整…(525)8.11传递函数的实验确定法…………………………………-…………(5288.11.1正弦信号产生器528)8112由伯德图求最小相位传递函教……………(529)3113非最小相位传递函数……………………………*……………(530)8.114关于实验确定传递区数的儿点说明例题和解答■血■■b■盘血■口■■自■■■■■■■■■_■■口■■盘血■鲁■1●和中血D命申…………(533习题(56第章控制系统设计的频率垧应法(571)引肓咖●?控制系统设计的频率响应法…………………(5719.1.2从开不頻率响应可以获得的信息…唱卓幽■“●■◆■■如····Paqq↓·甲571)9.1.3对严环频率响应的要求………44572)91.4超前滞后和滞后超前校正的基本特性…………∵(53)915本幸要点…573)92超前校正……………………………………"………………(573)92!超前校正装置的特性…………(573)922基于频率响应法的超前校正………………………(574)9.3滞后校币…■山d■…………581)9.3,1滞后校正装量的特性w■!冒■■昏『■晋■冒+『晶………"……(581)93,2幕于频卒响应法的滞后校正…………………………………(582)9.3.3关于滞居校正的一些说明588)9.4滞后-超前校正…………(589)94.1溢后-超前校正装置的特性·■曾曾4·冒骨T中·「冒曾雪『『會自口t曾會4會■■『中曾平·目自曾會日·日幽(589)9.4.2基于频率响应法的滞后超前校正……………·(5919.5结论5959.5.1超前、滞后和滞后超前校正的比较(595)95.2图形对比……………………………………………(595)9.53反馈校正………(5969.5.4不希望极点的抵消……………"(596)9.5.5不希望的共轭复数极点的抵消……(597)56结束语…鲁自日·b中日即………………………………………(598)例题和解答■唱↓昌郾昌■■昌昌■↓■■■〖』』晶■晷晶L■↓口■■↓■■■■即中↓(598)习题622)第10章PID控制与鲁控制………(625)10.!引言………(62510,2PI控制器的调节律………………………………………………(625)10.2.1控制对象的PD控制…………(625)10.2.2用来调整P控制器的齐格勒尼柯尔斯法则(62610.3.3第一种方法…………………………………………………(626)10.24第二种方珐……(628)102.5说明……"…10.3PI控制方案的变形(634)10.3.1PID控制63510.3.2IP担齣…………………(63610.3.3从IPD控制方案推广到带状态反馈的积分控制方案367)10.3.4二由度PD托制↓4吾··B日k日(6380.4二自由度痉制………,………(638)

【实例简介】利用python载入邻接矩阵绘制网络图，基于python语言的特点，对邻接矩阵加以处理后再进行应用，即将邻接矩阵去除第一列（节点序号列），复杂网络的基本拓扑结构可以用图论的方法表示成G =(V，E)，V中元素称为节点或顶点，E中元素称为边。在图论框架下, 可用不同的全局参量来表示复杂网络的特征, 最基本特征包括度分布 P(k)和聚类系数C。

自己总结的左神以及牛客网的算法pdf资料，算法是真的需要大家好好敲代码才行，加油，兄die们，共勉

多目标perota优化MATLAB算法，运用算法，通过MATLAB进行自动优化

zynq-7000 pcb design guide.pdf