第一篇：講稿版利用Matlab模擬點電荷電場的分布

利用MATLAB模擬點電荷電勢的分布

一、目的

1.熟悉單個點電荷及一對點電荷的電勢分布情況； 2.學會使用MATLAB進行數值計算，并繪出相應的圖形；

二、原理

根據庫侖定律：在真空中，兩個靜止點電荷之間的作用力與這兩個電荷的電量乘積成正比，與它們之間距離的平方成反比，作用力的方向在兩個電荷的連線上，?Q1Q2?R（式1）兩電荷同號為斥力，異號為吸引力，它們之間的力F滿足：F?k 2R?kQ?EE由電場強度的定義可知：?2R（式2）

R對于點電荷，根據場論基礎中的定義，有勢場E的勢函數為U?kQ（式3）R在MATLAB中，由以上公式算出各點的電勢U，可以用MATLAB自帶的庫函數繪出相應的電勢分布情況。

三、MATLAB基本語法

(一)標識符與數

標識符是標志變量名、常量名、函數名和文件名的字符串的總稱。(二)矩陣及其元素的賦值

賦值就是把數賦予代表常量或變量的標識符。MATLAB中的變量或常量都代表矩陣，標量應看作1×1價的矩陣。賦值語句的一般形式為 變量=表達式（或數）列如，輸入語句 a=[1 2 3;4 5 6;7 8 9] 則顯示結果為 a=1 2 3 5 6 8 9 輸入 x=[1 2 3 4 5 6 7 8 9] 結果為

x= 1 2 3 4 5 6 7 8 9 可以看出，矩陣的值放在方括號中，同一行中各元素之間以逗號或空格分開，不同行的元素以分號隔開。語句的結尾可用回車或逗號“，”，此時會立即顯示運算結果；如果不希望顯示結果，就以分號“；”結尾再回車，此時運算仍然執行，只是不作顯示。

變量的元素用圓括號“（）”中的數字（也稱為下標）來注明，一維矩陣（也稱數組）中的元素用一個下標表示，二維矩陣可有兩個下標數，以逗號分開。在MATLAB中可以單獨給元素賦值，例如，a(2,3)=6,x(2)=2等。(三)元素群運算

把n×m矩陣中的每個元素當作對象，成群地執行某種運算，稱為元素群運算。元素群運算能大大簡化編程，提高運算的效率，這是MATLAB優于其它許多語言的一個特色。

1、數組及其賦值

數組通常是指單行或單列的矩陣，一個N階數組就是1×N或N×1階矩陣。N階數組可以表示N維向量。

在求某些函數值或曲線時，常常要設定自變量的一系列值，例如，設間隔n在x軸上從-3到3之間，每隔1取一個點，共7個點，這是1×7階的數組。如果逐點給它賦值，將非常麻煩。MATLAB提供了兩種給等間隔數組賦值的簡易方法。（1）用兩個冒號組成等增量語句，其格式為x=[初值：增量：終值]。例如，鍵入

x=[-3：1：3] 得x=-3-2-1 0 1 2 3

當然增量為1時，這個增量值是可以略去的。

（2）linspace函數表述等距離分割，其格式為x=linspace(初值，終值，點數)。例如鍵入 x=linspace(-3,3,7)得x=-3-2-1 0 1 2 3

在x軸上-3和3實際上是一個點，所以這個命令是把x軸分為7份。第三個變元也可以不寫，此時取默認值100.2、元素群的四則運算和冪次運算 元素群運算也就是單個元素之間的運算。為了與矩陣作為整體的運算符號相區別，要在運算符“*”、“/”、“”、“^”前加一個點符號“.”,以表示進行元素群運算。矩陣的加減法本來就是對元素進行的，故不再有元素群運算符。參與元素群運算的兩個矩陣必須是同價的（標量除外）。

下列的例子可以說明利用元素群運算的優越性。例如，要求列出一個三角函數表。這在MATLAB中只要一下兩個語句

鍵入 x=[0：0.1：pi/4]’;[x, sin(x), cos(x), tan(x)] 第一條語句把數組x賦值，經轉置后成為一個列向量。因為sin, cos,tan函數都對元素群有效，得出的都是同階的列向量。第二條語句把4個列向量組成一個矩陣，進行顯示。

得 0

0

1.0000

0

0.1000

0.0998

0.9950

0.1003

0.2000

0.1987

0.9801

0.2027

0.3000

0.2955

0.9553

0.3093

0.4000

0.3894

0.9211

0.4228

0.5000

0.4794

0.8776

0.5463

0.6000

0.5646

0.8253

0.6841

0.7000

0.6442

0.7648

0.8432 第一列是x，以下各列依次是sin(x), cos(x), tan(x)。for語句

for語句的結構形式為

for k=初值：增量：終值

語句組A，end 即它把語句組A反復執行N次。在每次執行時程序中k值不同。可以算出循環次數為N=1+（終值—初值）/增量 用for 語句求三角函數表的程序為 for x=0：0.1：pi/4 disp([x, sin(x), cos(x), tan(x)]);end 所得結果將和上個例題中的答案相同。可以看出，MATLAB的元素群運算功能與一個for循環相當，由于它不需要每次檢驗表達式，運算速度比for語句快得多。for 語句可以嵌套使用。

四、MATLAB現用函數介紹

MATLAB語言的難點是函數較多，僅基本部分就有700多個，其中常用的近200個，要盡量多記少查，以提高編程效率，而且這是終生受益的。

1、線型分隔函數linspace（-xm,xm，n）：在-xm與xm之間均分地產生n個點值，形成1×n元向量。其中-xm是初值，xm是終值，n是點數。

2、兩個變量的標量指令 [X,Y]=meshgrid(x,y)：將向量x,y變換為數組X,Y，這樣就可以將兩個一維向量生成兩個二維矩陣。其中，數組X的各行是向量x的拷貝，總行數為y向量的元素個數；數組Y的各列是向量y的拷貝，總列數為x向量的元素個數。例如x=-3:3;y=1:5;[X,Y]=meshgrid(x,y)X=-3-2-1 0 1 2 3-3-2-1 0 1 2 3-3-2-1 0 1 2 3-3-2-1 0 1 2 3-3-2-1 0 1 2 3 Y= 1 1 1 1 1 1 1 2 2 2 2 2 2 2 3 3 3 3 3 3 3 4 4 4 4 4 4 4 5 5 5 5 5 5 5 Z=(X+Y).^2 Z= 4 1 0 1 4 9 16 1 0 1 4 9 16 25 0 1 4 9 16 25 36 1 4 9 16 25 36 49 4 9 16 25 36 49 64

3、二維空間中繪制等高線函數 contour(X,Y,U,u)：就是挑選出電勢相等的點，并在向量u指定的值上繪制等高線。注意u取值要恰當。

四、編程

1.畫單個點電荷的平面電場線與等勢線 clear %清除變量

E0=8.85e-12;%真空電介質常數 C0=1/4/pi/E0;%歸并常數 q=1.6*10^(-19);%元電荷電量 xm=2.5;%橫坐標范圍 ym=2;%縱坐標范圍

x=linspace(-xm,xm);%橫坐標向量 y=linspace(-ym,ym);%縱坐標向量 [X,Y]=meshgrid(x,y);%設置坐標網點 R=sqrt(X.^2+Y.^2);%點電荷到場點的距離 U=C0*q./R;%計算電勢

u=1e-9:0.5e-9:5e-9;%等電勢的電勢向量 figure %創建圖形窗口 contour(X,Y,U,u)%畫等高線 hold on %保持圖形

plot(0,0,'o','markersize',12)%畫點電荷 axis equal tight %使坐標刻度相等

title('單個點電荷的平面電場線與等勢線','fontsize',20)%顯示標題 xlabel('r','fontsize',16)%顯示橫坐標 ylabel('E(U)','fontsize',16)%顯示縱坐標

2.畫一對點電荷對的電場線和等勢線 clear %清除變量

E0=8.85e-12;%真空電介質常數 C0=1/4/pi/E0;%歸并常數 q=1.6*10^(-19);%元電荷電量 a=1;%電量比 xm=2.5;%橫坐標范圍 ym=2;%縱坐標范圍 x=linspace(-xm,xm);%橫坐標向量 y=linspace(-ym,ym);%縱坐標向量 [X,Y]=meshgrid(x,y);%設置坐標網點

R1=sqrt((X+1).^2+Y.^2);%第一個正電荷到場點的距離 R2=sqrt((X-1).^2+Y.^2);%第二個正電荷到場點的距離 U=C0*q./R1+C0*q./R2;%計算電勢 u=1e-9:0.5e-9:5e-9;%等電勢的電勢向量 figure %創建圖形窗口 contour(X,Y,U,u)畫等勢線 grid on %加網格

legend(num2str(u'))%圖例 hold on %保持圖像

plot(-1,0,'o','markersize',12)%畫第一個正電荷 plot(1,0,'o','markersize',12)%畫第二個正電荷 axis equal tight %使坐標刻度相等

title('等量同號點電荷的電場線和等勢線','fontsize',20)%顯示標題 xlabel('r','fontsize',16)%顯示橫坐標 ylabel('E(U)','fontsize',16)%顯示縱坐標

txt=['μ?oé±è:itQrm_2/itQrm_1=' num2str(a)];%電荷比文本 text(-xm,-ym-0.3,txt,'fontsize',16)%顯示電荷比

下列的例子可以說明利用元素群運算的優越性。例如，要求列出一個三角函數表。這在MATLAB中只要一下兩個語句 鍵入 x=[0：0.1：pi/4]’;[x, sin(x), cos(x), tan(x)] 第一條語句把數組x賦值，經轉置后成為一個列向量。因為sin, cos,tan函數都對元素群有效，得出的都是同階的列向量。第二條語句把4個列向量組成一個矩陣，進行顯示。得 0

0

1.0000

0

0.1000

0.0998

0.9950

0.1003

0.2000

0.1987

0.9801

0.2027

0.3000

0.2955

0.9553

0.3093

0.4000

0.3894

0.9211

0.4228

0.5000

0.4794

0.8776

0.5463

0.6000

0.5646

0.8253

0.6841

0.7000

0.6442

0.7648

0.8432 第一列是x，以下各列依次是sin(x), cos(x), tan(x)。for語句

for語句的結構形式為

for k=初值：增量：終值

語句組A，end 即它把語句組A反復執行N次。在每次執行時程序中k值不同。可以算出循環次數為N=1+（終值—初值）/增量 用for 語句求三角函數表的程序為 for x=0：0.1：pi/4 disp([x, sin(x), cos(x), tan(x)]);end 所得結果將和上個例題中的答案相同。可以看出，MATLAB的元素群運算功能與一個for循環相當，由于它不需要每次檢驗表達式，運算速度比for語句快得多。for 語句可以嵌套使用。

第二篇：電磁場與電磁波點電荷模擬實驗報告

重慶大學 電磁場與電磁波課程實踐報告 題 目：

點電荷電場模擬實驗 日 期：

2013 年 12 月 _J_ 日 N=28

《電磁場與電磁波》課程實踐 點電荷電場模擬實驗 1.實驗背景 電磁場與電磁波課程內容理論性強，概念抽象，較難理解。在電磁場教學中，各種點電荷的電場線成平面分布，等勢面通常用等勢線來表示。

MATLAB 是一種 廣泛應用于工程、科研等計算和數值分析領域的高級計算機語言，以矩陣作為數 據操作的基本單位，提供十分豐富的數值計算函數、符號計算功能和強大的繪圖 能力。為了更好地理解電場強度的概念，更直觀更形象地理解電力線和等勢線的 物理意義，本實驗將應用 MATLA 對點電荷的電場線和等勢線進行模擬實驗。

2.實驗目的 應用 MATLA 模擬點電荷的電場線和等勢線 3.實驗原理 根據電磁場理論，若電荷在空間激發的電勢分布為 V ,則電場強度等于電勢 梯度的負值，即：

r E V 真空中若以無窮遠為電勢零點，則在兩個點電荷的電場中，空間的電勢分布為：

q i q 2

V

y V 2

— 4 本實驗中，為便于數值計算，電勢可取為 V 蟲 R 4.實驗內容 應用 MATLA 計算并繪出以下電場線和等勢線，其中 q i 位于(-1,0,0), q 2 位于(1,0,0), n 為個人在班級里的序號：

(1)電偶極子的電場線和等勢線(等量異號點電荷對 q 2 ：

q i = 1，q 2 為負電荷)； ⑵

兩個不等量異號電荷的電場線和等勢線(q 2 ：

q 1 = 1 + n/2，q 2 為負電荷)；(3)兩個等量同號電荷的電場線和等勢線； 0 R 1 4

0 R 2 R 2

⑷

兩個不等量同號電荷的電場線和等勢線(q 2 ：

q 1 = 1 + n/2)；(5)三個電荷，q 1、q 2 為⑴中的電偶極子，q 3 為位于(0,0,0)的單位正電荷。、n=28(1)電偶極子的電場線和等勢線(等量異號 點電荷對 q 2 :q 1 = 1，q 2 為負電荷)； 程序 1 ：

clear all q=1;xm=2.5;ym=2;x=linspace(-xm,xm);y=linspace(-ym,ym);[X,Y]=meshgrid(x,y);R1=sqrt((X+1).A2+Y.A2);R2=sqrt((X-1).A2+Y.A2);U=1./R1-q./R2;u=-4:0.5:4;figure contour(X,Y,U,u, “--”);hold on plot(-1,0, “o” , “MarkerSize” ,12);plot(1,0, “o” , “MarkerSize” ,12);[Ex,Ey]=gradient(-U,x(2)-x(1),y(2)-y(1));dth1=11;th1=(dth1:dth1:360-dth1)*pi/180;r0=0.1;x1=r0*cos(th1)-1;y1=r0*sin(th1);streamline(X,Y,Ex,Ey,x1,y1);dth2=11;th2=(dth2:dth2:360-dth2)*pi/180;x2=r0*cos(th2)+1;y2=r0*sin(th2);streamline(X,Y,-Ex,-Ey,x2,y2);axis equal tight title(“ 卩 ？？？？？

Xo^ ? 卩 ？ 3???o ???>ntsize” ,12)

點偶極子的電場線和等勢線

(2)兩個不等量異號電荷的電場線和等勢線(q 2 ：

q i = 1 + n/2 , q 2 為負電荷)； 程序 2 : clear all q=15;xm=2.5;ym=2;x=li nspace(-xm,xm);y=li nspace(_ym,ym);[X,Y]=meshgrid(x,y);R1=sqrt((X+1).A 2+Y.A 2);R2=sqrt((X-1).A 2+Y.A 2);U=1./R1-q./R2;u=-4:0.5:4;figure con tour(X ,Y, U,u, “--”);hold on plot(-1,0, “o”，‘MarkerSize“ ,12);plot(1,0, ”o“ , ”MarkerSize“ ,12);[Ex,Ey]=gradie nt(-U,x(2)-x(1),y(2)-y(1));dth1=11;th1=(dth1:dth1:360-dth1)*pi/180;r0=0.1;x1=rO*cos(th1)-1;y1=r0*si n(th1);streamli ne(X,Y,Ex,Ey,x1,y1);dth2=11;

th2=(dth2:dth2:360-dth2)*pi/180;x2=r0*cos(th2)+1;y2=r0*si n(th2);streamli ne(X,Y,-Ex,-Ey,x2,y2);axis equal tight title(” 卩 ？？？？？

Xo^ ? 卩 ？ 3???o ???>ntsize“ ,12)點偶極子的電場線和等勢■線

-2.6-2-1.5-1-0.6 0 0.5 1 1.6 2

(3)兩個等量同號電荷的電場線和等勢線;程序 3 ：

clear all q=-1;xm=2.5;ym=2;x=li nspace(-xm,xm);y=li nspace(_ym,ym);[X,Y]=meshgrid(x,y);R1=sqrt((X+1).A 2+Y.A 2);R2=sqrt((X-1).A 2+Y.A 2);U=1./R1-q./R2;u=-4:0.5:4;figure con tour(X ,Y, U,u, ”--“);hold on plot(-1,0, ”o“，‘MarkerSize” ,12);

plot(1,0, “o” , “MarkerSize” ,12);[Ex,Ey]=gradie nt(-U,x(2)-x(1),y(2)-y(1));dth1=11;th1=(dth1:dth1:360-dth1)*pi/180;r0=0.1;x1=r0*cos(th1)-1;y1=r0*si n(th1);streamli ne(X,Y,Ex,Ey,x1,y1);dth2=11;th2=(dth2:dth2:360-dth2)*pi/180;x2=r0*cos(th2)+1;y2=r0*si n(th2);streamli ne(X,Y,Ex,Ey,x2,y2);axis equal tight title(“ 卩 ？？？？？

Xo^ ? 卩 ？ 3???o ???ntsize” ,12)點偶極子的電場線和等勢線

-2-15 「 1 0 0.5 1 1.,5 2

(4)兩個不等量同號電荷的電場線和等勢線(q 2 ：

q 1 = 1 + n/2);程序 4 : clear all q=-15;xm=2.5;ym=2;x=li nspace(-xm,xm);y=li nspace(_ym,ym);[X,Y]=meshgrid(x,y);R1=sqrt((X+1).A2+Y.A2);R2=sqrt((X-1).A 2+Y.A 2);U=1./R1-q./R2;

u=-4:0.5:4;figure con tour(X ,Y, U,u, “--”);hold on plot(-1,O, “o”，‘MarkerSize“ ,12);plot(1,O, ”o“ , ”MarkerSize“ ,12);[Ex,Ey]=gradie nt(-U,x(2)-x(1),y(2)-y(1));dth1=11;th1=(dth1:dth1:360-dth1)*pi/180;r0=0.1;x1=r0*cos(th1)-1;y1=r0*si n(th1);streamli ne(X,Y,Ex,Ey,x1,y1);dth2=11;th2=(dth2:dth2:360-dth2)*pi/180;x2=r0*cos(th2)+1;y2=r0*si n(th2);streamli ne(X,Y,Ex,Ey,x2,y2);axis equal tight title(” 卩 ？？？？？

Xo^ ? 卩 ？ 3???o ???ntsize“ ,12)點偶極子的電場線和等勢線

(5)三個電荷，q i、q 2 為(1)中的電偶極子，q 3 為位于(0,0,0)的單位正電荷 程序 5: clear all q=1;q3=-1;xm=2.5;ym=2;x=linspace(-xm,xm);

y=linspace(-ym,ym);[X,Y]=meshgrid(x,y);R1=sqrt((X+1).A 2+Y.A 2);R2=sqrt((X-1).A 2+Y.A 2);R3=sqrt(X.A2+Y.A2);U=1./R1-q./R2-q3./R3;u=-4:0.5:4;figure contour(X,Y,U,u, ”--“);hold on plot(-1,0, ”o“ , ”MarkerSize“ ,12);plot(1,0, ”o“ , ”MarkerSize“ ,12);[Ex,Ey]=gradient(-U,x(2)-x(1),y(2)-y(1));dth1=11;th1=(dth1:dth1:360-dth1)*pi/180;r0=0.1;x1=r0*cos(th1)-1;y1=r0*sin(th1);streamline(X,Y,Ex,Ey,x1,y1);dth2=11;th2=(dth2:dth2:360-dth2)*pi/180;x2=r0*cos(th2)+1;y2=r0*sin(th2);streamline(X,Y,-Ex,-Ey,x2,y2);dth3=11;th3=(dth3:dth3:360-dth3)*pi/180;x3=r0*cos(th3);y3=r0*sin(th3);streamline(X,Y,Ex,Ey,x3,y3);axis equal tight title(” 卩 ？？？？？

Xo^ ? 卩 ？ 3???o ???>ntsize" ,12)

點偶極子的電場線和等勢線

-1.6-1 白 0 05 1 15 2

從實驗過程中學習到的東西: 1.靈活學習，大膽求證，當不清楚 E1,E2 前面符號的正負時，隨便假設一 個，再根據電荷的正負關系，看得到的圖形是否正確，若不正確則再修 改符號 2.注意 q 的正負與兩電荷是否異號有關，異號與同號 q 的正負不同 3.學習初步使用 matlab 軟件，為以后的學習打好基礎 4.更加深入地了解電荷的電場線與等勢線

第三篇：利用多點地質統計學方法模擬巖相分布論文

傳統的地質統計學是利用變差函數描述地質變量的相關性和變異性，通過建立在某個方向上兩點之間的地質變量的變化關系來描述空間的變化特性.但是，建立在兩點統計關系上的變差函數本身在描述儲層非均質性上有很大的不足，它反映的僅僅是空間兩點之間的相關性，不能充分描述復雜幾何形狀砂體如河道砂體和沖積扇砂體空間的連續性和變異性.當井資料較少時，用于計算實驗變差函數的點對很少，它也就不能正確反應空間兩點之間的相關性.建立在兩點基礎上的變差函數在儲層地質建模中存在一定的不足，而多點地質統計學[14是建立在多個點的相關關系上，所以它在解決描述空間變量的連續性和變異性方面得到越來越廣泛的應用.斯坦福大學的Journel教授曾指出多點地質統計學是今后地質統計學發展的方向，它的優勢已越來越顯著.1多點地質統計學的原理

在闡述多點地質統計學之前，首先回顧一下變差函數的地質統計學方法是如何模擬儲層巖相分布的，以序貫指示模擬算法為例進行說明.該方法的基本原理簡述如下：

假設在模擬區域有k種巖相S1，S2…?》，對于模擬目標區域內的每一相，定義指示變量：

對于任一待模擬點，其出現第k種相的概為:P(Ik=1lz(u)SaVa)，a為待估點所包括的條件區域，利用兩點地質統計學方法計算該概率是采用克里格方法：

其中，&為克里格方法確定的權系數，它通過求解由變差函數或協方差函數建立的克里格方程組來確定.多點地質統計學與兩點地質統計學的主要區別在于上面的概率的確定方法不同，它首先引入一訓練圖像，通過在訓練圖像中尋找與待估點內條件數據分布完全相同的事件的個數來確定概率分布，因此它可以反映出多個位置的聯合變異性.例如，計算圖1()中u點的概率時，相應的條件數據場為da={Z(ul)Z(u2)，Z(U3)，Z(U4)}，其基本方法是首先要在訓練圖像(b)中尋找與圖(a)中數據分布完全相同的事件的個數，即要在訓練圖像中找出與圖（a)幾何完全相同的區域，同時在該區域中相同的位置處z(U1)，Z(u2)，z(u3)，z(u4)的值完全相同.在訓練圖像中一共找到4個既能滿足條件數據u1,u2，u3，u4數值，同時又能滿足它們分空間幾何形狀的事件，在這4個事件中，3個事件的u點的值為0，只有1個事件中u點值為1，因此u點巖相為1的條件概率為P{u=1Idn}=1/4,而P{u=0Idn}=3/4,這樣便可求出了u點的條件概率.因此u點巖相為1的條件概率為P{u=1Idn}=1/4,而P{u=0Idn}=3/4,這樣便可求出了u點的條件概率.上述方法不僅考慮了區域內條件數據的值而且也考慮了條件數據的幾何形狀.而兩點地質統計學只是依靠Z(u1)，Z(u2)，Z(u3)，Z(u4)的值及各點與u點距離通過求解克里格方程組來確定u點的概率，并沒有考慮dn的幾何形狀和各條件數據的配位關系.基于上述原理，SebastienStrebelle提出了snes-im模擬算法121，利用該算法可以快速、靈活地模擬巖相分布.該方法的具體步驟為：

(1)利用非條件模擬建立三維訓練圖像；

(2)定義通過所有待估結點的隨機路徑；

(3)對隨機路徑中的任意待估點/(=1，2，…，1):①定義查找范圍內的條件數據;②保留鄰區的數據點;③在訓練圖像中尋找與該區域內條件數據完全相同的事件，計算該點巖相的分布概率.④由MontoCarlo法得到位置處的一個模擬值;⑤將模擬結果歸入條件指示數據集中.(4)重復上一步模擬，直到所有的點全被模擬.訓練圖像既可以通過非條件模擬求出，也可以通過該地區的地質露頭資料分析得出.對訓練圖像的條件非模擬可以選擇非條件的布爾模擬方法，其方法和原理參見文獻.2實例分析

對于開發中后期的砂巖油藏儲層參數模擬采用兩階段模擬方法可以較為準確地反映儲層的非均質性，而“兩步建模”的第一步就是要建立儲層結構或流動單元模型，模擬沉積體在空間排列的復雜性;利用多點統計學模擬方法可以較好地完成砂體骨架模擬.模擬區域選擇我國東部某砂巖油藏第15小層，在該層一共有64口井，測井資料解釋結果表明有26口井鉆遇砂體，另外38口井鉆遇泥巖，砂體比例為40%.對巖相進行編碼，砂巖為1，泥巖為0，圖2為該層井位分布圖.采用上述方法模擬砂體的分布.首先建立訓練圖像，運用布爾模擬方法，把砂體比例40%輸入，為保證訓練圖像數據充足，網格劃分為250X250X1，一共由62500個模擬數據組成，布爾模擬結果見圖3.把條件數據和布爾模擬生成的訓練圖像，輸入到snesim模擬算法中進行模擬.根據該區域的特點，橢圓最大搜索半徑選為300m,搜索半徑內最多的條件數據設為30，搜索主方向選擇物源方向5°，得到該層的砂泥巖分布（圖4).從模擬結果看出，它很好地滿足了條件數據，即在各井點處的模擬結果與數據相一致，這表明該方法為條件模擬.同時，模擬的砂體展布方向和趨勢與依靠地質經驗手工繪制的砂體展布圖（圖5)比較吻合，在模擬的左下角與左上角砂體的展布與手工勾繪的幾乎完全一致，但該方法在局部區域表現出砂體展布的非均質性和不確定性，與手工勾畫砂體展布的平滑而唯一的表現是具有一定差別的，它充分體現了砂體局部的變異性和非均質性.3結論

(1)多點地質統計學是今后地質統計學發展的主要方向，它可以聯合反映空間多個位置點的幾何形狀和相互配位關系;在模擬具有復雜形狀地質體分布時，它比兩點地質統計學方法具有更大的優勢.(2)利用snesim模擬算法可以快速靈活地進行多點地質統計模擬，模擬的巖相展布圖具有一定的真實性，它為儲層參數的兩階段模擬奠定了基礎.

第四篇：中國廣東蜜源植物的分布及利用情況

中國廣東蜜源植物的分布及利用情況

廣東省蜜源植物分為野生和人工種植兩種。野生蜜源主要分布在山區，多達100多種，最主要的有鴨腳木、柃屬植物和山烏相等；人工種植植物的蜜源植物主要為果樹，面積超過100萬公頃，主要為荔枝、龍眼，還有柑橘、李、柚子、柿子、香蕉、菠蘿和雜果。此外，在全省還有大面積水稻，在雷州半島有大面積的核樹林以及一定面積的橡膠林和椰子樹，在粵北山區有少量的油菜。在全省各地還有一些人工種植的綠化花木等。

四季蜜源占有量分為，春季蜜源占30％，夏季蜜源占20％，秋季蜜源占10％以下，冬季蜜源占40％以上。一年只有三個生產季節，即春季的荔枝、龍眼，夏季的山烏柏和烏桕，冬季的鴨腳木和柃屬植物。秋季僅有一些輔助蜜源。

荔枝主要分布在北回歸線以南地區，以廣州地區，粵西的高州、化州、電白，珠江三角洲的東莞、深圳和粵東的普寧、惠來、惠東最多。龍眼以高州、化州和花都種植面積最大。柑梅全省都有種植，其中以潮汕平原、新會、四會、化州等老種植區，以及粵北、梅縣、惠陽等地種植面積最大。油菜在粵北的英德、韶關、翁源、乳源、曲江還有少量種植。香蕉主要分布在珠江三角洲，以東莞、順德及高州等地種植面積最多。桉樹全省各地都有種植，以雷州半島等地的面積最大，多為人工種植的防護林。

近年來，在韶關、梅州、清遠、湛江等地為主大量種植的尾葉桉，也稱藍桉，速生桉，是從澳洲及東南亞引進后經選育的新品系，主要用于造紙工業，其木材是造高級紙的用料，一般種植5～6年就可成材。種植4年后開花，每年9月中旬開花，花期一個月，與細葉桉花期相仿，尾葉核的蜂蜜，其特征性味道較其他核樹蜜差別不大，泌粉量不多。

在山區和丘陵地區，有大面積的次生林，生長著大量的山烏桕、野桂花和鴨腳木等主要蜜源植物。以廣州郊縣、惠州市、河源市、肇慶市、清遠市和英德市等地區分布最多。

廣東省有蜂群40多萬群，其中有30多萬群是在山區定地飼養，或小轉地飼養的中蜂；西方蜜蜂不足10萬群，主要為意大利蜂，多數流動放蜂。由于活框飼養技術的成功推廣，使中蜂成為廣東省的一種重要的授粉蜂種，主要用于對荔枝、龍眼等果樹的授粉，并取得明顯的增產效果。

全省養蜂數量和蜜源分布是不相稱的，一些地方由于交通不方便，蜂場無法到達，自然資源得不到充分利用。大面積的柑橘、水稻蜜源，由于開花期間施灑農藥，影響了蜂群的采集活動。人工種植的蜜源植物面積逐年增多，但山區的野生蜜源植物（如鴨腳木等），隨著山區資源的開發（鴨腳木可用來栽培食用菌）而遭到濫砍濫伐。近年來這種現象雖然引起政府有關部門的重視，采取措施加以限制，但要恢復到原來的水平，尚需一定時間。

每年都有2～5萬群的外地蜂群轉地來放蜂。原本是本省中蜂主要產品的冬蜜，部分地區出現大量西蜂進場采收、中蜂退居“二線”的局面，例如，在從化北部、佛岡、清遠、龍門等地鴨腳木花期，都有大量外來和本省的西蜂。本文摘自39蜂療網！

第五篇：《土地資源的利用和分布》教學設計

《土地資源的利用和分布》教學設計

教材分析

本節課是八年級地理上冊第三章“中國的自然資源”第二節“土地資源”,其教學重點是讓學生了解我國土地資源的概況和耕地、林地、草場、沙漠的分布特點,理解我國土地的基本國策。教材在地形、氣候之后安排土地資源,主要是考慮到地形、氣候是促進和影響土地資源形成與利用的主要自然地理環境要素。而土地資源的利用受人類長期的生產活動與社會經濟條件的影響,是人類賴以生存的最基本資源之一,是人類經濟活動(農業、工業、交通運輸等)的物質基礎。因此本節教學內容有明顯的承上啟下、前后關聯的特征。學情分析

初二學生思維活躍、感情豐富、求知欲強，在接受知識上帶有濃烈的感情色彩，樂于接受有趣的感情內容，所以在處理教材內容上，讓學生實際操作聯系實際生活，展開探究討論，使其在活動中發現問題、提出問題、分析并解決問題。

設計理念

本節課標要求“了解土地資源的利用類型和理解我國土地資源分布差異的原因”。讓學生充分結合學生日常生活中的所見所聞,積極主動地學習對生活有用的地理,實現學習方式的轉變,本節課就是以學生熟悉的被稱為“萬物之本”的土地為題,讓學生從自己的生活經歷和體驗中了解土地的用途、分布及利用中存在的問題、保護的措施等,以喚起學生節約和保護耕地的意識和行為。在學習探究的過程中擴展知識,培養能力,形成正確的人地觀,樹立可持續發展意識。

教學目標

知識與能力 1.使學生初步了解土地資源的利用類型和我國土地資源的基本特點。2.初步分析并理解我國土地資源分布差異的原因。

過程與方法 1.通過分析、理解我國土地資源分布差異的原因，提高學生分析問題、解決問題的能力。

情感態度與價值觀 通過對土地資源的利用類型、分布以及人類所面臨的世界性的土地問題，使學生真正認識到土地資源的重大意義，必須珍惜并加以合理利用，增強學生的資源觀。

教學重點、難點

重點：我國土地資源的分布及基本特點。難點：土地資源利用類型分布的原因

教學設計

[創設情景，導入新課] 教師：自然界中，以她最為寬廣深厚，她承載著自然界中的萬事萬物。一直以來，默默地以她仁慈寬厚的胸懷滋養著萬物生靈。她是山川之根、萬物之本、人類的衣食來源。這是什么？

學生：是土地。

總結：土地，它給人類提供了山川、萬物、衣食、財富；過去的一切得益于它，它更影響著人類的未來。這就是土地，這就是土地資源。

課堂活動： 一：知土地 教師：出示山地林區、草原牧區、平原耕作區、城市的景觀圖，討論幾幅景觀圖中土地自然環境和條件的差異以及利用方式的不同。

學生：山區──林地；高原──草地；平原──耕地：城市──建筑用地。

教師：從圖中可以看出，不同的土地類型，由于自然環境和自然條件的差異，人們對土地的利用方式和途徑就不同，根據土地的用途及土地利用的狀況，我們把土地資源分為耕地、林地、草地和建筑用地等，耕地、林地、草地為農業用地，建設用地為非農業用地。我國除了這些可以利用的土地外，還有很多難利用的土地，如荒漠、石山、永久冰雪地等。

[承接]我們知道了不同的土地類型，但你們知道祖國的土地資源有何特點嗎？

二：知特點 教師： 出示圖表、資料(部分國家人均耕地的比較圖以及農業用地總量和人均占有量數據)，同學們能得出什么樣的結論呢？

學生：由于我國國土遼闊，土地資源總量豐富，而且土地利用類型齊全。但我國人均土地資源占有量少。

問：我國土地資源的利用類型多種多樣，它為農業生產提供了哪些有利條件？

學生總結：有利條件：我國國土遼闊，土地資源總量豐富，而且土地利用類型齊全，這為我國因地制宜全面發展農、林、牧、副、漁業提供了有利條件。例如山區地形崎嶇，不宜耕作，但森林茂密，適宜發展林業；平原地勢平坦，土地肥沃，降水充足，適宜種植業發展。高原地處內陸，降水較少，草原廣闊，是畜牧業基地。

出示“我國土地利用類型的構成圖”。從這些資料中我們又能得出什么信息？

學生：耕地12.68%,最少；林地31.86%，較少；而難利用的土地15.72%多 學生：我國后備土地資源不足。

教師：我們十分需要的耕地、林地比重小，我們不需要的難利用的土地比重大，這樣的土地構成結構合理嗎？顯而易見是不合理的。我國耕地比重小，今后能不能得到大幅度的提高呢？

出示閱讀材料，學生朗讀。

教師：閱讀材料進一步說明了我國后備土地資源不足。大家可想而知，人與耕地的矛盾是很突出的，這些都是土地資源的劣勢。

教師總結：土地資源的構成特點可概括為“兩少、一大、一不足、一突出”，即：耕地、林地少，難利用土地多，后備土地資源不足，人與耕地的矛盾突出

承接：既然同學們了解了我國土地資源的利用情況，那么耕地、林地、草地又是如何分布的呢？

三：懂分布 出示“我國土地類型分布圖” 教師：：同學們，請你們讀圖，思考：我國的耕地、林地、草地和難以利用的土地分別是怎樣分布的。

師引生總結：我國東部地區以耕地和林地為主，西部地區主要是草地和難以利用的土地。這說明我國土地資源地區分布不均勻。

過渡：為什么不同地區的土地利用類型不相同？這與它們的自然環境有關。下面我們以小組合作的形式探討我國土地資源的分布與自然環境的關系。

學生合作探究，師生共同總結，完成板書

教師總結：不同地區，因為自然環境不同，所以土地的利用類型是不相同的。

教學反思

本節課所學內容是人教版“第三章中國的自然資源資源”中的第二節 “土地資源”。通過學習讓學生了解我國土地資源的利用類型和基本特點以及目前我國土地資源利用中存在的諸多問題，從而認識到土地資源的重要性，樹立合理用地、節約用地、保護耕地的資源觀，培養保護資源的意識。

在教學中，我采用了講解、探究、歸納的方法，達到了教學目的，重點突出。在解決難點時，設計了活動的方式，學生探究、交流，概括出我國土地資源的現狀，提高了學生分析、概括的能力。在教學設計上能夠突出學生的主體作用，學生通過探究活動認識了地理規律和原理，學會了分析解析地理現象。通過這一節課的教學，學生不僅學到了地理學科的結構知識，而且在一定程度上掌握了探究地理學習的思維方式。

存在的問題：

1、要找準切入。以點帶面，應選擇與本節課的內容聯系密切、真實有效、現實而熟悉的事物，2、教學中對學生預想不到位，在組織活動過程中分寸把握的不是很好，詳略處理的不是很好，3、本節課雖然選擇了大量的圖片和文字資料，對學生的學習起到了一定的指導作用，但在利用所選擇的材料利用的有效性還是欠缺。

4、教學經驗的積累和駕馭課堂能力還需要加強，雖然有一些客觀的原因，但作為一名合格的教師，應充分考慮到在課堂上有可能發生的問題包括知識和學生活動中的偶發事件。