Microsoft Word 00 a loinoidau TV (moi thang1 2016) docx ISSN 1859 1531 TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, SỐ 2(111) 2017 Quyển 1 119 HIỆU QUẢ KỸ THUẬT VÀ CÁC NHÂN TỐ ẢNH HƯỞNG ĐẾN HIỆU QUẢ[.]
ISSN 1859-1531 - TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, SỐ 2(111).2017-Quyển 119 HIỆU QUẢ KỸ THUẬT VÀ CÁC NHÂN TỐ ẢNH HƯỞNG ĐẾN HIỆU QUẢ KỸ THUẬT CỦA CÁC HỘ TRỒNG LÚA Ở THỊ XÃ HƯƠNG TRÀ, TỈNH THỪA THIÊN HUẾ TECHNICAL EFFICIENCY AND ITS DETERMINANTS IN RICE FARMING IN HUONG TRA TOWN, THUA THIEN HUE PROVINCE Trần Hạnh Lợi Trường Đại học Kinh tế, Đại học Huế; tranhanhloi@gmail.com Tóm tắt - Sử dụng Hàm sản xuất biên ngẫu nhiên phầm mềm Frontier 4.1, nghiên cứu nhằm xác định hiệu kỹ thuật yếu tố định đến suất hiệu kỹ thuật hộ trồng lúa thị xã Hương Trà, tỉnh Thừa Thiên Huế Nghiên cứu sử dụng số liệu bảng hai vụ lúa năm 2015 120 hộ lựa chọn phương pháp chọn mẫu ngẫu nhiên phân tầng Hiệu kỹ thuật hộgiao động từ 0,66 đến 1, trung bình đạt 0,795, khơng có khác biệt suất hiệu kỹ thuật trung bình vụ lúa Lượng phân kali, công lao động nước có tác động tích cực đến suất lúa lượng lúa giống, thuốc BVTV có ảnh hưởng trái chiều Kết phân tích hàm phi hiệu kỹ thuật cho thấy yếu tố vị trí mảnh ruộng so với hệ thống thuỷ lợi, số lao động nông nghiệp việc cung cấp nước lịch yếu tố có tác động tích cực tới hiệu kỹ thuật trồng lúa Abstract - Using Stochastic Production Frontier and Frontier 4.1 program, this study estimates technical efficiency and determinants of rice yield and determinants of technical efficiency of rice-farming households in Huong Tra town, Thua Thien Hue province The study uses panel data in two rice growing seasons of 2015, collected from 120 rice-farming households which are selected by stratified random sampling method Technical efficiency of surveyed households ranges from 0.66 to and attains 0.795 on average There are no statistically significant differences between means of rice yield and technical efficiency in two rice growing seasons Rice yield is positively impacted by kalium fertilizer, labour and water whereas seed and pesticide have reverse effects on the output of production The result of technical inefficiency model shows that tail-end location of the plot, the number of farming labor, and in-time irrigation service are the major factors that cause positive effects on households’ technical efficiency Từ khóa - hiệu kỹ thuật; lúa; hàm sản xuất biên ngẫu nhiên; Hương Trà; thuỷ lợi; frontier 4.1 Key words - technical efficiency, rice, stochastic production frontier; Huong Tra, irrigation; frontier 4.1 Đặt vấn đề Nông nghiệp lĩnh vực quan trọng Việt Nam, đóng góp đến 17% giá trị quốc nội tạo công ăn việc làm cho 44% lực lượng lao động Trong nông nghiệp, lúa trồng chính, chiếm đến 52,5% diện tích đất trồng trọt đóng góp 89,5% tổng sản lượng lương thực có hạt [1] Việt Nam trở thành nước xuất lúa gạo lớn kể từ năm 1989, nước xuất lúa gạo lớn thứ hai giới, sau Thái Lan Năm 2015, Việt Nam xuất 6,6 triệu gạo, chiếm 15% tổng lượng gạo giao dịch giới [2] Có kết nhờ thay đổi chế quản lý đất đai, chế thị trường, áp dụng kỹ thuật sản xuất lúa đại, sử dụng phân bón, giống lúa suất cao Thị xã Hương Trà nằm trung độ tỉnh Thừa Thiên Huế, khu vực từ lâu gắn liền với hoạt động sản xuất lúa tỉnh Thừa Thiên Huế Với diện tích 6000ha,đây vựa lúa quan trọng góp phần khơng nhỏ cho sản lượng lúa toàn tỉnh, đồng thời đem lại nguồn thu nhập đáng kể cho người nơng dân.Tuy nhiên, diện tích, suất, sản lượng lúa địa bàn thị xã năm gần có khuynh hướng giảm với cấu giống lúa gieo trồng nghèo nàn, chủ yếu Khang Dân 18 (KD) Hương Thơm số (HT1)[3] Do câu hỏi đặt liệu hộ trồng lúa sử dụng yếu tố đầu vào hiệu chưa để đạt suất tối đa? Có yêu tố ảnh hưởng đến suất hiệu kỹ thuật hộ trồng lúa? Hiện có hai cách để tiếp cận đo lường hiệu kỹ thuật: phương pháp tham số phương pháp phi tham số Trong phương pháp tham số dựa vào dạng hàm cụ thể phương pháp phi tham số không xây dựng dạng hàm cụ thể giả định số liệu đưa vào phân tích khơng chứa sai số ngẫu nhiên đo lường Đối với phương pháp tham số có mơ hình thường sử dụng hàm sản xuất biên xác định (Deterministic Production Frontier - DPF) hàm sản xuất biên ngẫu nhiên (Stochastic Production Frontier - SPF) Trong mơ hình DPF, sai lệch khỏi giới hạn sản xuất xem phi hiệu Phương pháp bình phương bé có hiệu chỉnh (COSL) sử dụng cho kiểu hàm Trong đó, mơ hình SPF cho phép sai số ngẫu nhiên xung quanh hàm sản xuất ước lượng Phương pháp hàm sản xuất biên ngẫu nhiên cho phép phân tách ảnh hưởng nhiễu ảnh hưởng phi hiệu [4] Do hoạt động sản xuất lúa nước ta ví dụ điển hình hoạt động sản xuất sử dụng nhiều đầu vào để tạo đầu nhất, khả có sai số đo lường cao việc ghi chép thơng tin nơng hộ cịn hạn chế nên sử dụng cách tiếp cận tham số, cụ thể hàm sản xuất biên ngẫu nhiên, để đo lường hiệu kỹ thuật hoạt động trồng lúa nông hộ phù hợp Bên cạnh đó, nghiên cứu hiệu kỹ thuật trồng lúa phương pháp SPF nước ta chưa nhiều, tập trung vùng đồng Sông Cửu Long đồng Sông Hồng Các nghiên cứu trước thường tập trung vào yếu tố đầu vào truyền thống giống, phân bón, thuốc BVTV, công lao động [5,6,7, 8] mà chưa đưa vào yếu tố đầu vào quan trọng thiếu cho hoạt động trồng lúa, thuỷ lợi Xuất phát từ thực tế trên, nghiên cứu sử dụng phương pháp SPF để tiến hành xác định hiệu kỹ thuật 120 Trần Hạnh Lợi trồng lúa phân tích yếu tố tác động đến hiệu kỹ thuật, có yếu tố thuỷ lợi Số liệu phân tích phần mềm Frontier 4.1 cho phép chạy kết bước Phương pháp nghiên cứu Farrel (1957) người đề cập đến khái niệm hiệu kinh tế, hiệu phân phối hiệu kỹ thuật Ông giải thích hiệu kỹ thuật khả đạt đến mức sản lượng tối đa từ tập hợp định yếu tố đầu vào cho trước Nghiên cứu Farrell áp dụng rộng rãi có nhiều nghiên cứu cải biến phát triển lên từ nghiên cứu Một số mơ hình hàm sản xuất biên ngẫu nhiên (SPF) cho phép xác định hiệu kỹ thuật hiệu kinh tế đơn vị sản xuất ước lượng hợp lý cực đại (Maximun Likelihood Estimation – MLE) Hàm SPF đề xuất Aigner, Lovell Schmidt [9], Meeusen va Broceck [10], sau có nhiều nghiên cứu cải biến mơ hình hàm Nghiên cứu sử dụng mơ hình Battese Coelli [11] Sử dụng liệu bảng, hàm sản xuất biên ngẫu nhiên có dạng sau: ( − ) = 1, , ; = 1, … , (1) = ( ; ) Mơ hình (1) viết lại dạng mơ hình CobbDouglas Translog Nghiên cứu sử dụng mơ hình Cobb-Douglas mơ hình sử dụng rộng rãi xác định hiệu sản xuất nông hộ, theo Kopp Smith [12] mơ hình có ảnh hưởng đến việc đo lường hiệu thực nghiệm Mơ hình (1) viết lại sau: dạng = + = ÷ ( + ) Hệ số γcó giá trị từ đến 1, thể đóng góp tương đối Vit Uit đến tổng sai số mơ hình (2), sử dụng để kiểm định phù hợp hàm SPF so với hàm OLS Nếu γ tiến đến nên dùng hàm sản xuất OLS, cịn γ tiến đến nên dùng hàm SPF Dựa vào mơ hình lý thuyết nêu trên, nghiên cứu thực nghiệm sử dụng biến độc lập cho mơ hình hàm sản xuất (2) 12 biến độc lập cho mơ hình hàm phi hiệu kỹ thuật (3) Thông tin biến sử dụng mơ hình (2) (3) mô tả bảng Nguồn số liệu sơ cấp lấy từ điều tra hộ thị xã Hương Trà Phường lựa chọn để điều tra dựa quy mơ diện tích suất lúa để đảm bảo tính đại diện cho tồn thị xã Đồng thời, để thấy ảnh hưởng thuỷ lợi đến hiệu kỹ thuật trồng lúa, nghiên cứu lựa chọn địa điểm điều tra dựa vị trí vùng canh tác lúa so với hệ thống kênh thuỷ lợi (đầu kênh, kênh cuối kênh) Thị xã Hương Trà có hệ thống thủy lợi Khe Ngang lớn Ba phường - Hương Hồ, Hương An, Hương Chữ - nằm vị trí đầu, cuối hệ thống thuỷ lợi nên lựa chọn để điều tra Phương pháp chọn mẫu ngẫu nhiên phân tầng Tổng số mẫu điều tra 120 hộ Bảng Các biến Mô tả đặc điểm biến Tên biến Mơ tả biến A Mơ hình hàm sản xuất biên ngẫu nghiên Năng suất lúa (Yit) Năng suất lúa (tạ/sào/vụ) t = : vụ Đông Xuân, t = : vụ Hè Thu Giống (X1it) Lượng giống lúa (kg/sào/vụ) Trong đó, Yit suất hộ i thời điểm t; Xjit yếu tố đầu vào j hộ i thời điểm t; T thời điểm (năm/vụ…); β 0, βj vàβ k+1 hệ số cần ước lượng; Vit sai số ngẫu nhiên giả định có phân phối chuẩn đối xứng N(0,σv2) độc lập với Uit; Uit sai số tác động phi hiệu kỹ thuật, giả định có phân phối chuẩn cụt với giá trị trung bình μ phương sai σu2,và Uit≥ Uit xác định dạnghàm, gọi hàm phi hiệu kỹ thuật Hàm sử dụng để giải thích yếu tố ảnh hưởng đến phi hiệu kỹ thuật mô tả sau: = +∑ + (3) Phân lân (X2it) Lượng phân đạm (kg/sào/vụ) Phân đạm (X3it) Lượng phân lân (kg/sào/vụ) Phân kali (X4it) Lượng phân kali (kg/sào/vụ) Thuốc BVTV (X5it) Lượng thuốc BVTV (1000đ/sào/vụ) Nước (X6it) Số lần tưới nước (lần/vụ) Công lao động (X7it) Số công lao động (công/sào/vụ) Vụ canh tác (T) T=1 vụ Đông Xuân, T=0 vụ Hè Thu Trong đó, Zjit yếu tố ảnh hưởng đến phi hiêu kỹ thuật; δ0 δj hệ số cần ước lượng; ωit sai số ngẫu nhiên Hiệu kỹ thuật (TE) tỷ số suất thực tếvà suất tối đa đạt đơn vị sản xuất thứ i thời điểm t, có giá trị từ TE tính sau: TEit = E(Yit|Uit>0,Xit)/E(Yit|Uit=0,Xit) = exp (-Uit ) (4) Phương pháp MLE sử dụng để xác định đồng thời tham số mơ hình (2) mơ hình (3) chương trình Frontier 4.1 Coelli (2007) Nó ước lượng tham số phương sai hàm khả (likelihood function) = + + + − (2) B Mơ hình hàm phi hiệu kỹ thuật Tuổi (Z1it) Tuổi chủ hộ (tuổi) Trình độ (Z2it) Trình độ học vấn chủ hộ (số năm học) Loại giống (Z3it) = HT1 ; = Khang Dân Số LĐNN (Z4it) Số lao động nông nghiệp (người) Dịch bệnh lúa (Z5it) = Có ; = Khơng Tập huấn (Z6it) = Có ; = Không Thu nhập từ lúa (Z7it) Thu nhập từ trồng lúa hộ (ngàn đồng/năm) Ruộng đầu kênh (Z8it) = Có ; = Khơng Ruộng cuối kênh (Z9it) = Có ; = Khơng Thuỷ lợi cung cấp đủ nước = khơng đủ, = bình thường, = đầy đủ (Z10it) ISSN 1859-1531 - TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, SỐ 2(111).2017-Quyển Thuỷ lợi cung cấp nước = khơng lịch, = bình lịch (Z11it) thường, = lịch Thuỷ lợi cấp nước linh = khơng linh động, = bình thường, = linh động động cần (Z12it) Ghi : sào = 500 m Kết thảo luận 3.1 Đặc điểm, đầu vào kết hoạt động sản xuất lúa hộ điều tra Số liệu thống kê mô tả biến quan trọng sử dụng hàm sản xuất (2) hàm phi hiệu kỹ thuật (3) mô tả bảng Năng suất lúa bình quân hộ đạt 2,8 tạ/sào, thấp so với năm suất lúa bình quân nước tỉnh năm 2015 [1], [13] Lượng giống lúa gieo sào 5,4 kg, cao mức khuyến cáo đến kg/sào Lượng phân bón tỷ trọng phân Đạm – Lân – Kali tương đối hợp lý Điều phù hợp với kết nghiên cứu Lê Khắc Phúc cộng [3] Hoạt động trồng lúa mang lại thu nhập bình quân hộ đạt 14,4 triệu đồng/năm, giao động khoảng rộng, từ 2,5 triệu đến 37,3 triệu đồng/năm Đa phần chủ hộ tuổi trung niên, trình độ học vấn khơng cao, bình quân có 6,6 năm đến trường Trong năm 2015, 1/5 số hố hộ tập huấn kỹ thuật 1/3 số hộ có ruộng lúa bị sâu bệnh Tỷ trọng hộ gieo trồng giống lúa KD HT1 tương đối cân Phần lớn hộ đánh giá dịch vụ thuỷ lợi tốt thông qua tiêu chí cung cấp đủ nước, cung cấp lịch linh động cần 121 Theo thang đo Thuỷ lợi cung cấp đủ nước (Z10it) 0,0 33,8 66,2 Thuỷ lợi cung cấp nước lịch (Z11it) 4,2 34,2 61,6 Thuỷ lợi cấpnước linh động cần (Z12it) 3,3 27,9 68,8 Năng suất lúa (Yit) 2,81 3,5 2,2 0,26 Giống (X1it) 5,42 0,54 3.2 Các nhân tố ảnh hưởng đến suất lúa Kết trình bày bảng cho thấy ảnh hưởng yếu tố đầu vào đến suất lúa hộ nông dân Công lao động có tác độngmạnh đến suất lúa với mức ý nghĩa 1% Khi tăng 1% cơng lao động suất lúa tăng thêm 0,211% Do đó, để tăng suất lúa, hộ nông dân nên đầu tư thêm công lao động Tiếp theo yếu tố nước, cụ thể nêu số lần tưới nước tăng lên 1% suất lúa tăng lên 0,123% Trong đó, biến thuốc bảo vệ thực vật lại có tác động tiêu cực đến suất lúa với mức ý nghĩa 1% Khi tăng thêm 1% chi phí cho thuốc bảo vệ thực vật suất lúa giảm xuống 0,089% Điều cho thấy hộ sử dụng thuốc BVTV nhiều mức tối ưu Trong ảnh hưởng phân kali mang dấu dương có ý nghĩa thống kê mức 1% Hệ số cho biết, tăng lượng phân kali lên 1% suất lúa tăng lên 0,037% Các biến giống, phân đạm phân lân khơng có ý nghĩa thống kê, hệ số ảnh hưởng biến giống phân lân mang dấu âm điều cho thấy hộ nông dân đầu tư lượng giống phân lân nhiều mức cần thiết Điều phù hợp với nhận định ban đầu dựa số liệu thống kê mô tả Hệ số biến giả T (vụ canh tác) khơng có ý nghĩa thống kê mức 5% Điều có nghĩa sai khác suất lúa vụ năm 2015 khơng có ý nghĩa thống kê Phân lân (X2it) 7,55 13 2,59 Bảng Kết ước lượng hàm sản xuất Phân đạm (X3it) 10,41 20 3,72 Phân kali (X4it) 5,34 11 2,25 1,89 91,89 150 60 20,27 7,4 52,72 Trình độ (Z2it) Số LĐNN (Z4it) Bảng Số liệu thống kê hộ trồng lúa Biến định lượng Thuốc BVTV (X5it) Nước (X6it) Công lao động (X7it) Tuổi (Z1it) Thu nhập từ lúa (Z7it) Bình quân Tối đa Tối thiểu Sai số chuẩn Ký hiệu biến Tên biến Hệ số Độ lệch Giá trị t β0 Hằng số 1,456*** 0,197 7,404 1,42 Ln X1 Giống -0,065ns 0,042 -1,553 0,98 Ln X2 Phân lân -0,010ns 0,031 -0,310 78 33 10,28 Ln X3 Phân đạm 0,032ns 0,025 1,277 6,61 12 3,77 1,88 0,38 Ln X4 Phân kali 0,037* 0,019 1,908 14363,7 37315 2499 7819,5 Ln X5 Thuốc BVTV -0,089** 0,032 -2,842 Tỷ trọng (%) Ln X6 Số lần nước tưới 0,123*** 0,009 13,030 Theo thang đo Ln X7 Công lao động 0,211*** 0,022 9,708 T Vụ canh tác 0,028ns 0,015 0,169 0,999*** 0,013 76,680 Biến định tính Loại giống (Z3it) 51,2 48,8 Dịch bệnh lúa (Z5it) 62,5 37,5 Tập huấn (Z6it) 71,7 21,3 Ruộng đầu kênh (Z8it) 66,7 33,3 Ruộng cuối kênh (Z9it) 66,7 33,3 γ (Nguồn: ước lượng từ số liệu điều tra) (Ghi chú: *, **, *** mức ý nghĩa hệ số hồi quy tương ứng với 10%, 5%, 1%; ns hệ số hồi quy khơng có ý nghĩa thống kê) 122 Trần Hạnh Lợi Hệ số gama (γ) 0,99 (gần tiến tới 1)và có ý nghĩa thống kê mức 1% Điều cho thấy yếu tố phi hiệu kỹ thuật có ảnh hưởng đáng kể đến kết hoạt động sản xuất Hay nói cách khác hoạt động sản suất hộ không ảnh hưởng bới việc sử dụng yếu tố đầu vào mà bị ảnh hưởng yếu tố kinh tế - xã hội hay gọi yếu tố phi hiệu kỹ thuật Như vậy, việc sử dụng mơ hình hàm sản xuất biên ngẫu nhiên hoàn toàn phù hợp với số liệu nghiên cứu 3.3 Hiệu kỹ thuật hộ trồng lúa Để thể thay đổi phân phối hiệu kỹ thuật hộ vụ lúa năm 2015, thực chạy tần suất hiệu kỹ thuật thể bảng Hiệu kỹ thuật bình qn đạt 0,79 vụ Đơng Xn 0,8 vụ Hè Thu Kết kiểm định cho thấy khơng có khác biệt giá trị TE mùa vụ với mức ý nghĩa 5% Đa số hộ đạt hiệu kỹ thuật 0,7, 50% số hộ có hiệu kỹ thuật từ 0,7 đến 0,8 Tuy nhiên số hộ nông dân đạt hiệu từ 0,9 trở lên ít, chiếm 7,5% 9,2% tổng số hộ vụ Đông Xuân vụ Hè Thu Bảng Hiệu kỹ thuật nông hộ Phân tổ theo TE Đông Xuân So sánh ĐX/HT % Số hộ % Số hộ % 3,3 3,3 0,0 0,7– 0,8 69 57,5 65 54,2 -4 -3,3 0,8– 0,9 38 31,7 40 33,3 1,7 Ký hiệu biến 0,9 – 7,5 11 9,2 1,7 Z0