Giới thiệu dự án

Thị trường dược liệu nội địa Việt Nam tiêu thụ khoảng 60.000 – 80.000 tấn/năm nhằm phục vụ sản xuất dược phẩm và thực phẩm bảo vệ sức khỏe. Trong đó, hạt cây bìm bìm biếc (Pharbitis nil (L.) Choisy, hay còn gọi là Khiên ngưu tử, Hắc/Bạch sửu) chứa hoạt chất Pharbitin (~2%) cùng các dẫn xuất Acid pharbitic C, D có hoạt tính giải độc gan, nhuận tràng, thông tiểu và kháng khuẩn mạnh. Hiện nay, nguồn dược liệu này tại Việt Nam chủ yếu thu hái tự nhiên hoặc nhập khẩu không chính ngạch từ Trung Quốc, dẫn đến biến động lớn về chất lượng và hàm lượng hoạt chất.

Đề tài “Nghiên cứu ảnh hưởng của liều lượng và thời điểm bón phân đến cây Bìm bìm trong vụ xuân năm 2021 tại Gia Lâm – Hà Nội” (tác giả Chu Thị Đào, hướng dẫn bởi ThS. Nguyễn Thị Thanh Hải – Bộ môn Cây Công nghiệp & Cây thuốc, Khoa Nông học, Học viện Nông nghiệp Việt Nam) giải quyết bài toán tiêu chuẩn hóa quy trình thâm canh nông học cho cây dược liệu Pharbitis nil.

Mục tiêu dự án

  1. Đánh giá ảnh hưởng của các công thức phân bón (NPK cơ bản, NPK bổ sung vi lượng YaraMila WINNER, NPK bổ sung chế phẩm nấm rễ nội cộng sinh AMF) và thời điểm bón thúc (1 tháng vs 2 tháng sau gieo) đến sinh trưởng, phát triển của cây bìm bìm biếc.
  2. Xác định biến thiên chỉ số sinh lý ($SPAD$, $F_v/F_m$, diện tích lá $LAI$, tích lũy chất khô) và mức độ nhiễm sâu bệnh hại theo tiêu chuẩn QCVN 01-38:2010/BNNPTNT.
  3. Xác định công thức phân bón và thời điểm bón thúc tối ưu cho năng suất thực thu ($NSTT$) cao nhất.

Phạm vi và Giới hạn

  • Phạm vi nghiên cứu: Giống bìm bìm biếc PN2 (lá xẻ thùy, hạt đen), thực hiện trong vụ Xuân 2021 (tháng 1 – tháng 9/2021) trên đất phù sa cổ tại Gia Lâm, Hà Nội.
  • Giới hạn: Nghiên cứu tập trung vào năng suất nông học và chỉ số sinh lý; chưa mở rộng đánh giá định lượng hoạt chất Pharbitin bằng phương pháp HPLC trên các lô thí nghiệm khác nhau.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Phương pháp canh tác Ưu điểm Nhược điểm Đánh giá khả thi
Thu hái tự nhiên Chi phí đầu tư ban đầu bằng 0 Sản lượng phân tán, chất lượng không đồng nhất, nguy cơ cạn kiệt nguồn gen Không đáp ứng tiêu chuẩn GACP-WHO
Bón vô cơ đơn thuần (Đỗ Ngọc Hoán, 2016) Năng suất hạt đạt 9,43 tạ/ha với công thức 60N : 120P₂O₅ : 60K₂O Đất dễ thoái hóa, chưa khai thác được hệ vi sinh vật bản địa Khả thi trung bình, hiệu suất thấp
Bón kết hợp NPK + Vi lượng / AMF (Đề tài đề xuất) Năng suất đạt 11,87 tạ/ha (+25,8%), rễ phát triển mạnh, cây chống chịu tốt Đòi hỏi kỹ thuật bón đúng thời điểm (sau gieo 1 tháng) Tối ưu cho quy mô công nghiệp

Yêu cầu hệ thống canh tác theo MoSCoW

  • Must have: Nền phân bón cân đối 90N : 120P₂O₅ : 90K₂O; bón thúc đúng pha sinh trưởng sau gieo 30 ngày (T1); mật độ trồng 150.000 cây/ha.
  • Should have: Bổ sung phân bón phức hợp chứa trung vi lượng (Mg, S, B, Mn, Zn) hoặc nấm rễ nội cộng sinh AMF (Glomus spp.).
  • Could have: Tích hợp hệ thống tưới nhỏ giọt kết hợp châm phân (Fertigation).
  • Won't have: Bón thừa đạm tự do (>120 kg N/ha) gây lốp đổ và ức chế phân hóa mầm hoa do dư ion $NH_4^+$.

Thiết kế hệ thống thí nghiệm (Split-plot Design)

Thí nghiệm 2 nhân tố được bố trí theo kiểu ô chính - ô phụ (Split-plot) với 3 lần nhắc lại:

  • Ô chính (Nhân tố A - Thời điểm bón thúc):
    • $T_1$: Bón thúc 50% N + 50% K₂O sau khi gieo 1 tháng.
    • $T_2$: Bón thúc 50% N + 50% K₂O sau khi gieo 2 tháng.
  • Ô phụ (Nhân tố B - Liều lượng và loại phân bón):
    • $L_1$: Nền NPK tiêu chuẩn (90N : 120P₂O₅ : 90K₂O).
    • $L_2$: Nền NPK + Phân bón YaraMila WINNER (15-9-20 + 1,8% MgO + 3,8% S + TE: B, Mn, Zn).
    • $L_3$: Nền NPK + Chế phẩm nấm rễ AMF (Glomus spp., Bacillus sp., Azotobacter sp., Burkholderia sp.).

Công cụ và Công nghệ phân tích (Technology Stack)

  • Phần mềm thống kê nông học: IRRISTAT version 5.0 (International Rice Research Institute) phục vụ phân tích ANOVA và Duncan test.
  • Phần mềm xử lý dữ liệu: Microsoft Excel 365 (Data Analysis Toolpak).
  • Thiết bị đo quang sinh học: Máy đo chỉ số diệp lục Minolta SPAD-502Plus và máy đo huỳnh quang diệp lục Handy PEA ($F_v/F_m$).
  • Tiêu chuẩn kiểm định sâu bệnh: Quy chuẩn kỹ thuật quốc gia QCVN 01-38:2010/BNNPTNT.

Implementation và kết quả

Quy trình canh tác và Kỹ thuật chi tiết

Quy trình canh tác chuẩn hóa:
1. Xử lý đất: Cày bừa, lên luống cao 30cm, rãnh rộng 30cm, xử lý nấm đất.
2. Xử lý hạt giống PN2: Chà cát nhẹ -> Ngâm nước ấm 45-50°C (2-3 giờ) -> Loại hạt lép.
3. Gieo trồng: Mật độ 150.000 cây/ha (hốc cách hốc 20cm, hàng cách hàng 40cm, 1-2 hạt/hốc).
4. Bón phân:
   - Bón lót: 100% P₂O₅ (Supe lân 16%) + 50% N (Ure 46%) + 50% K₂O (Kaliclorua 61%) + 2 tấn phân vi sinh Sông Gianh.
   - Bón thúc (T1: sau 30 ngày / T2: sau 60 ngày): 50% N còn lại + 50% K₂O còn lại (+ YaraMila / AMF theo lô).
5. Làm giàn: Cắm cọc chữ A hoặc hàng rào cao 2,0 - 2,5m tại thời điểm 20 ngày sau mọc.
6. Phòng trừ sâu bệnh: Pegasus 500SC (7-10ml/10L) trừ rệp; Accentan 50EC (15ml/10L) trừ bọ xít.
7. Thu hoạch: Sau 70-90 ngày khi quả chín 90%, phơi/sấy ở 65-70°C đến khi độ ẩm < 12%.

Mã nguồn thuật toán phân tích số liệu ANOVA Split-Plot

Dưới đây là mô hình xử lý thống kê được xây dựng tương thích với thuật toán của IRRISTAT 5.0 và biểu diễn qua R/Python:

import numpy as np
import pandas as pd
import statsmodels.api as sm
from statsmodels.formula.api import ols

def analyze_split_plot_agronomy(data_file: str):
    """
    Phân tích ANOVA mô hình Split-plot 2 nhân tố cho cây bìm bìm biếc
    Nhân tố chính: Main_Plot (Thời điểm bón T1, T2)
    Nhân tố phụ: Sub_Plot (Liều lượng L1, L2, L3)
    Khối: Block (Rep I, II, III)
    """
    df = pd.read_csv(data_file)
    
    # Mô hình Split-Plot ANOVA
    # Error A = Block:Main_Plot
    # Error B = Residual
    model = ols('Yield ~ C(Block) + C(Main_Plot) + C(Block):C(Main_Plot) + '
                'C(Sub_Plot) + C(Main_Plot):C(Sub_Plot)', data=df).fit()
    
    anova_table = sm.stats.anova_lm(model, typ=2)
    
    # Tính Năng suất lý thuyết (NSLT) và Năng suất thực thu (NSTT)
    # NSLT (tạ/ha) = (Số quả chắc/cây * Số hạt/quả * P1000 hạt * Mật độ) / 100.000
    df['Calculated_NSLT'] = (df['Pods_per_Plant'] * df['Seeds_per_Pod'] * 
                             df['P1000_weight_g'] * df['Plant_Density']) / 100000.0
    
    return anova_table, df.groupby(['Main_Plot', 'Sub_Plot'])['Yield'].mean()

# Thuật toán tính chỉ số hiệu suất huỳnh quang diệp lục:
def calculate_photosynthetic_efficiency(fo: float, fm: float) -> float:
    """Tính tỷ số Fv/Fm đo hiệu suất lượng tử cực đại của PSII"""
    fv = fm - fo
    return fv / fm

Kết quả thực nghiệm và Đo lường sinh trắc

Bảng tổng hợp chỉ tiêu sinh trưởng và năng suất giữa các công thức:
========================================================================================
Công thức   Chiều cao (cm)   Số lá/cây   SPAD    Fv/Fm    Số quả/cây   NSTT (tạ/ha)
----------------------------------------------------------------------------------------
L1T1        185,4            42,3        41,2    0,782    98,4         10,85
L1T2        172,1            38,6        39,5    0,765    91,2         10,26  (Thấp nhất)
L2T1        204,6            48,7        44,8    0,812    118,6        11,87  (Cao nhất)
L2T2        189,2            43,1        42,0    0,790    104,3        10,98
L3T1        218,5            52,4        46,5    0,828    112,1        11,54
L3T2        196,3            45,8        43,1    0,798    101,5        10,72
========================================================================================
LSD (0.05)   8,4              2,1         1,2     0,015    5,2          0,42
CV (%)       4,8              3,9         2,5     1,4      4,1          3,2
  • Sinh trưởng sinh dưỡng: Công thức L3T1 (NPK + AMF, bón thúc sau 1 tháng) đạt giá trị cao nhất về chiều cao cây (218,5 cm), số lá (52,4 lá/cây), chiều dài rễ và khả năng tích lũy chất khô nhờ mạng lưới sợi nấm Glomus spp. tăng cường hấp thu ion $H_2PO_4^-$ và nước.
  • Năng suất kinh tế: Công thức L2T1 (NPK + YaraMila WINNER, bón thúc sau 1 tháng) cho năng suất thực thu cao nhất đạt 11,87 tạ/ha, vượt trội hơn 15,7% so với đối chứng L1T2 (10,26 tạ/ha) với mức ý nghĩa thống kê $P < 0,05$.

Đổi mới và đóng góp

Đổi mới kỹ thuật

  1. Tích hợp nấm rễ cộng sinh AMF (Arbuscular Mycorrhizal Fungi): Đề tài là một trong những công trình tiên phong tại Việt Nam ứng dụng tổ hợp nấm Glomus spp. và vi khuẩn đối kháng rễ (Bacillus sp., Burkholderia sp.) trên đối tượng cây dược liệu dây leo ngắn ngày họ Convolvulaceae, giúp tăng chỉ số $F_v/F_m$ lên 0,828 và tăng sinh khối rễ vượt trội.
  2. Chuẩn hóa thời điểm bón thúc T1 (30 ngày sau gieo): Chứng minh bằng thực nghiệm rằng bón thúc sớm tại pha bắt đầu leo giàn đáp ứng đúng nhu cầu đồng hóa dưỡng chất trước khi phân hóa mầm hoa (tránh ức chế sinh thực do bón muộn ở 60 ngày).
  3. Phối chế phân phức hợp trung vi lượng: Bổ sung Magie (1,8%) và Bo (0,015%) thông qua phân YaraMila WINNER làm tăng tỷ lệ đậu quả và số hạt chắc/quả (đạt 118,6 quả/cây ở L2T1).

So sánh với các nghiên cứu trước đây

Tiêu chí Đỗ Ngọc Hoán (2016) Dương Thị Duyên & cs. (2017) Nghiên cứu này (2021)
Công thức phân bón 60N : 120P₂O₅ : 60K₂O đơn thuần 45N : 90P₂O₅ : 90K₂O + 2 tấn PVS 90N : 120P₂O₅ : 90K₂O + YaraMila / AMF
Thời điểm tác động Bón định kỳ thông thường Bón lót + thúc truyền thống Phân hóa rõ rệt 2 mốc thúc T1 (30d) vs T2 (60d)
Năng suất thực thu 9,43 tạ/ha 10,15 tạ/ha 11,87 tạ/ha (Tăng 17 – 25,8%)
Đo lường quang sinh học Không đo Đo diện tích lá cơ bản Đo trực tiếp $SPAD$, $F_v/F_m$, chỉ số $LAI$

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng sản xuất (Real-world Use Cases)

  • Liên kết chuỗi cung ứng dược liệu chuẩn GACP-WHO: Cung cấp nguồn nguyên liệu hạt Khiên ngưu tử đạt chuẩn hàm lượng cho các đơn vị sản xuất thuốc đông dược quy mô lớn (điển hình như Công ty Cổ phần Traphaco sản xuất thuốc giải độc gan nhuận mật).
  • Mô hình thâm canh chuyển đổi đất lúa kém hiệu quả: Triển khai tại các vùng chuyên canh dược liệu ven sông Hồng (Gia Lâm, Hưng Yên, Bắc Ninh) cho thu hoạch sau 75-90 ngày.

Phân tích Hiệu quả Kinh tế & ROI (Tính trên 1 ha canh tác)

Bảng hạch toán kinh tế (Đơn vị tính: VNĐ / ha / vụ):
----------------------------------------------------------------------------------------
1. TỔNG CHI PHÍ ĐẦU TƯ:                                             45.800.000 đ
   - Giống PN2 (25 kg @ 200.000 đ/kg):                                5.000.000 đ
   - Phân bón vô cơ + Phân vi sinh Sông Gianh:                        12.500.000 đ
   - Phân bón YaraMila WINNER / Chế phẩm AMF:                         4.800.000 đ
   - Cọc nứa, lưới làm giàn (khấu hao 2 vụ):                          8.500.000 đ
   - Công lao động (Làm đất, gieo, chăm sóc, thu hái):               15.000.000 đ
----------------------------------------------------------------------------------------
2. DOANH THU DỰ KIẾN (Công thức L2T1):
   - Năng suất thực thu: 1.187 kg/ha (11,87 tạ/ha)
   - Giá bán bao tiêu hạt khô: 90.000 đ/kg
   - Tổng doanh thu: 1.187 kg * 90.000 đ/kg =                         106.830.000 đ
----------------------------------------------------------------------------------------
3. LỢI NHUẬN RÒNG & CHỈ SỐ HOÀN VỐN:
   - Lợi nhuận thuần (Lãi ròng): 106.830.000 - 45.800.000 =           61.030.000 đ/ha/vụ
   - Tỷ suất sinh lời (ROI): (61.030.000 / 45.800.000) * 100% =      133,25%
   - Thời gian hoàn vốn: 01 chu kỳ vụ (90 ngày)
----------------------------------------------------------------------------------------

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Tính đại diện thời gian: Dữ liệu thu thập trong vụ Xuân 2021; cần kiểm chứng lặp lại trong vụ Thu Đông để xác định chính xác ảnh hưởng của biên độ nhiệt độ và quang chu kỳ.
  • Chỉ tiêu hoạt chất sâu: Chưa tiến hành sắc ký lỏng hiệu năng cao (HPLC) để định lượng tỷ lệ dẫn xuất Acid pharbitic C/D tương ứng với từng mức phân bón.

Hướng nghiên cứu tiếp theo

  1. Thử nghiệm kết hợp bón phân qua rễ (L2T1) với phun bổ sung phân bón lá chứa Bo - Canxi ở giai đoạn hoa rộ (45 ngày sau gieo).
  2. Xây dựng quy trình chuẩn hóa chiết xuất hoạt chất Pharbitin đạt tiêu chuẩn Dược điển Việt Nam V từ nguồn nguyên liệu trồng tập trung.

Đối tượng hưởng lợi

  • Sinh viên & Học viên cao học Nông học/Dược học: Tiếp cận phương pháp thiết kế thí nghiệm Split-plot, cách sử dụng phần mềm IRRISTAT 5.0 và kỹ thuật đo đạc quang sinh học ($SPAD$, $F_v/F_m$).
  • Kỹ sư nông nghiệp & Trang trại dược liệu: Có quy trình canh tác chi tiết từ mật độ gieo (150.000 cây/ha), loại phân bón (YaraMila WINNER, AMF) đến lịch trình bón thúc chính xác sau 30 ngày.
  • Doanh nghiệp Dược phẩm (Traphaco, Dược liệu TW): Tự chủ nguồn nguyên liệu sạch, đồng nhất chất lượng, giảm 100% phụ thuộc vào nguồn nhập khẩu trôi nổi.
  • Nhà nghiên cứu Sinh lý thực vật: Cung cấp bộ dữ liệu thực nghiệm về động thái trao đổi chất và phản ứng quang hợp của Pharbitis nil dưới tác động dinh dưỡng đa - trung - vi lượng.

Câu hỏi thường gặp

1. Yêu cầu kỹ thuật đất đai và tiểu khí hậu để triển khai quy trình là gì?

Cây bìm bìm biếc yêu cầu đất tơi xốp, thoát nước tốt, thành phần cơ giới nhẹ đến trung bình (đất phù sa, đất thịt nhẹ), tầng canh tác dày >25 cm, pH từ 5,5 – 7,0. Cây ưa sáng hoàn toàn và cần nhiệt độ trung bình từ 22 – 32°C trong giai đoạn phát triển thân lá.

2. Tại sao thời điểm bón thúc sau gieo 1 tháng (T1) lại vượt trội so với sau 2 tháng (T2)?

Giai đoạn 30 ngày sau gieo là thời điểm cây bắt đầu vươn ngọn leo giàn và phân cành cấp 1, chuẩn bị phân hóa mầm hoa. Cung cấp đủ dinh dưỡng ở pha này giúp tối ưu hóa diện tích lá ($LAI$) và hàm lượng diệp lục ($SPAD$). Nếu bón muộn sau 60 ngày (T2), cây đã bước vào giai đoạn nuôi quả rộ, việc bón đạm muộn dễ kích thích sinh trưởng sinh dưỡng thứ cấp, gây rụng hoa và lép hạt.

3. Nấm rễ AMF mang lại lợi ích gì so với bón phân khoáng thông thường?

Nấm nội cộng sinh AMF (Glomus spp.) thiết lập mạng lưới sợi nấm ngoại sinh trong vùng rễ, làm tăng diện tích tiếp xúc của hệ rễ lên gấp hàng chục lần, giải phóng enzyme phosphatase phân giải lân khó tiêu và kích thích cây tiết phytoalexin chống chịu bệnh lở cổ rễ (Rhizoctonia solani).

4. Cách phòng trừ sâu bệnh hại chính trên cây bìm bìm biếc?

  • Rệp muội chích hút: Thăm đồng giai đoạn cây con, phun Pegasus 500SC (7-10ml/bình 10-12 lít) khi mật độ vượt ngưỡng gây hại.
  • Bọ xít dài hại quả non: Phun Accentan 50EC (15ml/bình 10-12 lít) vào giai đoạn hoa rộ; đảm bảo thời gian cách ly tối thiểu 14 ngày trước thu hoạch.

5. Chi phí đầu tư 1 ha và thời gian thu hồi vốn thực tế ra sao?

Tổng chi phí đầu tư khoảng 45,8 triệu VNĐ/ha/vụ (đã bao gồm giống, phân bón, cọc giàn và nhân công). Với năng suất 11,87 tạ/ha và giá bán 90.000 đ/kg hạt khô, lợi nhuận ròng đạt trên 61 triệu VNĐ/ha chỉ sau một chu kỳ canh tác 3 tháng (ROI đạt 133,25%).

Kết luận

Đề tài “Nghiên cứu ảnh hưởng của liều lượng và thời điểm bón phân đến cây Bìm bìm trong vụ xuân năm 2021 tại Gia Lâm – Hà Nội” đã giải quyết thành công bài toán thâm canh cây dược liệu Pharbitis nil. Kết quả nghiên cứu xác định công thức tối ưu L2T1 (90N : 120P₂O₅ : 90K₂O kết hợp phân bón phức hợp YaraMila WINNER, bón thúc sau gieo 1 tháng) mang lại năng suất thực thu cao nhất đạt 11,87 tạ/ha, cùng công thức L3T1 bổ sung nấm rễ sinh học AMF giúp tối ưu hóa hệ rễ và chỉ số sinh lý quang hợp. Đây là cơ sở khoa học và thực tiễn vững chắc để chuyển giao gói kỹ thuật canh tác bìm bìm biếc đạt chuẩn GACP-WHO, góp phần nâng cao giá trị chuỗi sản xuất dược liệu Việt Nam.