Estimation of Soil Penetration Resistance Using Fractal Parameters of Particle and Aggregate Size Distributions

Document Type : Research Paper

Authors

Abstract

Penetration resistance (PR) is a highly variabe soil property and there is no accurate and representative PR data in most databases. On the other hand, direct measurement of the PR is costly and expensive. The objective of this study was prediction of PR from easily available soil properties using pedotransfer functions (PTFs). Fifty three soil samples were taken from dominant soil series of Hamadan province. Particle, macro and micro aggregate size distributions, water content and bulk density were measured in the laboratory. PR and soil structural qualitative parameters were measured in the field. In the first regression model, bulk density and water content as conventional predictors of PR were introduced to the model. Then, soil structural qualitative parameters (aggregates shape, type and size), primary particle, micro and macro aggregate size distribution fractal parameters were introduced to the models to predict PR. Finally the improvement of PR  prediction using various developed PTFs were evaluated by the statistical indices as R2, R2adj, root mean square error, Akaike  information criterion and relative improvement. The results showed that using fractal parameters as predictors along with the conventional predictors significantly improved the accuracy and reliability of the models. Using soil structural qualitative parameters along with fractal parameters of particle size distribution may significantly improve the PR predictions.

Keywords


برآورد مقاومت فروروی با بهره‌گیری از پارامتر‌های فراکتالی توزیع اندازه ذرات و خاکدانه‌ها

 

محمد جُرِّه، حسین بیات2*، علی اکبر صفری سنجانی3و ناصر دواتگر4

 

تاریخ دریافت: 03/11/90  تاریخ پذیرش: 06/08/91

1-دانشجوی سابق کارشناسی ارشد، گروه خاکشناسی دانشگاه بوعلی سینا همدان

2-استادیار گروه خاکشناسی دانشگاه بوعلی سینا همدان

3-استاد گروه خاکشناسی دانشگاه بوعلی سینا همدان

4-استادیار گروه خاکشناسی مؤسسه تحقیقات برنج کشور

* مسئول مکاتبه: h.bayat@basu.ac.ir: Email

 

چکیده

مقاومت فروروی (PR) خاک دارای تغییر‌پذیری بالایی است و در بسیاری از پایگاه‌های اطلاعاتی آگاهی چندانی درباره آن وجود ندارد. از سوی دیگر اندازه­گیری این پارامتر پرهزینه نیز می­باشد. بنابراین هدف این پژوهش برآورد PR از دیگر پارامتر­های زود یافت خاک با بهره­گیری از توابع انتقالی (PTFs) بود. تعداد 53 نمونه خاک از خاک­های شاخص بر پایه سری خاک­های غالب از استان همدان نمونه­برداری شد. توزیع اندازه ذرات (کوچک‌تر از 2 میلی‌متر)، توزیع اندازه خاکدانه­های درشت (25/0 تا 8 میلی‌متر) و ریز (کوچک‌تر از 2 میلی‌متر)، رطوبت خاک در حین اندازه‌گیری PR و جرم مخصوص ظاهری در آزمایشگاه اندازه­گیری شد. PR و پارامتر­های کیفی ساختمان خاک (شکل، نوع و اندازه خاکدانه) در محل نمونه برداری در دو لایه رویین (10 تا 35 سانتی‌متری) و زیرین (20 تا 45 سانتی‌متری) خاک اندازه­گیری شد. ابتدا جرم مخصوص ظاهری و رطوبت خاک که پارامتر­های مرسوم برآورد کننده PR می­باشند وارد مدل رگرسیونی شدند. سپس پارامتر­های کیفی ساختمان خاک (شکل، نوع و اندازه خاکدانه) و همچنین پارامتر‌های فراکتالی توزیع اندازه ذرات و خاکدانه‌های درشت و ریز نیز به پارامتر‌های پیشین افزوده شدند. در پایان میزان بهبود و برتری هر یک از مدل‌ها بر پایه آماره‌های ضریب تبیین،  ضریب تبیین تعدیل شده، مجذور میانگین مربعات خطا، معیار اطلاعات آکایک و بهبود نسبی ارزیابی شد. نتایج نشان داد که افزودن پارامتر‌های فراکتالی، به ورودی‌های مرسوم برآورد کننده PR باعث بهبود معنی‌دار معیار‌های ارزیابی صحت مدل‌ها شد. همچنین بکارگیری پارامتر­های کیفی ساختمان خاک به همراه پارامتر‌های فراکتالی توزیع اندازه ذرات باعث بهبود معنی‌دار پیش‌بینی PR شد.

 

واژهای کلیدی: بعد فراکتالی، توابع انتقالی،  مقاومت فروروی

 

 

 

 

Estimation of Soil Penetration Resistance Using Fractal Parameters of Particle and Aggregate Size Distributions

M Jorreh,1 H Bayat,2* A A Safari Sinejani3 and N Davatgar4

 

Received: 23 January 2012    Accepted: 26 November 2012

1-Former M.Sc. student of Soil Sci., Dept. of Soil Sci., Faculty of Agric., Bu Ali Sina Univ., Hamadan, Iran.

2-Assist. Prof., Dept. of Soil Sci., Faculty of Agric., Bu Ali Sina Univ., Hamadan, Iran.

3-Prof., Dept. of Soil Sci., Faculty of Agric., Bu Ali Sina Univ., Hamadan, Iran.

4-Assist. Prof., Dept. of Soil Sci., Rice Research Institute of Iran, Rasht.

*Corresponding Author Email: h.bayat@basu.ac.ir

 

Abstract

Penetration resistance (PR) is a highly variabe soil property and there is no accurate and representative PR data in most databases. On the other hand, direct measurement of the PR is costly and expensive. The objective of this study was prediction of PR from easily available soil properties using pedotransfer functions (PTFs). Fifty three soil samples were taken from dominant soil series of Hamadan province. Particle, macro and micro aggregate size distributions, water content and bulk density were measured in the laboratory. PR and soil structural qualitative parameters were measured in the field. In the first regression model, bulk density and water content as conventional predictors of PR were introduced to the model. Then, soil structural qualitative parameters (aggregates shape, type and size), primary particle, micro and macro aggregate size distribution fractal parameters were introduced to the models to predict PR. Finally the improvement of PR  prediction using various developed PTFs were evaluated by the statistical indices as R2, R2adj, root mean square error, Akaike  information criterion and relative improvement. The results showed that using fractal parameters as predictors along with the conventional predictors significantly improved the accuracy and reliability of the models. Using soil structural qualitative parameters along with fractal parameters of particle size distribution may significantly improve the PR predictions.

 

Key words: Fractal dimension, Pedotransfer functions, Penetration resistance

 


مقدمه

در خاک‌های سخت رشد ریشه و سبز شدن نهال با دشواری­های بسیاری روبرو می‌باشد. این خاک­ها هزینه­ها و انرژی کشت و کار را افزایش می‌دهند و در رژیم­های رطوبتی مناسب برای کشاورزی، این خاک­ها اعمال محدودیت می‌کنند (کمپل 1975). مقاومت فروروی ([1]PR) پارامتر[2] سودمندی در شناسایی فشردگی لایه‌های خاک و میزان نرم شدگی بین نمونه‌های خاک است. همچنین در شناسایی تغییرات مقاومت به فروروی در خاک با گذشت زمان و ارزیابی تأثیر مقاومت خاک بر رشد ریشه کاربرد دارد (اسالیوان و بال 1982). مقاومت به فروروی بستگی به نوع خاک، توزیع اندازه و شکل ذرات خاک، کانی شناسی رس، اندازه اکسید‌های بی شکل (آمورف)، درصد ماده آلی و شیمی محلول خاک دارد (جرارد 1965). به دلیل تغییرپذیری مکانی بالای PR، به دست آوردن مقادیر دقیق آن در کشتزارها دشوار است (وز و همکاران 2001). از اینرو پژوهشگران از توابع انتقالی (PTFs)[3] برای برآورد آن بهره‌ می‌برند. توابع انتقالی برای توصیف معادله‌هایی بکار می‌رود که وابستگی ویژگی‌های دیریافت خاک (مانند گنجایش تبادل کاتیونی و هدایت هیدرولیکی خاک) با متغیر‌های زودیافت خاک (مانند بافت، ساختمان، جرم مخصوص ظاهری و میزان مواد آلی) را نشان دهد (وستن و همکاران 2001). گزینش متغیرهای شایسته و بایسته به عنوان ورودی توابع انتقالی، در توانایی پیش‌بینی مدل از اهمیت ویژه‌ای برخوردار است (مسکینی 1389). نمس و همکاران (2003) دریافتند که اهمیت داده‌هایی که به عنوان ورودی توابع انتقالی بکار می‌رود بیشتر از چگونگی تهیه این توابع است. تاکنون از پارامتر‌های گوناگونی برای برآورد مقاومت فروروی در توابع انتقالی بهره‌گیری شده است که عبارتند از جرم مخصوص ظاهری (BD)[4] (پیدجن و ساون 1977 و هندرسون و همکاران 1988)، میزان رطوبت جرمی یا (GWC[5]) (فایر و داماتا 1994 و گرانوالد و همکاران a 2001)، بافت خاک (کوروپ و همکاران 1994 و پاپالا و همکاران 1995)، مکش ماتریک خاک (باسچر 1990)، عمق خاک (گرانوالد و همکاران b2001)، ماده آلی (تو و کی 2005)، تخلخل کل (TP[6]) و اشباع نسبی (θv/TP) (بیات و همکاران 2008).

جرم مخصوص ظاهری، رطوبت و ساختمان خاک بر مقاومت فروروی تأثیر زیادی دارند (بیرد و کاسل 1980). حتی اگر جرم مخصوص ظاهری خاک و رطوبت آن ثابت بماند باز مقاومت فروروی خاک با گذشت زمان تغییر می‌کند (پرفکت و همکاران 1992)، که شاید بتوان آن را به تغییرات تدریجی عوامل پیوندی خاکدانه‌ها در طی زمان مرتبط دانست. مقاومت فروروی با افزایش رطوبت خاک کاهش می‌یابد و با افزایش BD افزایش می‌یابد (کولن و کوپرس 1983). در برخی از خاک‌ها، تغییر در PR با BD و محتوای آب، خطی نیست (مارشال و همکاران 1996). مدل‌های گوناگون بسیاری (مانند مدل‌های خطی، توانی، چند جمله‌ای و نمایی) برای آزمون رابطه بینPR  با رطوبت و BD مورد استفاده قرار گرفته است (وز و همکاران 2001). آپادهایایا و همکاران (1982) با بررسی یک خاک لوم سیلتی، یک معادله نمایی را برای برآورد PR به عنوان تابعی از  GWCو BD ارائه کردند. همچنین آنها بررسی‌های بیشتری را برای ارزیابی درستی این مدل‌ها پیشنهاد کردند. داسیلوا و کی (1997) دریافتند که PR به مقدار رس و BD بستگی دارد. این پژوهشگران معادله‌ای را بدست آوردند که با ثابت بودن مقدار رس و با افزایش میزان آب خاک، اندازه پارامتر‌های BD و PR، کاهش می‌یابد. گرانوالد و همکاران (b2001) گزارش کردند که BD بیشترین تأثیر را در پیش بینی  PR دارد و پس از آن پارامتر GWC از اهمیت بیشتری برخوردار است. ساخت توابع انتقالی جدید دشوار است بنابراین منطقی است که از پارامتر‌هایی که در توابع پیشین ایجاد شده است، بهره‌گیری شود (مک­براتنی و همکاران 2002). بیات و همکاران (2008) با بررسی 381 نمونه خاک گردآوری شده از دو شهر همدان و مراغه، توانستند با بهره‌گیری از مدل‌های رگرسیونی و شبکه عصبی مصنوعی و به کمک پارامتر‌های BD،GWC ، TP وθv/TP  (درجه اشباع)، PR را برآورد کنند. آنها دریافتند که افزودن پارامترهای TP و θv/TP به مدل‌هایی که تنها دو پارامتر BD و GWC را دارا بودند، صحت مدل‌ها را بهبود بخشید. در این پژوهش تلاش شد تا افزون بر پارامتر‌های BD، رطوبت خاک و بافت خاک (که به گونه مرسوم برای برآورد مقاومت فروروی به کار می‌رود)، از پارامتر‌های فراکتالی توزیع اندازه ذرات و خاکدانه‌ها (خاکدانه‌های درشت (25/0 تا 8 میلی‌متر) و ریز (2میلی‌متر>) و همچنین از پارامتر‌های کیفی خاک مانند شکل، نوع و اندازه خاکدانه به عنوان ورودی‌های برآورد کننده مقاومت فروروی بهره‌گیری شود. 

به تازگی نظریه فراکتالی همانند ابزاری نیرومند برای مدل کردن پدیده‌های پیچیده بکار رفته است (تیلر و ویت کرفت 1992). این نگرش به پیشنهاد معادله‌هایی برای منحنی‌های مشخصه آب خاک،‌ هدایت هیدرولیکی غیر اشباع و اشباع خاک و توزیع اندازه منافذ درشت در خاک انجامیده است (شفیعی 1375). به تازگی بیات و همکاران (2011) به گونه کارآمدی از پارامتر‌های فراکتالی برای برآورد منحنی رطوبتی بهره‌گیری نمودند. مقیاس فراکتالی توزیع اندازه ذرات، بسیاری از فرایند‌های دینامیک و استاتیک در خاک، مانند انتقال آب و نمک­ها، گنجایش نگهداری آب، اندوخته گرمایی و هدایت گرمایی را کنترل می‌کند (ارساهین و همکاران 2006). احتمال می‌رود بهره‌گیری از پارامتر‌های فراکتالی باعث بهبود برآورد PR گردد. تا جایی که جستجو شد این پدیده تاکنون در هیچ پژوهشی بررسی و گزارش نشده است. بنابراین اهداف این پژوهش عبارتند از: 1. بررسی میزان بهبود برآورد PR با بهره‌گیری از پارامتر‌های کیفی خاک مانند شکل، نوع و اندازه خاکدانه 2. بررسی بهبود برآورد PR با بهره گیری از پارامتر‌های فراکتالی توزیع اندازه ذرات و توزیع اندازه خاکدانه‌ها (خاکدانه‌های ریز و درشت).

 

مواد و روش‌ها

در این پژوهش 53 نمونه خاک از خاک‌های شاخص بر پایه سری خاک‌های غالب از استان همدان نمونه برداری شد. نمونه‌های دست خورده از دو لایه رویین (10 تا 35 سانتی‌متری) و زیرین (20 تا 45 سانتی‌متری) خاک تهیه و پس از رساندن به آزمایشگاه، هوا خشک گردیده و برای برخی آزمایش‌های فیزیکی و شیمیایی آماده گردیدند. مقاومت فروروی به کمک دستگاه فروسنج دیجیتال مدل  Rimic CP20در محل نمونه برداری از خاک به گونه درجا اندازه‌گیری شد. نمونه جداگانه‌ای از خاک‌ها برای اندازه‌گیری رطوبت خاک بلافاصله پس از اندازه‌گیری PR برداشته شد. شکل، نوع و اندازه خاکدانه‌ها نیز در جایگاه نمونه‌برداری از خاک بر اساس راهنمای نمونه‌برداری و تشریح خاک‌ها (شونبرگر و همکاران 1998) یادداشت گردید. در آزمایشگاه میانگین وزنی قطر خاکدانه‌ها(MWD)[7] با روش الک تر (بیات و همکاران 1386)، توزیع اندازه خاکدانه‌های ریز (کوچکتر از 2 میلی‌متر) با روش هیدرومتر و الک تر (میلان و همکاران 2007) اندازه‌گیری شد. در واقع اندازه‌گیری خاکدانه‌های ریز مشابه توزیع اندازه ذرات (گی و ار 2002) است. با این تفاوت که از عوامل دیسپرس کننده استفاده نشده و به جای الک خشک برای تفکیک خاکدانه‌هایی در اندازه شن از الک تر استفاده می‌شود. توزیع اندازه ذرات (0 تا 2 میلی‌متر) به روش هیدرومتر و الک خشک (گی و ار 2002) و BD به روش کلوخه (کمپل و هنشال 2001) اندازه گیری شد. در این پژوهش از پارامتر‌های فراکتالی مدل‌های تیلر و ویت کرفت (1992)، بیرد و همکاران (2000) و مدل پریر و بیرد (2002 و 2003)، برای توزیع اندازه خاکدانه‌های درشت (25/0 تا 8 میلی‌متر)  و ریز (کوچک‌تر از 2 میلی‌متر) بهره‌گیری شد.

با توجه به اینکه تقسیم‌بندی استانداردی برای توزیع اندازه خاکدانه‌ها وجود ندارد، در این تحقیق بر اساس مطالعه میلان و همکاران (2007) خاکدانه‌های با قطر کمتر از 2 میلی‌متر به‌عنوان خاکدانه‌های ریز نامگذاری گردید. چون در این گروه، خاکدانه‌هایی در اندازه 2 میکرومتر نیز اندازه‌گیری می‌شوند. در اندازه‌گیری خاکدانه‌های درشت‌تر از روش کمپر و روزنو (1986) که توسط بیات و همکاران (1386) اصلاح شده بود و توزیع اندازه خاکدانه‌هایی به قطر 25/0 تا 8 میلی‌متر را اندازه‌گیری می‌کند، استفاده گردید. این گروه خاکدانه‌ها به‌صورت اختیاری تحت عنوان خاکدانه‌های درشت نامگذاری گردید.

برای توزیع اندازه ذرات نیز از مدل‌های فراکتالی بیرد و همکاران (2000) و مدل پریر و بیرد (2002 و 2003) بهره‌گیری شد.

شکل، نوع، اندازه و درجه وضوح خاکدانه‌ها که پارامتر‌هایی کیفی هستند نیز کمی شده و برای محاسبه دو پارامتر طول مؤثر مسیر انتشار[8] و طول مؤثر مسیر انتشار ASCALE[9] بهره‌گیری شدند (جارویس و همکاران 1997). به طور مختصر چگونگی محاسبه پارامترهای مذکور بدین ترتیب بود که سه پارامتر اندازه، شکل و درجه وضوح خاکدانه‌ها به کد‌هایی تبدیل شدند. به عنوان مثال برای خاکدانه‌های ریز عدد 1، متوسط عدد 2، درشت عدد 3 و خیلی درشت عدد 4 اختصاص داده شد. برای شکل و درجه وضوح خاکدانه‌ها نیز به همین ترتیب کد‌هایی داده شد. سپس از این کد‌ها استفاده کرده و دو پارامتر مذکور توسط نرم‌افزار ارائه شده به ‌وسیله جارویس و همکاران (1997) محاسبه گردید. پایداری خاک[10] بر اساس راهنمای نمونه‌برداری و تشریح خاک‌ها (شونبرگر و همکاران 1998) تعیین گردید. پایداری خاک یک ویژگی کیفی است و برای بررسی مقاومت خاک در مقابل تغییر شکل و گسستگی تعیین می‌گردد. با توجه به اینکه این ویژگی با میزان رطوبت خاک تغییر می‌کند می‌توان آن را در حالت‌های خشک، مرطوب و خیس تعیین نمود که در این تحقیق در حالت خشک تعیین گردید. در حالت خشک خاک‌ها از نظر مقاومت در 6 کلاس قرار می‌گیرند (شونبرگر و همکاران 1998): اگر خاک در برابر فشار هیچ‌گونه مقاومتی نداشته باشد در کلاس "سست[11]" قرار می‌گیرد. اگر خاک به آسانی و در برابر فشار جزئی به پودر و تک‌دانه تبدیل شود در کلاس "نرم[12]" قرار می‌گیرد. اگر خاک در بین انگشتان شست وسبابه به‌آسانی خرد شود در کلاس "کمی سخت[13]" قرار می‌گیرد. اگر خاک با دست‌ها به‌آسانی ولی در بین انگشتان شست و سبابه به سختی شکسته شود در کلاس "سخت[14]" قرار می‌گیرد. اگر خاک با فشار دست‌ها به‌ سختی شکسته شود در کلاس "خیلی سخت[15]" قرار می‌گیرد. اگر خاک با فشار دست‌ها شکسته نشود در کلاس "فوق‌العاده سخت[16]" قرار می‌گیرد. پس از تعیین کلاس پایداری هر خاک آن‌ها از طریق تبدیل به کد‌هایی کمی شدند. سپس دو پارامتر محاسبه شده به همراه پایداری خاک به عنوان ورودی برای برآورد مقاومت فروروی بکار رفتند.

محاسبه پارامتر‌های فراکتالی توزیع اندازه خاکدانه‌ها (خاکدانه‌های ریز و خاکدانه‌های درشت) و توزیع اندازه ذرات

با داشتن توزیع اندازه خاکدانه‌ها و توزیع اندازه ذرات از مدل‌های ذیل برای محاسبه ابعاد فراکتالی استفاده گردید. خاکدانه‌های درشت در پنج کلاس با اندازه 25/0 تا 5/0، 5/0 تا 1، 1 تا 2، 2 تا 75/4 و 75/4 تا 8 میلی متر تفکیک گردیدند.

مدل تیلر و ویت کرفت (1992)

[1]

 

 

Dm: بعد فراکتال جرمی، M(x<X): جرم تجمعی خاکدانه‌ها بر روی غربال‌ها با اندازه‌های کوچکتر از X، Mt: جرم کل خاکدانه‌ها (باقی مانده بر روی تمام غربال‌ها)، XL: بالاترین اندازه روزنه غربال که برابر با mm 8 برای خاکدانه‌های درشت و mm 2 برای خاکدانه‌های ریز است، x: میانگین قطر خاکدانه‌ها در هر کلاس.

Dm در فرمول فوق با استفاده از رگرسیون لگاریتمی داده‌های اندازه‌گیری شده از روش  الک خشک بدست می‌آید (گولسر 2006). ولی در این پژوهش توزیع اندازه خاکدانه‌ها که از روش الک تر بدست آمد در مدل فوق مورد استفاده قرار گرفت. در این مدل BD و شکل خاکدانه‌ها مستقل از مقیاس فرض می‌شود و قابل کاربرد بر ذرات خاک نیز می‌باشد.

مدل بیرد و همکاران (2000)

در این مدل، توزیع تجمعی جرم ذرات جامد مدل [17]PSF بر روی توزیع اندازه ذرات خاک [18](PSD) بکار رفت:

[2]

 

                            

: جرم یا حجم تجمعی عناصر از کوچک به بزرگ که‌ می‌تواند ذرات جامد و یا منافذ باشد. در مورد منافذ حجم تجمعی بکار می‌رود، di: حد بالایی اندازه عناصر در هر کلاس اندازه (منافذ یا ذرات)، D: بعد فراکتالی جرمی توزیع اندازه ذرات یا منافذ، c: ثابت مدل.

مدلPSF با فرض تناسب میان توزیع منافذ و ذرات جامد است و از آن برای محاسبه بعد فراکتالی توزیع اندازه ذرات و هم توزیع اندازه منافذ استفاده شده است. پریر و همکاران (1999) گزارش کردند که هم فاز جامد و هم منافذ را می‌توان با قانون توانی توزیع مدل‌سازی کرد.

مدل پریر و بیرد (2002 و 2003)

 

[3]                                      

 

M(x xi ) = (α D-3)(LD-3) xi3-D

: M(x xi ) جرم تجمعی خاک یا درصدجرمی خاک تشکیل شده از واحدهایی (مانند ذرات اولیه، ذرات خرد شده و خاکدانه‌های ریز) با قطر کوچکتر یا مساوی xi، α: ثابت، L: حد بالایی اندازه ذرات در هر کلاس، D: بعد فراکتال جرمی توزیع اندازه ذرات. یک نکته مهم اینکه، روابط اندازه و جرم فراکتالی بر پایه فرض یکسان بودن BD می‌باشند (میلان و همکاران 2007). توصیف کلیه علائم استفاده شده در این تحقیق در جدول 1 ذکر گردید.

در این پژوهش برای ساخت توابع انتقالی از روش رگرسیون غیر خطی بهره‌گیری شد. در آغاز همه متغیرها از دیدگاه نرمال بودن داده‌ها بررسی شدند و آنهایی که توزیع نرمال نداشتند نرمال و سپس استاندارد شدند. برای بررسی اثر پارامتر‌های فراکتالی ذرات، خاکدانه‌های ریز و درشت و پارامتر‌های ساختمانی چندین مدل ساخته شد که ورودی‌های مدل‌ها در جدول 2 آورده شده است.

آماره‌های ارزیابی توابع انتقالی

در پایان برتری و صحت هر یک از مدل‌ها بر پایه آماره‌هایR2 (ضریب تبیین) و R2adj(ضریب تبیین تعدیل شده)، مجذور میانگین مربعات خطا (RMSE[19])، معیار اطلاعات آکایک و ضریب بهبود نسبی (RI[20]) ارزیابی شد.

 

ضریب بهبود نسبی به صورت زیر محاسبه شد

 

[4]                    

RI: ضریب بهبود نسبی، RMSE1: مجذور میانگین مربعات خطای مدل پایه، RMSE2: مجذور میانگین

مربعات خطای مدل مورد نظر

معیار اطلاعات آکایک[21]، AIC، (آکایک 1974 )


 [5]                             

        

معیار اطلاعات آکایک تعداد پارامترهای مدل که باید برآورد شوند، np، را نیز وارد محاسبات می‌نماید. N تعداد نمونه‌های خاک، PRm مقاومت فروروی اندازه‌گیری شده و PRp مقاومت فروروی براورد شده می‌باشد. هرچه AIC کوچک‌تر باشد کیفیت برازش بالاتر است. اگر انطباق کامل مقادیر مشاهده‌ای و برآوردی رخ دهد AIC به سمت منفی بی نهایت میل خواهد کرد.

روش مقایسه با بهره‌گیری از آماره AIC به این گونه می‌باشد که نخست یکی از مدل‌ها همانند مدل پایه گزینش می‌شود و مدل‌های دیگر نسبت به آن سنجیده می‌شود. اندازه AIC برای مدل پایه و مدل‌های دیگر تعیین می‌شود. اگر اندازه به دست آمده از سنجش مدل پایه با مدل مقایسه شونده از 95/0 اندازه AIC به دست آمده برای مدل پایه کمتر باشد، به این معنی است که مدل مقایسه شونده با اطمینان 95 درصد بهتر از مدل پایه است و در غیر این صورت مدل پایه ترجیح داده می‌شود (هوانگ و همکاران 2002).

 

نتایج و بحث

کلیه علائم اختصاری استفاده شده د این تحقیق به همراه توصیف آن‌ها در جدول 1 آورده شده است. با توجه به اینکه در این تحقیق 18 تابع انتقالی ایجاد گردید، متغیر‌های ورودی هر تابع انتقالی در جدول 2 نشان داده شده است. متغیر وابسته برای کلیه توابع ایجاد شده PR بود. ویژگی‌های آماری کلیه متغیرها در جدول شماره 3 نشان داده شده است. این جدول شامل تمام متغیرهایی است که مستقیماً اندازه‌گیری شده و یا از طریق مدل‌های فراکتالی محاسبه گردیده‌اند. کلیه متغیرهای مربوط به بافت و ساختمان خاک و نیز پارامتر‌های فراکتالی محاسبه شده دارای دامنه تغییرات بالایی هستند.  

همبستگی خطی میان PR و متغیر‌های ورودی، در جدول 4 آورده شده است. بین PR و W، همبستگی منفی در سطح احتمال1%  وجود داشت. بین  PRو  BD همبستگی هر چند  مثبت بود اما معنی دار و چشم گیر نبود. بیات و همکاران (2008) نیز گزارش کردند که PR همبستگی مثبت و خیلی چشم گیر با BD در خاک‌های همدان داشت اما برای خاک‌های مراغه این همبستگی چشم گیر نبود. دلیل همبستگی مثبتPR  با BD و همبستگی منفی PR با  Wاین است که PR با افزایش محتوای آب خاک کاهش می‌یابد و با افزایش BD افزایش می‌یابد (کولن و کوپرس 1983).

مقاومت فروروی همبستگی منفی و چشم گیری در سطح احتمال 5% با شن و همبستگی مثبت و چشم گیری در سطح احتمال5% با میزان سیلت داشت (جدول 4). در یک مقدار ثابت از اندازه آب و BD، PR با کاهش اندازه ذرات تمایل به افزایش دارد (مارشال و همکاران 1996). دلیل این افزایش بالاتر بودن تنش مؤثر[22] در نتیجه پایین بودن پتانسیل ماتریک در خاکی با بافت ریزتر می‌باشد (مولینز و همکاران 1987).

نکته بارز اینکه PR همبستگی مثبت و چشم گیری نیز با ابعاد فراکتالی توزیع اندازه ذرات مدل‌های بیرد و پریر در سطح احتمال 1% داشت (جدول 4). یکی از دلایل همبستگی این پارامتر با ابعاد فراکتالی توزیع اندازه ذرات را می‌توان اینگونه توجیه کرد که با ریزتر شدن ذرات خاک بعد فراکتالی افزایش می‌یابد (تیلور و ویت کرفت 1992). در مقابل PR همبستگی چشم گیری با ابعاد فراکتالی توزیع اندازه خاکدانه‌های ریز و خاکدانه‌های درشت نداشت (جدول 4). شاید غیر خطی بودن همبستگی آنها یکی از دلایل یافته بالا باشد.

 

 

جدول 1-  توصیف علائم استفاده شده

متغیر

توصیف

متغیر

توصیف

Z1

طول مؤثر مسیر انتشار

C Psd B

پارامتر فراکتالی توزیع اندازه ذرات مدل بیرد

Z2

طول مؤثر مسیر انتشار ASCALE

D Psd B

بعد فراکتالی توزیع اندازه ذرات مدل بیرد

Z3

پایداری خاک

D micro M

بعد فراکتالی خاکدانه‌های ریز مدل تیلر و ویت کرافت

W

رطوبت جرمی در حین اندازه‌گیری مقاومت فروروی

α micro MI

پارامتر فراکتالی توزیع اندازه خاکدانه‌های ریزمدل پریر و بیرد

PR

مقاومت فروروی

D micro MI

بعد فراکتالی توزیع اندازه خاکدانه‌های ریز مدل پریر و بیرد

% Clay

درصد رس

C micro B

پارامتر فراکتالی توزیع اندازه خاکدانه‌های ریزمدل بیرد

% Silt

درصد سیلت

D micro B

بعد فراکتالی توزیع اندازه خاکدانه‌های ریز مدل بیرد

% Sand

درصد شن

D MWD M

بعد فراکتالی خاکدانه‌های درشت مدل تیلر و ویت کرافت

MWD

میانگین وزنی قطر خاکدانه‌ها

C MWD B

پارامتر فراکتالی توزیع اندازه خاکدانه‌های درشت مدل بیرد

BD

جرم مخصوص ظاهری

D MWD B

بعد فراکتالی توزیع اندازه خاکدانه‌های درشت مدل بیرد

α Psd MI

پارامتر فراکتالی اندازه ذرات مدل پریر و بیرد

α MWD MI

پارامتر فراکتالی  خاکدانه‌های درشت مدل پریر و بیرد

D Psd MI

بعد فراکتالی اندازه ذرات مدل پریر و بیرد

D MWD MI

بعد فراکتالی  خاکدانه‌های درشت مدل پریر و بیرد

 

 

در جدول 5 نتایج مدل‌های رگرسیونی آورده شده است. به طور کلی این یافته‌ها نشان داد که افزودن پارامتر‌های فراکتالی توزیع اندازه ذرات و توزیع اندازه خاکدانه‌های درشت و ریز می‌تواند باعث بهبود مدل‌های مرسوم برآورد کننده PR شود. در مدل‌هایی که از پارامتر‌های فراکتالی توزیع اندازه ذرات بهره‌گیری شد، اندازهR2 و R2adj افزایش یافته است، که این نتیجه می‌تواند به دلیل همبستگی بالای ابعاد فراکتالی توزیع اندازه ذرات با PR باشد (جدول 4). همانگونه که در جدول 5 دیده می‌شود در مدل چهارم با افزودن پارامتر‌های فراکتالی توزیع اندازه ذرات و پارامتر‌های کیفی به مدل، اندازه R2 و R2adj نسبت به مدل 1 افزایش معنی‌داری نشان داد (به ترتیب از 30/0 به 38/0 و از 18/0 به 28/0)، ضمن اینکه اندازه RMSE کاهش یافت (از 82/0 به 77/0).

هر چند که تا کنون در مورد استفاده از پارامتر‌های فراکتالی توزیع اندازه ذرات برای تخمین PR گزارشی منتشر نشده است، ولی هورن (1984) گزارش کرد که PR تحت تاثیر توزیع اندازه ذرات خاک قرار دارد. بنابراین با توجه به کمی شدن توزیع اندازه ذرات خاک توسط پارامتر‌های فراکتالی، نتیجه حاصله مورد انتظار بود.

 

 

 

 

 

 

 

 

 

 

ورودی‌های مدل

مدل

ورودی‌های مدل

مدل

تمام ابعاد مدل‌های پیشین به علاوه رطوبت و جرم مخصوص ظاهری                                                    (D micro M ,DMWD M, W, Pb ,C Psd B, D Psd B)

10

رطوبت و جرم مخصوص ظاهری±(W, Pb)      

1

 بعد و پارامتر فراکتالی پریر و بیرد توزیع اندازه ذرات و ورودی‌های مدل اول                                                               (W, Pb , α Psd MI, DPsd MI )                                                                                                                                                                 

11

پارامتر‌های فراکتالی توزیع اندازه ذرات مدل بیرد (C Psd B, D Psd B )

2

غالب ابعاد فراکتالی مدل‌های قبلی و ورودی‌های مدل اول

(Dmicro M, D MWD M, C Psd B, W, Pb, a Psd MI, D Psd MI)

12

رطوبت ،جرم مخصوص ظاهری، پارامتر‌های کیفی ( Z1 ,Z2, Z3, W, Pb)          

3

پارامتر‌های فراکتالی توزیع اندازه ذرات و خاکدانه‌های درشت مدل پریر و بیرد و پارامتر‌های مدل1 (α Psd MI , D psd MI,   α MWD MI,D MWD MI ,W.Pb)

13

تمام پارامتر‌های مدل 2و3                           ( Z1 ,Z2, Z3, W, Pb ,C Psd B, D Psd B)              

4

پارامتر‌های فراکتالی توزیع اندازه ذرات و خاکدانه‌های ریز مدل پریر و بیرد و پارامتر‌های مدل1

(α psd MI , D psd MI ,α micro MI,D micro MI  ,W.Pb )

14

رطوبت، جرم مخصوص ظاهری، شن و رس                             (W, Pb ,Clay,Sand)

5

پارامتر‌های فراکتالی توزیع اندازه خاکدانه‌های درشت و خاکدانه‌های ریز مدل پریر و بیرد و پارامتر های مدل1    (α micro MI,D micro MI,   α MWD MI,D MWD MI,W.Pb)

15

بعد فراکتالی خاکدانه‌های ریز  مدل تیلر و ویت کرفت و پارامتر‌های مدل1                                            (D micro M , W, Pb )

6

پارامتر‌های فراکتالی توزیع اندازه ذرات و خاکدانه‌های درشت مدل بیرد و پارامتر‌های مدل1

(α MWD MI, D MWD MI, W, Pb, α Psd MI’ D Psd MI)

16

  بعد فراکتالی تیلر و ویت کرفت خاکدانه‌های درشت، MWD و پارامتر‌های مدل 1                                             ( W, Pb ,MWD,  DMWD M  ) 

7

پارامتر‌های فراکتالی توزیع اندازه ذرات و خاکدانه‌های ریز مدل بیرد و پارامتر‌های مدل1

(α psd B , D psd B,α micro B,D micro B,W.Pb)

17

بعد فراکتالی تیلر و ویت کرفت خاکدانه‌های ریز و تمام پارامتر‌های مدل 1و 2                                         (D micro M, W, Pb ,C Psd B, D Psd B)

8

پارامتر‌های فراکتالی توزیع اندازه خاکدانه‌های درشت و خاکدانه‌های ریز مدل بیرد و پارامتر‌های مدل 1(αmicro B,D micro B,αMWD B,D MWD B,W.Pb )

18

بعد فراکتالی تیلر و ویت کرفت خاکدانه‌های درشت و تمام پارامتر‌های مدل 1و 2                                         (DMWD M, W, Pb ,C Psd B, D PsdB) 

9

جدول 2- ورودی مدل‌ها

 ± .توصیف علائم در جدول 1 ذکر شده است.

 

 

 

در مدل دهم افزودن پارامتر‌های فراکتالی توزیع اندازه ذرات (مدل بیرد) و توزیع اندازه خاکدانه‌های ریز و خاکدانه‌های درشت (مدل تیلر و ویت کرافت) نیز توانسته میزان R2 و  R2adjرا نسبت به مدل 1به خوبی افزایش دهد (به ترتیب از 30/0 به 38/0 و از 18/0 به 30/0) که می‌تواند وابسته به رابطه بعد فراکتالی خاکدانه‌ها با فرایند خاکدانه‌سازی باشد. در واقع پارامتر‌های فراکتالی توزیع اندازه خاکدانه‌ها (ریز و درشت) ساختمان خاک را کمی می‌کنند. بنگوق و همکاران (2001) نیز گزارش کردند که PR تحت تاثیر ساختمان خاک و BD است. بنابراین نتایج این تحقیق در راستای نتایج محققان فوق است . چرا که با نگاهی به آماره آکایک در جدول 5 می توان این نکته را دریافت که با افزودن ابعاد فراکتالی توزیع اندازه ذرات و خاکدانه‌ها به مدل پایه، در اکثر مدل‌ها بهبود برآورد PR در سطح احتمال 5% معنی دار­ بوده است. در مدل‌های دوم، چهارم، هشتم، نهم، دهم، دوازدهم، سیزدهم، چهاردهم، شانزدهم و هفدهم، افزوده شدن ابعاد فراکتالی به مدل پایه باعث بهبود معنی دار دقت مدل‌ها در سطح احتمال 5% شده است. در مدل‌های چهارم، دهم و دوازدهم به ترتیب میزان آماره آکایک از 39/12- در مدل اول به 60/18- ، 60/18- و 84/19- کاهش یافته است که بهبود قابل توجه و معنی‌دری را نشان می‌دهد. زیرا افزایش بعد فراکتالی میزان رس بالای خاک را نشان می‌دهد که خود این رس جزء ترکیبات مهم در پیدایش خاکدانه‌ها است. ضریب بهبود نسبی (RI) نیز در مدل‌های شماره دو، چهار، هشت، نه، ده، دوازدهم، سیزده، چهارده، شانزده و هفده افزایش داشت که نشان از کاهش مجذور میانگین مربعات خطا در این مدل‌ها نسبت به مدل اول دارد. در مدل چهارم ضریب بهبود نسبی (7/5%) نشان از بهبود این مدل با افزودن پارامترهای کیفی و پارامتر‌های فراکتالی مدل بیرد دارد.

 

 

جدول 3-  ویژگی­های متغیرها

 

متغیر

میانگین

میانه

انحراف استاندارد

حد اقل

حداکثر

متغیر

میانگین

میانه

انحراف استاندارد

حد اقل

حداکثر

Z1±

45/5

5

74/1

2

10

C Psd B

09/1

1/1

07/0

95/0

22/1

Z2

43/7

5

34/4

2

20

D Psd B

76/2

78/2

1/0

51/2

92/2

Z3

4/3

3

08/1

1

5

D micro M

78/2

8/2

07/0

58/2

9/2

W ( g g-1)

63/6

97/5

79/3

76/1

7/16

α micro MI

5/0

48/0

14/0

28/0

94/0

PR ( kPa )

31/1295

1210

43/635

11/201

2610

D micro MI

73/2

75/2

1/0

48/2

88/2

Clay (%)

06/22

20

2/12

89/5

9/63

C micro B

03/1

02/1

08/0

87/0

28/1

Silt (%)

37/38

26/42

03/16

95/0

9/61

D micro B

72/2

73/2

1/0

48/2

88/2

Sand (%)

61/39

42/35

33/19

61/9

7/81

D MWD M

7/3

8/3

29/0

72/2

8/3

MWD (mm)

6/2

22/2

7/1

21/0

54/5

c MWD B

14/0

15/0

1/0

001/0

34/0

BD ( Mg m-3)

6/1

61/1

1/0

36/1

76/1

D MWD B

28/2

67/2

88/0

24/0

96/2

α Psd MI

36/0

34/0

14/0

15/0

94/0

α MWD MI

71/7

68/1

27/13

11/0

12/11

D Psd MI

76/2

78/2

1/0

51/2

92/2

D MWD MI

28/2

67/2

88/0

23/0

95/2

.± توصیف علائم در جدول 1 ذکر شده است.

 

 

 

از مدل‌های 13 تا 18 این نکته دریافت شد که هرگاه از پارامتر‌های توزیع اندازه ذرات به همراه پارامتر‌های فراکتالی توزیع اندازه خاکدانه‌های ریز و درشت با هم، برای برآورد PR بهره‌گیری شود میزان آماره‌هایR2، R2adj و RI افزایش می‌یابد. اما هرگاه از پارامتر‌های فراکتالی توزیع اندازه خاکدانه‌های درشت و ریز با یکدیگر، برای برآورد PR بهره‌گیری شود اندازه آماره‌های یاد شده بهبودی نیافتند. شاید دلیل این نتیجه این باشد که PR متأثر از بافت و ساختمان خاک است (بنگوق و همکاران 2001). به همین دلیل استفاده همزمان از پارامتر‌های فراکتالی توزیع اندازه ذرات و خاکدانه‌ها موجب بهبود تخمین PR گردید.

پارامتر‌های فراکتالی توزیع اندازه ذرات و توزیع اندازه خاکدانه‌ها موجب افزایش صحت و قابلیت اعتماد مدل‌ها می‌گردند و کاهش خطای برآورد را به همراه دارند. احتمالاً دلیل آن قابلیت مدل‌های فراکتالی در مدل‌سازی PSD و توزیع اندازه خاکدانه‌ها است. داسیلوا و کی (1997) نیز PR را از طریق یک رابطه رگرسیونی به‌صورت تابعی از BD و رطوبت خاک به‌دست آوردند. آن‌ها ضرایب مدل خود را توسط توابع انتقالی با استفاده از جزء رس و ماده آلی تخمین زدند. بنابراین همانگونه که در این تحقیق به‌دست آمد روابط آن‌ها حاکی از آن است که بافت و ساختمان خاک عوامل اصلی مؤثر بر PR می‌باشند. بهره‌گیری از پارامتر‌های فراکتالی PSD و توزیع اندازه خاکدانه‌ها برای برآورد PR توصیه‌ می‌گردد. با مقایسه آماره‌های مدل 5 (جدول 5) که شامل اجزاء بافت خاک هستند با آماره‌های مدل‌هایی که از پارامتر‌های فراکتالی توزیع اندازه ذرات بهره‌گیری شده است (مانند مدل‌های 2، 4 و 14) این نکته دریافت می‌شود که پارامتر‌های فراکتالی توزیع اندازه ذرات، اطلاعات کامل‌تری از بافت خاک را در اختیار مدل قرار می‌دهد و نسبت به اجزاء بافت خاک پارامتر قابل اطمینان تری می‌باشند. یوسوویکز و لیپیک (2009) نیز تاثیر تراکم بر توزیع مکانی PR را توسط ابعاد فراکتالی مورد بررسی قرار دادند و مشاهده کردند که در جهت عمودی با افزایش تراکم بعد فراکتالی نیز افزایش می‌یابد. بنابراین نتیجه‌گیری کردند که بعد فراکتالی شاخص مناسبی برای ارزیابی سطح تراکم خاک می‌باشد.

مدل تیلر و ویت کرفت به تنهایی باعث بهبود اندازه آماره‌های یاد شده نشد و نتوانست شاخص خوبی برای توزیع اندازه خاکدانه‌های درشت و ریز باشد، که شاید یکی از دلایل این یافته، تک پارامتر‌ی بودن این مدل نسبت به دو مدل دیگر (که هر دو دارای دو پارامتر برآورد شده هستند)، باشد.

 

نتیجه‌گیری کلی

 افزودن پارامتر‌های فراکتالی توزیع اندازه ذرات مدل بیرد و مدل پریر و بیرد به ورودی‌های مرسوم برآورد کننده PR باعث بهبود معنی‌دار معیارهای ارزیابی صحت مدل‌ها مانند R2، RI،  R2adjو RMSE گردید. افزودن پارامتر‌های کیفی خاک به همراه پارامتر‌های توزیع اندازه ذرات نیز می‌تواند آماره‌های یاد شده را بطور معنی‌داری بهبود دهد. همچنین  می‌توان نتیجه گرفت که هرگاه از پارامتر‌های توزیع اندازه ذرات با پارامتر‌های توزیع اندازه خاکدانه‌های درشت یا ریز به عنوان ورودی بهره‌گیری شود، برآورد PR بهبود معنی‌داری پیدا می‌کند. با بررسی و مقایسه نوع مدل‌های بکار رفته در این پژوهش نیز می‌توان دریافت که پارامتر‌های فراکتالی توزیع اندازه ذرات و توزیع اندازه خاکدانه‌های ریز و درشت مدل بیرد، نسبت به مدل پریر و بیرد در بیشتر مدل‌ها برتری محسوسی دارد. با افزودن مدل تیلر و ویت‌کرفت بهبودی در مدل‌ها ایجاد نشد و این مدل برای برآورد PR پیشنهاد نمی‌شود. از آن جایی که مدل‌های فراکتالی بسیاری در پژوهش‌های گوناگون برای توزیع اندازه ذرات و خاکدانه‌ها ارائه شده است، بنابراین بهتر است که از دقیق‌ترین این مدل‌ها در چنین مطالعاتی بهره‌گیری شود.


 

جدول 4-همبستگی خطی مقاومت فروروی با 23 متغیر مدل­های رگرسیونی

±% Sand

% Silt

% Clay

W

Z3

Z2

Z1

PR

 

*33/0-

*30/0

12/0

*30 /0-

13/0

05/0-

08/0

-

PR

BD

D Psd B

C Psd B

D micro M

DMWD M

α Psd MI

MWD

DPsd MI

 

14/0

**36/0

18/0

07/0

14/0

25/0

06/0

**36/0

PR

D MWD MI

α  MWD MI

D MWD B

C MWD B

D micro B

C micro B

D micro MI 

α micro MI

 

04/0-

15/0-

04/0-

11/0-

23/0

15/0

25/0

08/0

PR

±کلیه علائم اختصاری در جدول 1 آورده شده است.

و  ** به ترتیب معنی دار در سطوح احتمال 5 و 1 درصد.

 

جدول 5- مدل‌های برازش شده بر پارامتر‌های ورودی و آماره های مربوط به هر مدل

 

 

AIC

39/12-

93/17-*

34/29

60/18- *

39/12-

46/2-

35/4-

55/13-*

65/15-*

60/18-*

62/11-

84/19-*

65/15-*

33/14-*

22/9-

98/16-*

65/15-*

76/6-

±ترتیب متغیرهای x1 تا xn در مدل‌ها، در جدول 2 به ترتیب از چپ به راست، برای هر شماره مدل آورده شده است.

 

 

R I

 

5%

7/1-%

7/5%

1/0-%

8/9-%

7-%

1/1%

1/3%

6/5%

8/0-%

7/6%

5/2%

2%

3-%

4%

2/3%

5-%

آمارهها

RMSE

825/0

783/0

839/0

778/0

825/0

906/0

890/0

816/0

80/0

778/0

831/0

769/0

80/0

81/0

85/0

79/0

80/0

87/0

 

Ra2

181/0

226/0

207/0

287/0

248/0

113/0

141/0

250/0

272/0

30/0

238/0

285/0

25/0

25/0

17/0

25/0

27/0

14/0

 

R2

307/0

375/0

283/0

383/0

306/0

164/0

190/0

322/0

342/0

38/0

297/0

395/0

34/0

33/0

26/0

36/0

35/0

24/0

مدل‌ها (روابط رگرسیونی)±

Model

a+b/x1+c/x12+d/x13+e/x14+f×x2+g×x22+h×x23+i×x24

a+b/x1+c/x12+d/x13+e/x14+ f/x15+g×x2+h×x22+i×x23+j×x24+k×x25

exp(a×x1+b×x2+c×x3+d×x4+e×x5+f)

a×x1+b×x2+c×x3+d×x4+e×x5+f×x6+g×x7+h

a×x1+b×x2+c×x3+d×x4+e

a×x1+b×x2+c×x3+d

a×x1+b×x2+c×x3+d×x4+e

a×x1+b×x2+c×x3+d×x4+e×x5+f

a×x1+b×x2+c×x3+d×x4+e×x5+f

a×x1+b×x2+c×x3+d×x4+e×x5+f×x6+g

a×x1+b×x2+c×x3+d×x4+e

a×x1+b×x2+c×x3+d×x4+e×x5+f×x6+g×x7 +h

a×x1+b×x2+c×x3+d×x4+e×x5+f×x6+g

exp(a×x1+b×x2+c×x3+d×x4+e×x5+f×x6+g)

a×x1+b×x2+c×x3+d×x4+e×x5+f×x6+g

a×x1+b×x2+c×x3+d×x4+e×x5+f×x6+g

a×x1+b×x2+c×x3+d×x4+e×x5+f×x6+g

a×x1+b×x2+c×x3+d×x4+e×x5+f×x6+g

 

k

 

19/0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

80/2-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

11/1-

17/15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

29/2

95/37-

 

83/0-

 

 

 

 

 

 

 

76/0

 

 

 

 

 

 

 

g

19/1-

18/43

 

40/0

 

 

 

 

 

62/5-

 

09/0-

53/9-

43/9-

23/9-

88/12-

16/10-

28/9-

مدل‌ها

f

001/0-

0005/0-

59/1-

03/0-

 

 

 

58/4-

27/1-

39/0

 

27/0

12/0

18/0

09/0

14/0

16/0

12/0

ضرایب

e

001/0-

01/0

35/0

09/0

88/0

 

04/0-

39/0

19/0

16/0

83/0-

29/0-

37/0-

29/0-

41/0-

37/0-

35/0-

36/0-

 

d

04/0

01/0

46/1-

44/0-

58/0-

97/2-

18/0

09/0

19/0

16/0

32/0

44/0-

13/0-

39/3-

06/0

16/0

06/4-

14/0

 

c

04/0

42/0-

16/1

29/0-

30/0-

22/0

05/0

14/0

15/0

40/0-

15/0-

22/2

001/0

33/1

00/0

97/3-

91/3-

73/2-

 

b

31/0-

51/0

17/0-

29/0-

20/0

35/0-

23/0

36/0-

37/0-

26/0

33/0-

42/0

53/3

36/6

61/3

38/3

61/7

06/3

 

a

01/0-

42/17-

38/0

27/0

36/0-

15/0

34/0-

17/0

23/0

21/0

33/0-

11/0

92/0-

46/1-

10/0-

77/3

12/4

35/1

 

شماره مدل

1

2

 

3

 

 

 

 

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

منابع مورد استفاده

بیات ح، محبوبی ع ا، حاج عباسی م ع و مصدقی م ر، 1386. اثر سیستم‌های خاک ورزی و انواع ماشین‌های کشاورزی بر جرم مخصوص ظاهری، شاخص مخروطی و پایداری ساختمان یک خاک لوم شنی. علوم و فنون کشاورزی و منابع طبیعی. شماره 42 (ب). صفحه‌های 451 تا 461.

شفیعی ع، 1375. مطالعه فراکتال‌ها و کاربرد آنها. پایان نامه کارشناسی ارشد، دانشکده مهندسی مکانیک دانشگاه صنعتی اصفهان.

مسکینی ف، 1389. تجزیه علیت خصوصیات خاک موثر بر هدایت هیدرولیکی اشباع خاک‌های شالیزاری. نشریه آب و خاک، جلد 24، شماره6، صفحه‌های 1246 تا 1253.

Akaike H, 1974. A new look at the statistical model identification. IEEE Transactions on Automatic Control AC-19:716-723.

Bayat H,  Neyshabouri MR,  Hajabbasi MA,  Mahboubi AA and Mosaddeghi MR,  2008. Comparing neural networks, linear and nonlinear regression techniques to model penetration resistance. Turk J Agric For  32: 425-433.

Bayat H, Neyshabouri MR, Mohammadi K and Nariman-Zadeh N, 2011. Estimating water retention with pedotransfer functions using multi-objective group method of data handling and ANNs. Pedosphere 21: 107-114.

Bengough AG, Campbell DJ and Osallivan MF, 2001. Penetrometer techniques in relation to soil compaction and root growth. Pp 377-404. In: Smith KA and Mullins CA (eds). Soil and EnvironmentalAnalysis. Marcel Dekker, Inc. New York • Basel.

Byrd CW and Cassel DK, 1980. The effect of sand content upon cone index and selected physical properties. Soil Sci 129 : 197-204.

Bird NRA, Perrier E and Rieu M, 2000. The water retention function for a model of soil structure with pore and solid fractal distributions. Eur.J Soil Sci 51: 55-63.

Busscher WJ, 1990. Adjustment of flat-tipped penetrometer resistance data to a common water content. Trans  ASAE  33: 519-524.

Campbell D J, 1975. Liquid limit determination of arable topsoils using a drop-cone penetrometer.J Soil Sci 26:234 –240.

Campbell DJ and Henshall JK, 2001. Bulk Density. Pp 315-348. In: Smith KA and Mullins CE (eds). Soil and Environmental Analysis. Marcel.

Dekker, Inc. New York

Da Silva AP and Kay  BD, 1997. Estimating the least limiting water range of soils from properties and management. Soil Sci  Soc  Am  J  61: 877-883.

Ersahin S, Gunal H, Kutlu T, Yetgin B and  Coban S, 2006. Estimating specific surface area  and cation exchange capacity in soils using fractal dimension of particle-size distribution. Geoderma 136: 588-597

Faure AG and  da Mata  JDV, 1994. Penetration resistance value along compaction curves. J. Geotech Eng 120: 46-59.

Gee GW and Or D,  2002. Particle – Size Analysis. Pp. 255-295. In: Warren AD (ed). Methods of Soil Analysis. Part 4. Physical Methods. SSSA Inc, Madison WT.

Gerard C J, 1965. The influence of soil moisture, soil texture, drying conditions and exchangeable cations on soil strength. Soil Sci Soc Am Proc 29: 641– 645.

Grunwald S, Lowery B, Rooney DJ and McSweeney K, 2001a. Profile cone penetrometer data used to distinguish between soil materials. Soil Till Res 62: 27-40.

Grunwald S, Rooney DJ, McSweeney K., and Lowery  B, 2001b. Development of pedotransfer functions for a profile cone penetrometer. Geoderma 100: 25-47.

Gulser C, 2006. Effect of forage cropping treatment on soil structure and  relationship with fractal  dimensions. Geoderma131:33-44.

Henderson C, Levett A, and Lisle D, 1988. The effect of soil water content and bulk density on the compactability and soil penetration resistance of some Western Australian sandy soils. Aust J Soil Res 26: 391-400.

Horn R, 1984. The prediction of the penetration resistance of soils by multiple regression analysis. Z Kulturtechn. Flurbereinig 25: 377–380.

Hwang SIl, Kwang PL, Dong SL and Susan EP, 2002. Models for estimating soil particle-size distributions. SSAJ. 66:1143–1150.

Jarvis NJ, Hollis  JM, Nicholls PH, Mayer  T and Evans SP, 1997. Macro – DB : a decision – support tool for assessing pesticide fate and mobility in soils. Environmental Modeling and Software 12: 251-265.

Kemper WD and Rosenau RC, 1986. Aggregate stability and size distribution. Pp: 425-442. In: Klute A (ed). Methods of Soil Analysis. Part 1. 2nd ed ASA and 385A Madison WI. 

Koolen  AJ and  Kuipers  H, 1983. Agricultural Soil Mechanics. Heidelberg: Springer- Verlag.

Kurup PU, Voyiadjis  GZ and  Tumay  MT, 1994. Calibration chamber studies of piezocone test in cohesive soils. J Geotech Eng 120: 81-107.

Marshall TJ, Holmes JW and Rose CW, 1996. Soil Physics. 3rd ed. Cambridge, Cambridge University Press.

McBratney AB, Minasny  B, Cattle SR and Vervoort  RW, 2002. From pedotransfer functions to soil inference systems. Geoderma 109:41-73.

Millan H, Gonzalez-Posada M, Morilla AA and Perez, E, 2007. Self-similar organization of Vertisol microstructure: A pore–solid fractal interpretation. Geoderma 138:  185–190.

Mullins CE, Young IM, Bengough AG, and Ley GJ, 1987. Hard-setting soils. Soil Use Manag. 3:79–83.

Nemes A, Schaap MG and Wosten JHM,  2003. Functional evaluation of pedotransfer functions derived from different scales of data collection. Soil Sci Soc Am J 67: 1093-1102.

Osullivan MF and Ball BC, 1982. A comparison of five instruments for measuring soil strength in cultivated and uncultivated cereal seedbeds. Soil Sci J 33 : 597-608.

Perrier E, Bird N and  Rieu M, 1999. Generalizing the fractal model of soil structure: the pore-solid fractal approach. Geoderma 88: 137-164.

Perrier E, Bird N, 2002. Modelling soil fragmentation: the PSF approach. Soil Till Res 64: 91–99.

Perrier EMA, Bird NRA, 2003. The PSF model of soil structure: a multiscale approach. Pp. 1–18. In: Pachepsky Ya, Radcliffe DE, Selim HM (eds). Scaling Methods in Soil Physics. CRC Press, Boca Raton, FL,

Perfect E, Rasiah V and  Kay BD, 1992. Fractal dimensions of soil aggregate-size distributions  calculated by number and mass. Soil Sci Soc  Am  J 56: 1407-1409.

Pidgeon JD and Soane  BD, 1977. Effects of tillage and direct drilling on soil properties during the growing season in a long-term barley mono-culture system. J Agric Sci 88: 431-442.

Puppala AJ, Acar  YB and Tumay MT, 1995. Cone penetration in very weakly cemented sand. J Geotech  Eng  121: 589-600.

Schoeneberger PJ, Wysocki DA, Benham EC and Brooderson DE, 1998. Field book for describing and sampling soils. Natural Resourses Conservation Service, USDA, Natural Soil Survey Center, Lincoln, NE.

To J and Kay  BD, 2005. Variation in penetrometer resistance with soil properties: the contribution of effective stress and implications for pedotransfer functions. Geoderma 126: 261-276.

Tyler SW and Wheatcraft SW, 1990. Fractal process in soil water retention. Water Resour Res 26: 1047-1054.

Tyler SW and Wheatcraft SW, 1992. Fractal scaling of soil particle-size distribution: analysis and limitations. Soil Sci  Soc  Am  J  56: 362-369.

Upadhyaya SK, Kemble  LJ and  Collins  NE, 1982. Cone index prediction equations for Delaware soils. ASAE paper. 82: 1452-1456.

Usowicz B, and Lipiec J, 2009. Spatial distribution of soil penetration resistance as affected by soil compaction: The fractal approach. Ecological Complexity 6: 263–271.

Vaz CMP, Luis HB and  Hopmans  JW, 2001. Contribution of water content and bulk density to field soil penetration resistance as measured by a combined cone penetrometer-TDR probe. Soil Till Res 60: 35-42.

Wosten JHM, Pachepsky Y and Rawls WJ, 2001. Pedotransfer functions: bridging the gap between available basic soil data and missing soil hydraulic characteristics. J Hydrol 251: 123-150.

 

 



[1] Penetration  resistance

[2] Parameter

[3] Pedotransfer functions

[4] Bulk density

[5] Gravimetric water content

[6] Total porosity

[7]  Mean weight diameter

1 Effective diffusion pathlength (mm)

[9] Effective diffusion pathlength ASCALE (mm)

[10] Soil consistency

[11]  Loose

[12]  Soft

[13]  Slightly hard

[14]  Hard

[15]  Very hard

[16]  Extremely hard

[17]Pore-solid fractal

[18] Particle size distribution

[19] Root mean square error

[20] Relative improvement

[21] Akaike  information criterion

[22] Effective stress

بیات ح، محبوبی ع ا، حاج عباسی م ع و مصدقی م ر، 1386. اثر سیستم‌های خاک ورزی و انواع ماشین‌های کشاورزی بر جرم مخصوص ظاهری، شاخص مخروطی و پایداری ساختمان یک خاک لوم شنی. علوم و فنون کشاورزی و منابع طبیعی. شماره 42 (ب). صفحه‌های 451 تا 461.
شفیعی ع، 1375. مطالعه فراکتال‌ها و کاربرد آنها. پایان نامه کارشناسی ارشد، دانشکده مهندسی مکانیک دانشگاه صنعتی اصفهان.
مسکینی ف، 1389. تجزیه علیت خصوصیات خاک موثر بر هدایت هیدرولیکی اشباع خاک‌های شالیزاری. نشریه آب و خاک، جلد 24، شماره6، صفحه‌های 1246 تا 1253.
Akaike H, 1974. A new look at the statistical model identification. IEEE Transactions on Automatic Control AC-19:716-723.
Bayat H,  Neyshabouri MR,  Hajabbasi MA,  Mahboubi AA and Mosaddeghi MR,  2008. Comparing neural networks, linear and nonlinear regression techniques to model penetration resistance. Turk J Agric For  32: 425-433.
Bayat H, Neyshabouri MR, Mohammadi K and Nariman-Zadeh N, 2011. Estimating water retention with pedotransfer functions using multi-objective group method of data handling and ANNs. Pedosphere 21: 107-114.
Bengough AG, Campbell DJ and Osallivan MF, 2001. Penetrometer techniques in relation to soil compaction and root growth. Pp 377-404. In: Smith KA and Mullins CA (eds). Soil and EnvironmentalAnalysis. Marcel Dekker, Inc. New York • Basel.
Byrd CW and Cassel DK, 1980. The effect of sand content upon cone index and selected physical properties. Soil Sci 129 : 197-204.
Bird NRA, Perrier E and Rieu M, 2000. The water retention function for a model of soil structure with pore and solid fractal distributions. Eur.J Soil Sci 51: 55-63.
Busscher WJ, 1990. Adjustment of flat-tipped penetrometer resistance data to a common water content. Trans  ASAE  33: 519-524.
Campbell D J, 1975. Liquid limit determination of arable topsoils using a drop-cone penetrometer.J Soil Sci 26:234 –240.
Campbell DJ and Henshall JK, 2001. Bulk Density. Pp 315-348. In: Smith KA and Mullins CE (eds). Soil and Environmental Analysis. Marcel.
Dekker, Inc. New York
Da Silva AP and Kay  BD, 1997. Estimating the least limiting water range of soils from properties and management. Soil Sci  Soc  Am  J  61: 877-883.
Ersahin S, Gunal H, Kutlu T, Yetgin B and  Coban S, 2006. Estimating specific surface area  and cation exchange capacity in soils using fractal dimension of particle-size distribution. Geoderma 136: 588-597
Faure AG and  da Mata  JDV, 1994. Penetration resistance value along compaction curves. J. Geotech Eng 120: 46-59.
Gee GW and Or D,  2002. Particle – Size Analysis. Pp. 255-295. In: Warren AD (ed). Methods of Soil Analysis. Part 4. Physical Methods. SSSA Inc, Madison WT.
Gerard C J, 1965. The influence of soil moisture, soil texture, drying conditions and exchangeable cations on soil strength. Soil Sci Soc Am Proc 29: 641– 645.
Grunwald S, Lowery B, Rooney DJ and McSweeney K, 2001a. Profile cone penetrometer data used to distinguish between soil materials. Soil Till Res 62: 27-40.
Grunwald S, Rooney DJ, McSweeney K., and Lowery  B, 2001b. Development of pedotransfer functions for a profile cone penetrometer. Geoderma 100: 25-47.
Gulser C, 2006. Effect of forage cropping treatment on soil structure and  relationship with fractal  dimensions. Geoderma131:33-44.
Henderson C, Levett A, and Lisle D, 1988. The effect of soil water content and bulk density on the compactability and soil penetration resistance of some Western Australian sandy soils. Aust J Soil Res 26: 391-400.
Horn R, 1984. The prediction of the penetration resistance of soils by multiple regression analysis. Z Kulturtechn. Flurbereinig 25: 377–380.
Hwang SIl, Kwang PL, Dong SL and Susan EP, 2002. Models for estimating soil particle-size distributions. SSAJ. 66:1143–1150.
Jarvis NJ, Hollis  JM, Nicholls PH, Mayer  T and Evans SP, 1997. Macro – DB : a decision – support tool for assessing pesticide fate and mobility in soils. Environmental Modeling and Software 12: 251-265.
Kemper WD and Rosenau RC, 1986. Aggregate stability and size distribution. Pp: 425-442. In: Klute A (ed). Methods of Soil Analysis. Part 1. 2nd ed ASA and 385A Madison WI. 
Koolen  AJ and  Kuipers  H, 1983. Agricultural Soil Mechanics. Heidelberg: Springer- Verlag.
Kurup PU, Voyiadjis  GZ and  Tumay  MT, 1994. Calibration chamber studies of piezocone test in cohesive soils. J Geotech Eng 120: 81-107.
Marshall TJ, Holmes JW and Rose CW, 1996. Soil Physics. 3rd ed. Cambridge, Cambridge University Press.
McBratney AB, Minasny  B, Cattle SR and Vervoort  RW, 2002. From pedotransfer functions to soil inference systems. Geoderma 109:41-73.
Millan H, Gonzalez-Posada M, Morilla AA and Perez, E, 2007. Self-similar organization of Vertisol microstructure: A pore–solid fractal interpretation. Geoderma 138:  185–190.
Mullins CE, Young IM, Bengough AG, and Ley GJ, 1987. Hard-setting soils. Soil Use Manag. 3:79–83.
Nemes A, Schaap MG and Wosten JHM,  2003. Functional evaluation of pedotransfer functions derived from different scales of data collection. Soil Sci Soc Am J 67: 1093-1102.
Osullivan MF and Ball BC, 1982. A comparison of five instruments for measuring soil strength in cultivated and uncultivated cereal seedbeds. Soil Sci J 33 : 597-608.
Perrier E, Bird N and  Rieu M, 1999. Generalizing the fractal model of soil structure: the pore-solid fractal approach. Geoderma 88: 137-164.
Perrier E, Bird N, 2002. Modelling soil fragmentation: the PSF approach. Soil Till Res 64: 91–99.
Perrier EMA, Bird NRA, 2003. The PSF model of soil structure: a multiscale approach. Pp. 1–18. In: Pachepsky Ya, Radcliffe DE, Selim HM (eds). Scaling Methods in Soil Physics. CRC Press, Boca Raton, FL,
Perfect E, Rasiah V and  Kay BD, 1992. Fractal dimensions of soil aggregate-size distributions  calculated by number and mass. Soil Sci Soc  Am  J 56: 1407-1409.
Pidgeon JD and Soane  BD, 1977. Effects of tillage and direct drilling on soil properties during the growing season in a long-term barley mono-culture system. J Agric Sci 88: 431-442.
Puppala AJ, Acar  YB and Tumay MT, 1995. Cone penetration in very weakly cemented sand. J Geotech  Eng  121: 589-600.
Schoeneberger PJ, Wysocki DA, Benham EC and Brooderson DE, 1998. Field book for describing and sampling soils. Natural Resourses Conservation Service, USDA, Natural Soil Survey Center, Lincoln, NE.
To J and Kay  BD, 2005. Variation in penetrometer resistance with soil properties: the contribution of effective stress and implications for pedotransfer functions. Geoderma 126: 261-276.
Tyler SW and Wheatcraft SW, 1990. Fractal process in soil water retention. Water Resour Res 26: 1047-1054.
Tyler SW and Wheatcraft SW, 1992. Fractal scaling of soil particle-size distribution: analysis and limitations. Soil Sci  Soc  Am  J  56: 362-369.
Upadhyaya SK, Kemble  LJ and  Collins  NE, 1982. Cone index prediction equations for Delaware soils. ASAE paper. 82: 1452-1456.
Usowicz B, and Lipiec J, 2009. Spatial distribution of soil penetration resistance as affected by soil compaction: The fractal approach. Ecological Complexity 6: 263–271.
Vaz CMP, Luis HB and  Hopmans  JW, 2001. Contribution of water content and bulk density to field soil penetration resistance as measured by a combined cone penetrometer-TDR probe. Soil Till Res 60: 35-42.
Wosten JHM, Pachepsky Y and Rawls WJ, 2001. Pedotransfer functions: bridging the gap between available basic soil data and missing soil hydraulic characteristics. J Hydrol 251: 123-150.